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

TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29347‑29642
                                                                                                                                                        |                                                                                                                                                  
ttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAatat                                                                                                                                                    <  1:466255/150‑1 (MQ=255)
 tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                                   <  1:431582/150‑1 (MQ=255)
  aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                                 <  1:462129/151‑1 (MQ=255)
     gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaaa                                                                                                                                               >  2:800600/1‑150 (MQ=255)
       aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATc                                                                                                                                            <  1:592883/151‑1 (MQ=255)
        aaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAt                                                                                                                                             <  1:271255/149‑1 (MQ=255)
         aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                                          >  2:572552/1‑151 (MQ=255)
         aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                                          <  2:471574/151‑1 (MQ=255)
             tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                       >  2:80816/1‑150 (MQ=255)
              gACCATTTGGTCCACTTTTTTCTGCTCTTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGTTTATATGATCTTTTTGTCCCTGAATGCCGGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCa                                                                                                                                      <  1:1092451/150‑1 (MQ=255)
                ccATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCAtt                                                                                                                                    >  2:977391/1‑150 (MQ=255)
                 cATTTGGGCCACTTTTTTCTGCCCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAAAACAGGTTAGATGATCTTTTTTTCGCTTAATGCCTGCAAAACTTGCATGAGCCACAAAATAATATAAAAAATCCCGCCa                                                                                                                                      <  1:600388/147‑1 (MQ=255)
                  aTTTGGTCCACTTTTTTATGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                   >  2:579875/1‑149 (MQ=255)
                     tGGTCCTTTTTTTTGTTTTCTTTTTTATTTCATGCAATCTTCTTGCTGCGCAGGGTTTTTCAAAAAAAGGTTATATGATTTTTTTTTCGTTTAATCCCTAAAAAACAGGTATGACCCAAAAAAAAATAAAAAAATCCCCGCCATTAAGtt                                                                                                                               <  1:826078/150‑1 (MQ=255)
                       gTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGAc                                                                                                                            >  2:505746/1‑151 (MQ=255)
                         ccctttttttCTGCTCTTTTTTTTTTCATGCATTTTTCTTGCTGCGCAAGCTTTTTCCAAAACAGGTTAGATGATTTTTTTGTCGTTAAATCCCTGAAAAACATGCAGGACCCAAAAAATAATATAAAAAATCCCGCCATTAAGtt                                                                                                                               <  1:148416/142‑1 (MQ=255)
                            tttttttCTCCTCTTTTTTTTTTCAGGAATTTTTCTGGCGGCCAAAGGTTTTCCCAAAAAAGGTTAGATGATTTTTTTGTCGCTAAAGGCTTAAAAAACAGGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                         <  2:383856/148‑1 (MQ=255)
                             tttttttcccTCGTTTTTATTTCATGCATTCTTCTTGCGGCGAAGGTTTTTCAAAAAACGGGTAAGAGTATTTTTTTGACGTTTAATGCGTCTAAAACATCCATGACCCAAAAAATAAAAAAAAAATTCCCGCCATTAAGTTGACTTTTa                                                                                                                       <  1:201734/141‑1 (MQ=255)
                             ttttttCTGCTCTTTTTTATTTCATGCATTCTTCTTGCGGCGCAAGCGTTTCCAAAAACGGTTGAGATAATTTTTTTGTCGTTAAATGCCTAAAAAACATGCAGGACCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTa                                                                                                                       <  2:822192/150‑1 (MQ=255)
                             ttttttCTGCTCTTTTTTATTTCATGCAATCTTCTTGCTGCGCAGGCTTTTTCCAAAACAGGTTAGATGATTTTTTTGTCGCTTAATGCCTGAAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                        <  2:645916/149‑1 (MQ=255)
                              tttttCTTCTTTTTTTTTTTTCATGCAATCTTCTTGCTGGCCAGCCGTTTACAAAAACAGGTTAAATAATCTTTTTGTCACTAATTGCCAGTAAAACATGCATCGGCCAAAAAATAATATAAAAACTCCCGCCATTAAGTTGACTTTTAg                                                                                                                      <  1:450310/150‑1 (MQ=255)
                                   cTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                                 >  1:1013627/1‑150 (MQ=255)
                                   cTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGAGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                                 >  1:279985/1‑150 (MQ=255)
                                      cttttttttATTTCATGCAATTTTTTTGCTGCGCAGGCGTTTCCCAAAAGAGTTTAGATGTTTTTTTTGTCGCTTAATCCCTATAAACCATGCAGGAGCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCa                                                                                                                <  2:811196/144‑1 (MQ=255)
                                       tttttttttttTCATGCATTCTTTTTGCTGCGCAAGCTTTTCCAAGAACAGTTTAGATGTTTTTTTTGCCGTTTAATGCCTGAAAAACATGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGcc                                                                                                                  <  2:621081/142‑1 (MQ=255)
                                        cGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                           >  2:1021362/1‑151 (MQ=255)
                                            tttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                         >  2:645828/1‑149 (MQ=255)
                                             ttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                         >  2:257185/1‑148 (MQ=255)
                                                 ttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgc                                                                                                  >  2:329795/1‑151 (MQ=255)
                                                        aaTCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                           >  2:93910/1‑151 (MQ=255)
                                                          tttttttGCTGCCAAGGCTTTTTCAAGAACGGTTTAGATGTTTTTTTTTTCGTTTAATCCTTGAAAAAAATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                              <  2:175417/141‑1 (MQ=255)
                                                          tcttctTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTg                                                                                         >  1:1095997/1‑151 (MQ=255)
                                                            ttcttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                       >  2:877301/1‑151 (MQ=255)
                                                              cttGCTGCGCAAGCGTTTGCAAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                       >  1:796909/1‑149 (MQ=255)
                                                                tgctgcGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAg                                                                                     >  1:528855/1‑149 (MQ=255)
                                                                        aaGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                           >  2:411459/1‑151 (MQ=255)
                                                                         aagcttttCCAAAAAAGGTTAAGATGTTTTTTTTGTCGCTAAATGCTTGAAAAACATGCAGGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAtt                                                                                <  2:614062/142‑1 (MQ=255)
                                                                         aGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgt                                                                             >  1:29645/1‑148 (MQ=255)
                                                                           cGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATGATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGt                                                                         >  2:195984/1‑150 (MQ=255)
                                                                             ttttCCAAAACAGGTAAGATTATCTTTTTGCCTATTACTGCTTTTAAAACAGGTATGGCCAACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                               <  2:514668/142‑1 (MQ=255)
                                                                             tttaCCAAAACGGTTTAGATGTTTTTTTTGTCGTTTAATGCCTAAAAAACGTGCATGCCCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                        <  1:787602/149‑1 (MQ=255)
                                                                              tttCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                     >  2:39047/1‑151 (MQ=255)
                                                                              ttctcaaaaaCAGGTTAGATATTTTTTTTGCCGTTAATTCCTTGAAAAACAGTCATGACCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTaa                                                                       <  1:150737/142‑1 (MQ=255)
                                                                                tCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                     >  1:96817/1‑149 (MQ=255)
                                                                                     aaaaaGGTTAAGATGTTTTTTTTTTCTTTAAATCCCTAAAAAAAATGAAGGAGCAAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCa                                                                    <  1:728059/143‑1 (MQ=255)
                                                                                      aaaaGGTTAAATGATTTTTTTGTCGTTTAATCCTTAAAAAACATGCATGACCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                               <  1:147737/147‑1 (MQ=255)
                                                                                      aaCAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                              >  2:302682/1‑151 (MQ=255)
                                                                                       aCAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                              >  1:682118/1‑150 (MQ=255)
                                                                                         aGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAt                                                           >  1:246087/1‑151 (MQ=255)
                                                                                           gTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                         >  2:813151/1‑151 (MQ=255)
                                                                                            ttAGATTTTTTTTTTGTCGTTAAATGCTTAAAAAACATGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                                          <  2:717232/149‑1 (MQ=255)
                                                                                                  tgatTTTTTTGTCGCTTAATGCTTGAAAAACATGCATGACCCAAAAAATAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                                    <  2:746155/150‑1 (MQ=255)
                                                                                                     ttttttttTCGTTTAATCTCTAAAAAACATCAATGACACAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                                    <  2:33863/145‑1 (MQ=255)
                                                                                                       tttttGTCTATAAATGCTAATAAAATATGAAGGACCAAAAAAAGAATAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                <  1:338535/149‑1 (MQ=255)
                                                                                                       tttttGCCCCTAAATTCCTAAAAAATATACATGACCAACAAAATAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtca                                               <  2:625693/150‑1 (MQ=255)
                                                                                                              cgtttAATGCCTGAAAAACATGAATGACCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                         <  1:296924/146‑1 (MQ=255)
                                                                                                                tttAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTACGTTGACTTTTAGCGCCCATATCTCCAGAATGCTGCCGTTTGTCAGAAATTCGCCGGCAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                       >  1:141081/2‑149 (MQ=255)
                                                                                                                  tAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                   >  2:363484/1‑151 (MQ=255)
                                                                                                                    agtgctgcaaaaaCAGGAATGCGCCAAAAAATAAAAAAAAAACGCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                        <  2:302668/136‑1 (MQ=255)
                                                                                                                    aggCCTGAAAAACATGAAGGACCAAAAAAAAAATAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                                  <  1:723198/148‑1 (MQ=255)
                                                                                                                      gCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                                >  1:631085/1‑150 (MQ=255)
                                                                                                                       ccttgaaaaaCAAGAATGAACCAAAAAAAAAAAAAAAAACTTCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                                  <  2:1085561/142‑1 (MQ=255)
                                                                                                                        cTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                             >  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)
                                                                                                                                                        |                                                                                                                                                  
TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  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: