Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091329,4960TA43.4% 95.9 / 60.2 83intergenic (+301/‑155)dapB/carAdihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (39/8);  new base A (0/36);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 4.41e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29347‑29642
                                                                                                                                                        |                                                                                                                                                  
ttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAatat                                                                                                                                                    <  1:466255/150‑1 (MQ=255)
 tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                                   <  1:431582/150‑1 (MQ=255)
  aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                                 <  1:462129/151‑1 (MQ=255)
     gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATaaaaa                                                                                                                                               >  2:800600/1‑150 (MQ=255)
       aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATc                                                                                                                                            <  1:592883/151‑1 (MQ=255)
        aaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAAt                                                                                                                                             <  1:271255/149‑1 (MQ=255)
         aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                                          >  2:572552/1‑151 (MQ=255)
         aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                                          <  2:471574/151‑1 (MQ=255)
             tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                       >  2:80816/1‑150 (MQ=255)
              gACCATTTGGTCCACTTTTTTCTGCTCTTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGTTT‑ATATGATCTTTTTGTCCCTGAATGCC‑GGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCa                                                                                                                                      <  1:1092451/150‑1 (MQ=255)
                ccATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCAtt                                                                                                                                    >  2:977391/1‑150 (MQ=255)
                 cATTTGGGCCACTTTTTTCTGCCCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAAAACAGGTT‑AGATGATCTTTTTTTCGCTTAATGCC‑TGCAAAACTTGCATGAGCCACAAAATAATATAAAAAATCCCGCCa                                                                                                                                      <  1:600388/147‑1 (MQ=255)
                  aTTTGGTCCACTTTTTTATGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                   >  2:579875/1‑149 (MQ=255)
                     tGGTCCTTTTTTTTGTTTTCTTTTTTATTTCATGCAATCTTCTTGCTGCGCAGGGTTTTTC‑AAAAAAAGGTT‑ATATGATTTTTTTTTCGTTTAATCCC‑TAAAAAACAGGTATGACCCAAAAAAAAATAAAAAAATCCCCGCCATTAAGtt                                                                                                                               <  1:826078/150‑1 (MQ=255)
                       gTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGAc                                                                                                                            >  2:505746/1‑151 (MQ=255)
                         ccctttttttCTGCTCTTTTTTTTTTCATGCATTTTTCTTGCTGCGCAAGCTTTTTC‑CAAAACAGGTT‑AGATGATTTTTTTGTCGTTAAATCCC‑TGAAAAACATGCAGGACCCAAAAAATAATATAAAAAATCCCGCCATTAAGtt                                                                                                                               <  1:148416/142‑1 (MQ=255)
                            tttttttCTCCTCTTTTTTTTTTCAGGAATTTTTCTGGCGGCCAAAGGTTTTCC‑CAAAAAAGGTT‑AGATGATTTTTTTGTCGCTAAAGGCT‑TAAAAAACAGGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                         <  2:383856/148‑1 (MQ=255)
                             tttttttcccTCGTTTTTATTTCATGCATTCTTCTTGCGGCG‑AAGGTTTTTCAAAAAACGGGTA‑AGAGTATTTTTTTGACGTTTAATGCG‑TCTAAAACATCCATGACCCAAAAAATAAAAAAAAAATTCCCGCCATTAAGTTGACTTTTa                                                                                                                       <  1:201734/141‑1 (MQ=255)
                             ttttttCTGCTCTTTTTTATTTCATGCATTCTTCTTGCGGCGCAAGCGTTTCC‑AAAAACGGTTG‑AGATAATTTTTTTGTCGTTAAATGCC‑TAAAAAACATGCAGGACCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTa                                                                                                                       <  2:822192/150‑1 (MQ=255)
                             ttttttCTGCTCTTTTTTATTTCATGCAATCTTCTTGCTGCGCAGGCTTTTTC‑CAAAACAGGTT‑AGATGATTTTTTTGTCGCTTAATGCC‑TGAAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                        <  2:645916/149‑1 (MQ=255)
                              tttttCTTCTTTTTTTTTTTTCATGCAATCTTCTTGCTGGCCAGCCGTTTAC‑AAAAACAGGTT‑AAATAATCTTTTTGTCACTAATTGCC‑AGTAAAACATGCATCGGCCAAAAAATAATATAAAAACTCCCGCCATTAAGTTGACTTTTAg                                                                                                                      <  1:450310/150‑1 (MQ=255)
                                   cTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                                 >  1:1013627/1‑150 (MQ=255)
                                   cTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGAGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                                 >  1:279985/1‑150 (MQ=255)
                                      cttttttttATTTCATGCAATTTTTTTGCTGCGCAGGCGTTTCC‑CAAAAGAGTTT‑AGATGTTTTTTTTGTCGCTTAATCCC‑TATAAACCATGCAGGAGCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCa                                                                                                                <  2:811196/144‑1 (MQ=255)
                                       tttttttttttTCATGCATTCTTTTTGCTGCGCAAGCTTTTCC‑AAGAACAGTTT‑AGATGTTTTTTTTGCCGTTTAATGCC‑TGAAAAACATGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGcc                                                                                                                  <  2:621081/142‑1 (MQ=255)
                                        cGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                           >  2:1021362/1‑151 (MQ=255)
                                            tttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                         >  2:645828/1‑149 (MQ=255)
                                             ttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                         >  2:257185/1‑148 (MQ=255)
                                                 ttCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgc                                                                                                  >  2:329795/1‑151 (MQ=255)
                                                        aaTCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                           >  2:93910/1‑151 (MQ=255)
                                                          tttttttGCTGCCAAGGCTTTTTC‑AAGAACGGTTT‑AGATGTTTTTTTTTTCGTTTAATCCT‑TGAAAAAAATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                              <  2:175417/141‑1 (MQ=255)
                                                          tcttctTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTg                                                                                         >  1:1095997/1‑151 (MQ=255)
                                                            ttcttGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                       >  2:877301/1‑151 (MQ=255)
                                                              cttGCTGCGCAAGCGTTTGC‑AAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                       >  1:796909/1‑149 (MQ=255)
                                                                tgctgcGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAg                                                                                     >  1:528855/1‑149 (MQ=255)
                                                                        aaGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                           >  2:411459/1‑151 (MQ=255)
                                                                         aagcttttC‑CAAAAAAGGTTAAGATGTTTTTTTTGTCGCTAAATGCT‑TGAAAAACATGCAGGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAtt                                                                                <  2:614062/142‑1 (MQ=255)
                                                                         aGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgt                                                                             >  1:29645/1‑148 (MQ=255)
                                                                           cGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATGATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGt                                                                         >  2:195984/1‑150 (MQ=255)
                                                                             ttttC‑CAAAACAGGTA‑AGATTATCTTTTTGCCTATTACTGCT‑TTTAAAACAGGTATGGCCAACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                               <  2:514668/142‑1 (MQ=255)
                                                                             tttaC‑CAAAACGGTTT‑AGATGTTTTTTTTGTCGTTTAATGCC‑TAAAAAACGTGCATGCCCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                        <  1:787602/149‑1 (MQ=255)
                                                                              tttC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                     >  2:39047/1‑151 (MQ=255)
                                                                              ttctcaaaaaCAGGTT‑AGATATTTTTTTTGCCG‑TTAATTCCTTGAAAAACAGTCATGACCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTaa                                                                       <  1:150737/142‑1 (MQ=255)
                                                                                tC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                     >  1:96817/1‑149 (MQ=255)
                                                                                     aaaaaGGTTAAGATGTTTTTTTTTTCTTTAAATCCC‑TAAAAAAAATGAAGGAGCAAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCa                                                                    <  1:728059/143‑1 (MQ=255)
                                                                                      aaaaGGTT‑AAATGATTTTTTTGTCGTTTAATCCT‑TAAAAAACATGCATGACCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                               <  1:147737/147‑1 (MQ=255)
                                                                                      aaCAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                              >  2:302682/1‑151 (MQ=255)
                                                                                       aCAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                              >  1:682118/1‑150 (MQ=255)
                                                                                         aGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAt                                                           >  1:246087/1‑151 (MQ=255)
                                                                                           gTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                         >  2:813151/1‑151 (MQ=255)
                                                                                            tt‑AGATTTTTTTTTTGTCGTTAAATGCT‑TAAAAAACATGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                                          <  2:717232/149‑1 (MQ=255)
                                                                                                  tgatTTTTTTGTCGCTTAATGCT‑TGAAAAACATGCATGACCCAAAAAATAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                                    <  2:746155/150‑1 (MQ=255)
                                                                                                     ttttttttTCGTTTAATCTC‑TAAAAAACATCAATGACACAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                                    <  2:33863/145‑1 (MQ=255)
                                                                                                       tttttGTCTATAAATGCT‑AATAAAATATGAAGGACCAAAAAAAGAATAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                <  1:338535/149‑1 (MQ=255)
                                                                                                       tttttGCCCCTAAATTCC‑TAAAAAATATACATGACCAACAAAATAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtca                                               <  2:625693/150‑1 (MQ=255)
                                                                                                              cgtttAATGCC‑TGAAAAACATGAATGACCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                         <  1:296924/146‑1 (MQ=255)
                                                                                                                tttAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTACGTTGACTTTTAGCGCCCATATCTCCAGAATGCTGCCGTTTGTCAGAAATTCGCCGGCAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                       >  1:141081/2‑149 (MQ=255)
                                                                                                                  tAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                   >  2:363484/1‑151 (MQ=255)
                                                                                                                    agtgctgcaaaaaCAGGAATGCGCCAAAAAATAAAAAAAAAACGCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                        <  2:302668/136‑1 (MQ=255)
                                                                                                                    aggCC‑TGAAAAACATGAAGGACCAAAAAAAAAATAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                                  <  1:723198/148‑1 (MQ=255)
                                                                                                                      gCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                                >  1:631085/1‑150 (MQ=255)
                                                                                                                       ccttgaaaaaCAAGAATGAACCAAAAAAAAAAAAAAAAACTTCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                                  <  2:1085561/142‑1 (MQ=255)
                                                                                                                        c‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                             >  2:85648/1‑151 (MQ=255)
                                                                                                                             aaaaCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                          >  2:1136404/1‑150 (MQ=255)
                                                                                                                             aaaaCATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                                <  2:415755/144‑1 (MQ=255)
                                                                                                                             aaaaCAGGCAGGCCCCAAAAAAAAATAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                              <  1:394623/146‑1 (MQ=255)
                                                                                                                              aaaCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                        >  2:1137046/1‑151 (MQ=255)
                                                                                                                              aaaCATGAATGACCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                        <  2:848496/151‑1 (MQ=255)
                                                                                                                                     catgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                  >  2:308029/1‑150 (MQ=255)
                                                                                                                                     catgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt                   >  2:810296/1‑149 (MQ=255)
                                                                                                                                      atgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                  >  2:552986/1‑149 (MQ=255)
                                                                                                                                      atgACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                  <  2:303532/149‑1 (MQ=255)
                                                                                                                                             acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat           >  2:781823/1‑149 (MQ=255)
                                                                                                                                             aaaaaaaaATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat         <  2:156893/149‑1 (MQ=255)
                                                                                                                                             aaaaaaTAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat         <  2:1037826/149‑1 (MQ=255)
                                                                                                                                             aaaaaaTAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAata          <  1:854789/148‑1 (MQ=255)
                                                                                                                                             aaaaaaTAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat           <  1:643421/147‑1 (MQ=255)
                                                                                                                                              aaaaaaaATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt        <  1:329979/150‑1 (MQ=255)
                                                                                                                                               aaaaaaATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAgtgagtga             <  1:416278/145‑1 (MQ=255)
                                                                                                                                               aaaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc       >  1:1058882/1‑151 (MQ=255)
                                                                                                                                               aaaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc       >  2:236700/1‑151 (MQ=255)
                                                                                                                                                  ataataAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct    <  1:1067761/151‑1 (MQ=255)
                                                                                                                                                  aaaataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct    <  1:111976/149‑1 (MQ=255)
                                                                                                                                                    aataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg  >  1:272519/1‑151 (MQ=255)
                                                                                                                                                    aataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg   <  1:80813/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTC‑CAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29347‑29642

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