Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,421,4430AC38.4% 138.4 / 12.4 73Y40S (TAC→TCC) folXD‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base A (44/1);  new base C (0/28);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 2.37e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTG  >  NC_000913/2421296‑2421589
                                                                                                                                                       |                                                                                                                                                  
aTCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTAccc                                                                                                                                                 >  1:1064794/1‑151 (MQ=255)
aTCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTAccc                                                                                                                                                 >  2:1169916/1‑151 (MQ=255)
   cTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCccgc                                                                                                                                              >  1:642297/1‑151 (MQ=255)
           aacaacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATaa                                                                                                                                        >  1:1019465/1‑149 (MQ=255)
              aacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgc                                                                                                                                   >  2:191896/1‑151 (MQ=255)
                cGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgcgc                                                                                                                                 >  2:1198787/1‑151 (MQ=255)
                 ggAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgcgc                                                                                                                                 >  1:1332918/1‑150 (MQ=255)
                    gaACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑TTATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                              >  1:713549/1‑150 (MQ=255)
                         tGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                         >  1:1018144/1‑150 (MQ=255)
                           cAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGc                                                                                                                           >  1:993574/1‑146 (MQ=255)
                            aaTGCCACACCCTGCCTTTTTTATTCTAAAAAAGACCCTTCGTTGGCTTACGTTTAT‑GGAAATAAAGA‑AAAAAAAAATAAA‑CACCCGTAAGG‑ATTTTGTTATCAAGGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAg                                                                                                                       <  2:941322/149‑1 (MQ=255)
                                 ccaccccTGCCTTTTTTTTTCGTATAAAGACCTTTCTTTTGCGTAGTTTTAT‑CGAAATGAAGA‑AAAAAAAATTAAA‑CACCCGAAGGG‑ATTGTTTTATCAAGGTCTGGACCCTCACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                                                                                 <  2:426548/145‑1 (MQ=255)
                                 cacaACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                                                                                 >  1:311794/1‑150 (MQ=255)
                                    aaCCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                              >  2:1225113/1‑150 (MQ=255)
                                      ccTGCCTTTTTTATTCGTAAAAAAACCCTTCTTTTGCGTACTTTTAT‑GGAAATAAAGA‑AAAAAAAAATTAA‑CACCCGTCGGG‑ATATTGTTATAAATGGGACGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                            <  2:1404546/150‑1 (MQ=255)
                                      ccTGCCGTTTTTATTGGTAAAAACACCCTTCTTTTGCATGTGTTTAT‑CAGAATAAAGG‑AAAAACAAATAAA‑CAACCACCAGG‑ATTTTTTTAAAAAGCGGCCATTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                                                                               <  2:1267676/147‑1 (MQ=255)
                                        tGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                         >  1:532094/1‑151 (MQ=255)
                                         gCCGCTATTATCCGTATAAAGTACCTTCGTTTGCGTACGTTTAT‑CGGAATTAATG‑AAGAAGAAATTAA‑CAACTGAGAGG‑TTATTGTTATCAATGTGACTATCCTTTACCCCGCCCATAATGCGCGTTCTTGTGAAGATATCAACTAt                                                                                                            >  1:57786/1‑147 (MQ=255)
                                          ccGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                       >  2:1321138/1‑151 (MQ=255)
                                           cGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                      >  2:114928/1‑151 (MQ=255)
                                           cGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                      >  1:1204046/1‑151 (MQ=255)
                                             cTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGa                                                                                                     >  1:585952/1‑150 (MQ=255)
                                              tattatTCATAAAAAGACTCTTTTTTTGCGTTTTTTTAT‑CGAAATAAAGA‑AAAAAAAAATAAA‑CACCCAGCGGG‑ATTTTTTTATAAATGTGACACCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGa                                                                                                             <  1:223891/141‑1 (MQ=255)
                                                 tattCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                                                                >  1:616147/1‑151 (MQ=255)
                                                    tCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGATATCAATGTGACGATCTACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                                                               >  1:919234/1‑149 (MQ=255)
                                                     cGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                             >  1:1245376/1‑150 (MQ=255)
                                                       tataAAGACCCTTTTTTGGCGT‑CGTTTTTGCGAAATTAAGA‑AAGAAAAAATAAA‑CAGCCGTAAGG‑ATTTTTTTACCATTGTGGGCACCCTCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                                                               <  2:1285051/146‑1 (MQ=255)
                                                           aaaacctttCTTTTGCGTCCTTTTAT‑GGAAATAAAGA‑AAAAAAAAATAAA‑CACCCGTCAGG‑ATTTTGTTACAAATGTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTaa                                                                                      <  2:134932/144‑1 (MQ=255)
                                                             gAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                    >  1:1085886/1‑151 (MQ=255)
                                                               aCCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                                   >  1:305011/1‑150 (MQ=255)
                                                                  ttCGTTTGCTTACTTTTAT‑CGAAATAAAGA‑AAAAAAAAATAAA‑CACCCGTCAGG‑ATATTTTTATCATTGTGACGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                    <  2:567245/146‑1 (MQ=255)
                                                                  ttCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                                >  1:1326424/1‑150 (MQ=255)
                                                                    cGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                              >  2:668011/1‑150 (MQ=255)
                                                                    cGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                              >  1:1202912/1‑150 (MQ=255)
                                                                      tttGCGTACTTTTAT‑GGAAATAAAGA‑AAAAAAAAATTAA‑CACCCGTCAGA‑ATTTTGTTATAAATGTGCCGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                              <  2:1389197/148‑1 (MQ=255)
                                                                      ttggCGTACGTTTAT‑CGAAATAAAGG‑AAAAAAAAATTAA‑CACCCGTAAGG‑ATTTTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                              <  2:675794/145‑1 (MQ=255)
                                                                           gcttttttttatagAAATAAAG‑AAAAAAAAATTAACCACCCGTGAGT‑TTTTTTTTATAAATGTTCCACCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                                                                          <  1:399676/136‑1 (MQ=255)
                                                                            tACTTTTAT‑CGAAATAAAGA‑AAAAAAAAATAAA‑CACCCGTCAGG‑TTATTTTAAACAATGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                         <  2:533315/147‑1 (MQ=255)
                                                                            tACTTTTAG‑CGAAAAAAAGG‑AAAAAAAAATAAA‑CACCCGTCAGG‑ATTTTTTAATAAATGTCACTCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                                                                          <  1:1161198/146‑1 (MQ=255)
                                                                            tACGTTTAT‑GAAAATAAAGA‑AAAAAAAAATAAA‑CACCCGTCAGG‑ATATTGTTAAAATTTTGACCCCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCa                                                                      <  1:615707/150‑1 (MQ=255)
                                                                                 tttatgaaaaTAAAGA‑AAAAAAAAATAAA‑CACCCGGCGGT‑ATTTTGTTATCAATGGGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                  <  2:455155/142‑1 (MQ=255)
                                                                                 ttAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                >  1:1345804/1‑151 (MQ=255)
                                                                                   aT‑CGAAATAAAGA‑AAAAAGAAAATAC‑CAACCGCCAGG‑ATATTTTTACCAATGTGACGCTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                                                               <  2:625397/150‑1 (MQ=255)
                                                                                      cGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                              >  1:1204884/1‑149 (MQ=255)
                                                                                          aatAAGGAAAGAAAAAATTAA‑CACCCGTCGGG‑ATTTTGTTACAAATGTGACCACCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                           <  2:1399035/147‑1 (MQ=255)
                                                                                             aaGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                     >  1:658624/1‑151 (MQ=255)
                                                                                                g‑aagaagaaATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGtt                                                  >  2:1265869/1‑151 (MQ=255)
                                                                                                g‑aagaagaaATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGtt                                                  >  1:1027337/1‑151 (MQ=255)
                                                                                                      agaaATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                             >  2:855018/1‑151 (MQ=255)
                                                                                                       gaaATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                            >  1:880081/1‑151 (MQ=255)
                                                                                                       aaaaaTAA‑CAACCCTTAGGAATTTTTTTATCAAGGTGACTACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                             <  1:29075/150‑1 (MQ=255)
                                                                                                         aaTTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                          >  2:554138/1‑151 (MQ=255)
                                                                                                         aaTTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                          >  1:600501/1‑151 (MQ=255)
                                                                                                             aa‑caccTGTAAGG‑TTATTGTTTAAAATGTGACAACCCACACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                       <  2:395540/150‑1 (MQ=255)
                                                                                                               accaccctatagtttatTTTTATCAATGTGACGCCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                       <  1:1123156/136‑1 (MQ=255)
                                                                                                                  acCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGCTGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCGGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                                   >  2:1301316/1‑150 (MQ=255)
                                                                                                                   ccGTAAGG‑TTTTTTTTATAAATGAGACAATCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaa                                  <  1:1411358/150‑1 (MQ=255)
                                                                                                                          g‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                          >  2:1098348/1‑151 (MQ=255)
                                                                                                                              atTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAgg                       >  2:183962/1‑151 (MQ=255)
                                                                                                                               ttgtaaaaattgtgACGACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTa                              <  2:304726/133‑1 (MQ=255)
                                                                                                                               ttGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAgg                       >  2:120768/1‑150 (MQ=255)
                                                                                                                               ttGTTAACAATGTCACCCTCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAgg                       <  2:40766/150‑1 (MQ=255)
                                                                                                                                  ttATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAt                     <  2:1019529/149‑1 (MQ=255)
                                                                                                                                   tATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg                  >  1:1307225/1‑151 (MQ=255)
                                                                                                                                    agaaaTGGGATGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtg                    <  1:432355/145‑1 (MQ=255)
                                                                                                                                      aaaTGCGACCACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCt                <  1:436590/149‑1 (MQ=255)
                                                                                                                                          gtgACGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGat            <  2:1244136/150‑1 (MQ=255)
                                                                                                                                           tgccGATCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGa             <  1:369308/145‑1 (MQ=255)
                                                                                                                                           tgACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGatat          >  2:368974/1‑151 (MQ=255)
                                                                                                                                            gACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATc         >  1:1170646/1‑151 (MQ=255)
                                                                                                                                               caTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCa      <  1:1287238/150‑1 (MQ=255)
                                                                                                                                                  ccACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGt   >  2:269571/1‑151 (MQ=255)
                                                                                                                                                   cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg  >  2:1287390/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
ATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTAT‑CGGAATTAAGG‑AAGAAGAAATTAA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTG  >  NC_000913/2421296‑2421589

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