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

ATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTG  >  NC_000913/2421296‑2421589
                                                                                                                                                       |                                                                                                                                                  
aTCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAccc                                                                                                                                                 >  1:1064794/1‑151 (MQ=255)
aTCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAccc                                                                                                                                                 >  2:1169916/1‑151 (MQ=255)
   cTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgc                                                                                                                                              >  1:642297/1‑151 (MQ=255)
           aacaacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATaa                                                                                                                                        >  1:1019465/1‑149 (MQ=255)
              aacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgc                                                                                                                                   >  2:191896/1‑151 (MQ=255)
                cGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgcgc                                                                                                                                 >  2:1198787/1‑151 (MQ=255)
                 ggAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgcgc                                                                                                                                 >  1:1332918/1‑150 (MQ=255)
                    gaACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGTTATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                              >  1:713549/1‑150 (MQ=255)
                         tGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                         >  1:1018144/1‑150 (MQ=255)
                           cAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGc                                                                                                                           >  1:993574/1‑146 (MQ=255)
                            aaTGCCACACCCTGCCTTTTTTATTCTAAAAAAGACCCTTCGTTGGCTTACGTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGTAAGGATTTTGTTATCAAGGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAg                                                                                                                       <  2:941322/149‑1 (MQ=255)
                                 ccaccccTGCCTTTTTTTTTCGTATAAAGACCTTTCTTTTGCGTAGTTTTATCGAAATGAAGAAAAAAAAATTAAACACCCGAAGGGATTGTTTTATCAAGGTCTGGACCCTCACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                                                                                 <  2:426548/145‑1 (MQ=255)
                                 cacaACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                                                                                 >  1:311794/1‑150 (MQ=255)
                                    aaCCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                              >  2:1225113/1‑150 (MQ=255)
                                      ccTGCCTTTTTTATTCGTAAAAAAACCCTTCTTTTGCGTACTTTTATGGAAATAAAGAAAAAAAAAATTAACACCCGTCGGGATATTGTTATAAATGGGACGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                            <  2:1404546/150‑1 (MQ=255)
                                      ccTGCCGTTTTTATTGGTAAAAACACCCTTCTTTTGCATGTGTTTATCAGAATAAAGGAAAAACAAATAAACAACCACCAGGATTTTTTTAAAAAGCGGCCATTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                                                                               <  2:1267676/147‑1 (MQ=255)
                                        tGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                         >  1:532094/1‑151 (MQ=255)
                                         gCCGCTATTATCCGTATAAAGTACCTTCGTTTGCGTACGTTTATCGGAATTAATGAAGAAGAAATTAACAACTGAGAGGTTATTGTTATCAATGTGACTATCCTTTACCCCGCCCATAATGCGCGTTCTTGTGAAGATATCAACTAt                                                                                                            >  1:57786/1‑147 (MQ=255)
                                          ccGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                       >  2:1321138/1‑151 (MQ=255)
                                           cGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                      >  2:114928/1‑151 (MQ=255)
                                           cGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                      >  1:1204046/1‑151 (MQ=255)
                                             cTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGa                                                                                                     >  1:585952/1‑150 (MQ=255)
                                              tattatTCATAAAAAGACTCTTTTTTTGCGTTTTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCAGCGGGATTTTTTTATAAATGTGACACCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGa                                                                                                             <  1:223891/141‑1 (MQ=255)
                                                 tattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                                                                >  1:616147/1‑151 (MQ=255)
                                                    tCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGATATCAATGTGACGATCTACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                                                               >  1:919234/1‑149 (MQ=255)
                                                     cGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                             >  1:1245376/1‑150 (MQ=255)
                                                       tataAAGACCCTTTTTTGGCGTCGTTTTTGCGAAATTAAGAAAGAAAAAATAAACAGCCGTAAGGATTTTTTTACCATTGTGGGCACCCTCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                                                               <  2:1285051/146‑1 (MQ=255)
                                                           aaaacctttCTTTTGCGTCCTTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGTCAGGATTTTGTTACAAATGTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTaa                                                                                      <  2:134932/144‑1 (MQ=255)
                                                             gAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                    >  1:1085886/1‑151 (MQ=255)
                                                               aCCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                                   >  1:305011/1‑150 (MQ=255)
                                                                  ttCGTTTGCTTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTCAGGATATTTTTATCATTGTGACGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                    <  2:567245/146‑1 (MQ=255)
                                                                  ttCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                                >  1:1326424/1‑150 (MQ=255)
                                                                    cGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                              >  2:668011/1‑150 (MQ=255)
                                                                    cGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                              >  1:1202912/1‑150 (MQ=255)
                                                                      tttGCGTACTTTTATGGAAATAAAGAAAAAAAAAATTAACACCCGTCAGAATTTTGTTATAAATGTGCCGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                              <  2:1389197/148‑1 (MQ=255)
                                                                      ttggCGTACGTTTATCGAAATAAAGGAAAAAAAAATTAACACCCGTAAGGATTTTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                              <  2:675794/145‑1 (MQ=255)
                                                                           gcttttttttatagAAATAAAGAAAAAAAAATTAACCACCCGTGAGTTTTTTTTTATAAATGTTCCACCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                                                                          <  1:399676/136‑1 (MQ=255)
                                                                            tACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTCAGGTTATTTTAAACAATGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                         <  2:533315/147‑1 (MQ=255)
                                                                            tACTTTTAGCGAAAAAAAGGAAAAAAAAATAAACACCCGTCAGGATTTTTTAATAAATGTCACTCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                                                                          <  1:1161198/146‑1 (MQ=255)
                                                                            tACGTTTATGAAAATAAAGAAAAAAAAAATAAACACCCGTCAGGATATTGTTAAAATTTTGACCCCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCa                                                                      <  1:615707/150‑1 (MQ=255)
                                                                                 tttatgaaaaTAAAGAAAAAAAAAATAAACACCCGGCGGTATTTTGTTATCAATGGGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                  <  2:455155/142‑1 (MQ=255)
                                                                                 ttATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                >  1:1345804/1‑151 (MQ=255)
                                                                                   aTCGAAATAAAGAAAAAAGAAAATACCAACCGCCAGGATATTTTTACCAATGTGACGCTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                                                               <  2:625397/150‑1 (MQ=255)
                                                                                      cGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                              >  1:1204884/1‑149 (MQ=255)
                                                                                          aatAAGGAAAGAAAAAATTAACACCCGTCGGGATTTTGTTACAAATGTGACCACCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                           <  2:1399035/147‑1 (MQ=255)
                                                                                             aaGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                     >  1:658624/1‑151 (MQ=255)
                                                                                                gaagaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGtt                                                  >  2:1265869/1‑151 (MQ=255)
                                                                                                gaagaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGtt                                                  >  1:1027337/1‑151 (MQ=255)
                                                                                                      agaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                             >  2:855018/1‑151 (MQ=255)
                                                                                                       gaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                            >  1:880081/1‑151 (MQ=255)
                                                                                                       aaaaaTAACAACCCTTAGGAATTTTTTTATCAAGGTGACTACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                             <  1:29075/150‑1 (MQ=255)
                                                                                                         aaTTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                          >  2:554138/1‑151 (MQ=255)
                                                                                                         aaTTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                          >  1:600501/1‑151 (MQ=255)
                                                                                                             aacaccTGTAAGGTTATTGTTTAAAATGTGACAACCCACACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                       <  2:395540/150‑1 (MQ=255)
                                                                                                               accaccctatagtttatTTTTATCAATGTGACGCCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                       <  1:1123156/136‑1 (MQ=255)
                                                                                                                  acCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGCTGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCGGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                                   >  2:1301316/1‑150 (MQ=255)
                                                                                                                   ccGTAAGGTTTTTTTTATAAATGAGACAATCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaa                                  <  1:1411358/150‑1 (MQ=255)
                                                                                                                          gATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                          >  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)
                                                                                                                                                       |                                                                                                                                                  
ATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTG  >  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: