Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,474,138 (A)6→7 coding (15/594 nt) bdcR → transcriptional repressor for divergent bdcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,474,1321.A98.7% 244.9 / ‑3.8 76coding (9/594 nt)bdcRtranscriptional repressor for divergent bdcA
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (40/35);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

ACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGG  >  NC_000913/4474009‑4474263
                                                                                                                            |                                                                                                                                    
aCTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AAAGAAACAATCTcg                                                                                                                       >  2:1089956/1‑139 (MQ=255)
    ccACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGtt                                                                                                                   >  1:1421521/1‑139 (MQ=255)
          aGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTcg                                                                                                                       >  2:673193/1‑129 (MQ=255)
          aGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTcg                                                                                                                       <  1:673193/129‑1 (MQ=255)
           ggATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGgtc                                                                                                            >  1:393036/1‑139 (MQ=255)
                       tCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGtt                                                                                                >  1:1477472/1‑139 (MQ=255)
                          tACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCg                                                                                             <  1:865267/139‑1 (MQ=255)
                            ccTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCt                                                                                           >  1:1286835/1‑139 (MQ=255)
                               gTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCCTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTcg                                                                                                                       <  2:882834/108‑1 (MQ=255)
                               gTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTcg                                                                                                                       >  1:882834/1‑108 (MQ=255)
                                 aaaaGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGa                                                                                      >  2:1176663/1‑139 (MQ=255)
                                 aaaaGCGCCCATAAGGACTCCTTGATTTATTATGTAACAGGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGa                                                                                      <  2:338263/139‑1 (MQ=255)
                                   aaGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGc                                                                                    <  2:1421521/139‑1 (MQ=255)
                                      cgcCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTcg                                                                                                                       <  1:1614212/101‑1 (MQ=255)
                                      cgcCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTcg                                                                                                                       >  2:1614212/1‑101 (MQ=255)
                                         ccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                          <  2:1364476/127‑1 (MQ=255)
                                         ccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                          >  1:1364476/1‑127 (MQ=255)
                                                   tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGc                                                                    >  2:227587/1‑139 (MQ=255)
                                                                   tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAgg                                                                                                              >  2:322155/1‑81 (MQ=255)
                                                                   tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAgg                                                                                                              <  1:322155/81‑1 (MQ=255)
                                                                      cATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAgg                                                                                                              >  1:1414953/1‑78 (MQ=255)
                                                                      cATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAgg                                                                                                              <  2:1414953/78‑1 (MQ=255)
                                                                        tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTcg                                                                                                                       <  1:371602/67‑1 (MQ=255)
                                                                        tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTcg                                                                                                                       >  2:371602/1‑67 (MQ=255)
                                                                        tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCtttt                                                        <  1:1313288/130‑1 (MQ=255)
                                                                        tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCtttt                                                        >  2:1313288/1‑130 (MQ=255)
                                                                        tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                               >  2:1051170/1‑139 (MQ=255)
                                                                        tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                               >  1:1666198/1‑139 (MQ=255)
                                                                        tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                               >  2:1315600/1‑139 (MQ=255)
                                                                        tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                               >  1:1156794/1‑139 (MQ=255)
                                                                         gCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaaa                                              >  2:1598287/1‑139 (MQ=255)
                                                                           aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTcg                                                                                                                       <  1:778797/64‑1 (MQ=255)
                                                                           aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTcg                                                                                                                       >  2:778797/1‑64 (MQ=255)
                                                                           aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAgg                                            >  2:1538150/1‑139 (MQ=255)
                                                                                 aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAgg                                                                                                              <  2:630561/67‑1 (MQ=255)
                                                                                 aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAgg                                                                                                              >  1:630561/1‑67 (MQ=255)
                                                                                 aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGa                                      <  2:1286835/139‑1 (MQ=255)
                                                                                 aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGa                                      <  2:294869/139‑1 (MQ=255)
                                                                                 aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGa                                      <  1:1598287/139‑1 (MQ=255)
                                                                                 aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGa                                      <  1:1051170/139‑1 (MQ=255)
                                                                                  aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGAt                                     <  1:227587/139‑1 (MQ=255)
                                                                                    cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTcg                                                                                                                       <  1:341415/55‑1 (MQ=255)
                                                                                    cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTcg                                                                                                                       >  2:341415/1‑55 (MQ=255)
                                                                                     tGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCaaa                                                                 <  2:1364396/108‑1 (MQ=255)
                                                                                     tGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCaaa                                                                 >  1:1364396/1‑108 (MQ=255)
                                                                                       ttttAACTTTCTGTCGACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCGGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                                >  1:1115812/1‑139 (MQ=255)
                                                                                       ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCttttt                                                       >  2:692250/1‑116 (MQ=255)
                                                                                       ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCttttt                                                       <  1:692250/116‑1 (MQ=255)
                                                                                       ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGttt                                         <  1:1498786/130‑1 (MQ=255)
                                                                                       ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGttt                                         >  2:1498786/1‑130 (MQ=255)
                                                                                         ttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCttttt                                                       >  2:586630/1‑114 (MQ=255)
                                                                                         ttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCttttt                                                       <  1:586630/114‑1 (MQ=255)
                                                                                          tAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAg                             <  1:1538150/139‑1 (MQ=255)
                                                                                               ttCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTg                        >  2:936520/1‑139 (MQ=255)
                                                                                                tCTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGc                       <  2:1156794/139‑1 (MQ=255)
                                                                                                 cTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCt                                                                                           <  1:781348/70‑1 (MQ=255)
                                                                                                 cTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCt                                                                                           >  2:781348/1‑70 (MQ=255)
                                                                                                 cTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                     <  2:189148/76‑1 (MQ=255)
                                                                                                 cTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                     >  1:189148/1‑76 (MQ=255)
                                                                                                 cTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCt                      <  2:1531018/139‑1 (MQ=255)
                                                                                                 cTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCt                      >  1:1386402/1‑139 (MQ=255)
                                                                                                 cTGTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCt                      >  1:1531018/1‑139 (MQ=255)
                                                                                                   gTCAACAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGa                    >  1:748820/1‑139 (MQ=255)
                                                                                                        cAGGTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTa               <  1:936520/139‑1 (MQ=255)
                                                                                                           gTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCt                                                                                        <  1:534544/63‑1 (MQ=255)
                                                                                                           gTTTTGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCt                                                                                        >  2:534544/1‑63 (MQ=255)
                                                                                                            ttttGTAATGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa           <  2:393036/139‑1 (MQ=255)
                                                                                                                   aTGGTCACTAAAA‑AAATAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCt                                  >  1:27966/1‑109 (MQ=255)
                                                                                                                   aTGGTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCt                                  <  2:27966/109‑1 (MQ=255)
                                                                                                                     ggTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCt                                                                                        <  2:1214942/53‑1 (MQ=255)
                                                                                                                     ggTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCt                                                                                        >  1:1214942/1‑53 (MQ=255)
                                                                                                                     ggTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTgg  >  2:1136757/1‑139 (MQ=255)
                                                                                                                      gTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAg                             <  1:1227347/111‑1 (MQ=255)
                                                                                                                      gTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAg                             >  2:1227347/1‑111 (MQ=255)
                                                                                                                      gTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCt     <  2:1263971/135‑1 (MQ=255)
                                                                                                                      gTCACTAAAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCt     >  1:1263971/1‑135 (MQ=255)
                                                                                                                            |                                                                                                                                    
ACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AAA‑AAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGG  >  NC_000913/4474009‑4474263

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑