Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,421,4430AC43.5% 152.7 / 29.8 93Y40S (TAC→TCC) folXD‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base A (43/9);  new base C (0/40);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 1.69e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGT  >  NC_000913/2421295‑2421588
                                                                                                                                                      |                                                                                                                                                 
cATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAcc                                                                                                                                                 <  1:1467885/151‑1 (MQ=255)
   ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccg                                                                                                                                              >  1:1017648/1‑151 (MQ=255)
      aaGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgccg                                                                                                                                           >  1:806961/1‑151 (MQ=255)
      aaGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgccg                                                                                                                                           >  1:179667/1‑151 (MQ=255)
         tAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGAt                                                                                                                                         >  2:1164785/1‑150 (MQ=255)
                    agaCCCTTCAAGGCCAAACCCCGCCGTTATTATTGTTAAAAAAACCCTTCTTTTGTATTCTTTTATCGTAATTAAGAAAAAAAAAATAAACACCCGTCAGGTTATTTTTATCATTGTGCCTACCCACACCCCCGCCGATAAAGCGCGCACt                                                                                                                             <  1:767485/151‑1 (MQ=255)
                    aaaaCCTGCAAGGGCACACCCTGCCGTTATTATTCGAATAAAAAACCTTCTTTGGCGTACTTTTATCGAAATTAGGAAAAAAAAAATAAACACCCGCCAGGATATTTTTACCAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACt                                                                                                                             <  1:237779/149‑1 (MQ=255)
                      aaCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                             >  2:233019/1‑149 (MQ=255)
                       aCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGc                                                                                                                          >  2:1095255/1‑151 (MQ=255)
                        cctaaaaTGGCAAACCCTGCCGTTATTATTCGTATAAAGACCTTTCGTTTGCGTACTTTTATGGAAATAAAGGAAAAAAAAATAAACACCCGTAAGGATATTTTTACAAATGTGACCACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                         <  2:208467/146‑1 (MQ=255)
                           gaaaggccacaCCCTGCCTTTTTTATTCGTATAAAAACCCTTCTTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTCAGGATATTTTTATCAGTGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAg                                                                                                                      <  1:841283/144‑1 (MQ=255)
                           gaaaTGCCACACCCTGCCTTTATTTTTCGTAAAAAGACCCTTCGTTTGCGTACGTTTATCGAAATAAAGAAAGAAAAAATAACCACCCGTCAGAATATTTTTATCAAGGGGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAg                                                                                                                      <  2:747010/149‑1 (MQ=255)
                           gCAAGGGCACACCCTGCCATTATTATTCTTAAAAAAAACCTTTGTTGGCGTACGTTTTTCGGAATTAGGGAAAAAGAAATTAACAACCGTAAGGATATTTTTACTAATTGGGGATTCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                        <  1:1399229/149‑1 (MQ=255)
                            aaaTGGCACACCCTGCAAATTTTATTTGAAAAAAAAACCTTGTTTTGCTTGCGTTTATAGAAATAAAGAAAAAAGAAAATAACAACCGCCAGTTTGTTTTTATAATGGTGACGATCCGCAACCCCGCCGATAAAGCGCGCACTAGCGAAGa                                                                                                                     <  1:686188/150‑1 (MQ=255)
                             aagggCACACCCTGCCGTTTTTATTCGTAAAAAGAACTTTCTTTTGCGTCCTTTTACCGAAATAAAGAAAAAAAAAATAAACACCCGTCAGTATTTTTTTACAAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGat                                                                                                                    <  1:1198556/148‑1 (MQ=255)
                             aaTGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                        >  1:690444/1‑147 (MQ=255)
                                gccacaCCCTCCCGTTATTATTCGTAAAAAAACCCTTCTTTTGCGTCCGTTTATCGAAATAAAGAAAAAAAAAATTAACACCCGTCAGGATATTGTTATAAATGGGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATc                                                                                                                 <  2:987070/149‑1 (MQ=255)
                                      aCCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                           >  2:202606/1‑151 (MQ=255)
                                          gCCGTTATTATTCGTAAAAAAAACTTTCTTTTGCTTACTTTTATCGAAATAAAGAAAGAAAAAATAAACACCCGTAAGGATATTGTTATCAATGTGCCGACCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                                                         <  1:829634/149‑1 (MQ=255)
                                              ttattatTCTAAAAAAGACCTTTCTTTTGCGTACTTTTAGCGAAATTAAGAAAGAAAAAATAAACACCCGTAAGGATTTTGTTAAAAATGTGACTACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGa                                                                                                    <  2:583961/149‑1 (MQ=255)
                                               tattatTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                                                                                  >  1:720511/1‑151 (MQ=255)
                                                 ttattCAAAAAAAGACCTTTCTTTTGGTTACTTTTACCGAAATTAAGAAAAAAAAAATAAAAACCCGCAAGAATATTGTTATAAAGGAGACGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                                 <  2:297407/150‑1 (MQ=255)
                                                  tattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                                                               >  2:145261/1‑151 (MQ=255)
                                                     tCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                            >  1:229976/1‑151 (MQ=255)
                                                        aaaaaaaacccTTCTTTTGCGTACTTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGTCAGGATTTTTTAATAAAGGTGAGGTCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAcc                                                                                         <  1:746260/142‑1 (MQ=255)
                                                         ataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                        >  1:617746/1‑151 (MQ=255)
                                                          taAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                       >  2:246318/1‑151 (MQ=255)
                                                           aaaGACCCTTCTTTGGCGTACGTTTATGGAAATAAAGAAAAAAAAAATTAACACCCGTAAGGATATTGTTATCAATGTGACGATCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                      <  1:106815/151‑1 (MQ=255)
                                                           aaaGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                       >  2:550154/1‑150 (MQ=255)
                                                             aaaccttttttttGCGTACGTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGCAAGAATATTGTTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                    <  1:778209/141‑1 (MQ=255)
                                                               acccTTCTTTTGCGTACTTTTATCGAATTTAAGAAAAAAAAAATAAACACCCGCCAGGATTTTTTTATCAATGTGCCGTTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                                  <  1:892345/149‑1 (MQ=255)
                                                                 ccttgttttGCGTAATTTTCTCAAAATAAAGAAAGAAGAAATTAACAACCGTCAGTTTATTATTATAAATGTGCCGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                                <  1:736226/145‑1 (MQ=255)
                                                                        ttGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTCAGGATTTTTTTATAAAGGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                          <  1:877512/150‑1 (MQ=255)
                                                                        ttGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                          >  2:461358/1‑150 (MQ=255)
                                                                         tgggtccttttaaaagAATTAAGAAAGAAGAAATCAACACCCAGAGGGTTATTGTTATAAATGGGACCCGCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                                                                         <  2:1458177/135‑1 (MQ=255)
                                                                         tGCTTACTTTTAGCGAAATAAAGAAAAAAAAAAAAAACACCCGTGAGAATATTGTTACAATTGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                                                                         <  1:580849/150‑1 (MQ=255)
                                                                          gCGTTCTTTTATGGAAATAAAGGAAAAAAAAATAAAAACCCGTCAGGATATTTTTATCAATGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                        <  2:646864/150‑1 (MQ=255)
                                                                          gCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattatt                                                                       >  2:685845/1‑151 (MQ=255)
                                                                          gCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattatt                                                                       >  2:540330/1‑151 (MQ=255)
                                                                          gCGTACGTTTATCGAAATAAAGGAAAAAAAAATAAACAACCGTCAGAATTTTTTTATCAATGTGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                        <  1:1091217/150‑1 (MQ=255)
                                                                           cGTACTTTTATCGAAATAAAGGAAAAAAAAATTAACACCCGTCAGGATATTGTTATCAATGGGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                      <  1:681015/151‑1 (MQ=255)
                                                                            gTATGTTTATCGGAATTAAGAAAAAAAAAATTACCAACCGTCAGAATTTTTTTATCAAGGGCACGTCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCa                                                                     <  1:570597/151‑1 (MQ=255)
                                                                            gTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                      >  2:634621/1‑150 (MQ=255)
                                                                             tactttgatggaaaTTAAGGAAGAAAAAATAACCACCCGAAAGGATTTTTTTCAAAGGGTGACGTTCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                      <  1:1147496/137‑1 (MQ=255)
                                                                              ccttttATCGAAATAAAGAAAAAAAAAATAAACACCCTCCAGTATTTTTTAATCAATGTGACGACCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg                                                                    <  1:572665/147‑1 (MQ=255)
                                                                              acttttATGGAAATAAAGAAAGAAAAAATAAACACCCGTCAGAATATTTTTATAAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                   <  1:1303819/148‑1 (MQ=255)
                                                                               tttttATCGAAATAAAGAAAAAAAAAATAAACAACCGTAAGGTTTTTGTTATAAAGGTGACGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                   <  2:330727/148‑1 (MQ=255)
                                                                                gTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                 >  1:1025037/1‑151 (MQ=255)
                                                                                gTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                 >  2:118820/1‑151 (MQ=255)
                                                                                 tttAGCGAAATAAGGAAAAAAAAAATAAACAACCGTAAGGATATTGTTATCAATGGGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                 <  1:730656/150‑1 (MQ=255)
                                                                                  ttATGGAAATAAAGGAAAAAAAAATAAACACCCGCAAGGATTTTTTTATAATTGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg                                                                <  1:1086470/150‑1 (MQ=255)
                                                                                  ttATGGAAATAAAGAAAAAAAAAATAACCAACCGTAAGGATTTTGTTATAAATGTGACCACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg                                                                <  2:1279238/150‑1 (MQ=255)
                                                                                  ttATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                               >  1:1085137/1‑151 (MQ=255)
                                                                                   attCGTAATGAAGGAAAAAAAAATGAACACCCAGAAGTATATTTTTATAAATTTGATGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                                                              <  2:1375290/149‑1 (MQ=255)
                                                                                    atggaaaTAAAGAAAAAAAAAATAAACACCCGTCAGGATATTTTTACAAATGTGACACCCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                                                              <  1:674634/145‑1 (MQ=255)
                                                                                    aTCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                               >  2:1332920/1‑149 (MQ=255)
                                                                                    aTCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                               >  2:798481/1‑149 (MQ=255)
                                                                                     tCGAAATAAAGAAAAAAAAAATTACCACCCGTAAGTATTTTTTTATAAATGTGACCACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                             <  2:810148/150‑1 (MQ=255)
                                                                                     agcaaaatAAAGAAAAAAAAAATTACCACCGTAAAGGATATTTTTATAAATGTGACGTTCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                                                              <  2:1104558/143‑1 (MQ=255)
                                                                                      ggaaaTAAAGAAAAAAAAAATAAACAACCGTAAGAATTTTTTTATAAATGGGCCATCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTaga                                                            <  2:1120905/147‑1 (MQ=255)
                                                                                      cGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                             >  1:479862/1‑149 (MQ=255)
                                                                                       ggAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                             >  1:1449603/1‑148 (MQ=255)
                                                                                         aaTAAAGAAAAAAAAAAAAAACACCCGTCAGTATTTTTTTAAAAATGTGACTCCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAAt                                                        <  2:171619/151‑1 (MQ=255)
                                                                                          ataaaGAAAAAAAAAATTAACACCCGTAAGGTTTTTGTTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa                                                       <  2:169714/148‑1 (MQ=255)
                                                                                           ttAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa                                                       >  2:761609/1‑150 (MQ=255)
                                                                                            tAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                     >  2:1202713/1‑151 (MQ=255)
                                                                                            aaaGAAAAAAAAAATAAACACCCGTCAGGATATTTTTACAAATTTGACACCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATaa                                                      <  2:647230/149‑1 (MQ=255)
                                                                                             aaGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                    >  1:1127021/1‑151 (MQ=255)
                                                                                                aaagaaaaaATAAACACCCGGCAGAATATTTTTATAAATGGGACCACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGtt                                                 <  1:1276589/150‑1 (MQ=255)
                                                                                                  aaaaaaaaTTAACAACCTCGAGGTTGTTTTTACAAATGTGAGTACCCCCAGCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtc                                               <  2:1265544/146‑1 (MQ=255)
                                                                                                   gaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtct                                              >  2:208632/1‑151 (MQ=255)
                                                                                                       aaaTTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                          >  1:858243/1‑151 (MQ=255)
                                                                                                       aaaTTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                          >  2:314687/1‑151 (MQ=255)
                                                                                                         aataaccaccggtaaGGATTTTTTTATCAATGTGACGCTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                         <  2:424806/137‑1 (MQ=255)
                                                                                                            aacaacCGTGAGTTTATTTTAATAATTGTGCCAATCCACCACGCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                     <  1:1316521/151‑1 (MQ=255)
                                                                                                                 ccGTCAGGATATTTTTATAAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaaa                                <  1:1335147/151‑1 (MQ=255)
                                                                                                                       aGGATATTTTTAACAATGTGCCGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAAc                           <  1:734284/151‑1 (MQ=255)
                                                                                                                       aGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAAc                           >  1:836806/1‑151 (MQ=255)
                                                                                                                         aatttttttAACAATTTGACGACCCACGCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                         <  1:897561/144‑1 (MQ=255)
                                                                                                                              tGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAt                    >  1:597178/1‑151 (MQ=255)
                                                                                                                               tttATAAATGGGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAt                    <  1:300617/149‑1 (MQ=255)
                                                                                                                               gTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtg                   >  2:1186108/1‑151 (MQ=255)
                                                                                                                                 tATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg                 >  2:320989/1‑151 (MQ=255)
                                                                                                                                  aTCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGc                >  1:988072/1‑151 (MQ=255)
                                                                                                                                  aTCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGc                >  1:993568/1‑151 (MQ=255)
                                                                                                                                      aTGTGACGATCCACTACCCCGCCGATAAAGCGCGGACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg             >  2:366764/1‑150 (MQ=255)
                                                                                                                                      aTGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg             >  1:1071165/1‑150 (MQ=255)
                                                                                                                                          gaagATCCACTCCCCCGCCGATACAGCGCTCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATACCCGTCTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGCTATc        <  2:789117/148‑1 (MQ=255)
                                                                                                                                           ccGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCg       <  2:229980/150‑1 (MQ=255)
                                                                                                                                           aCGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCg       >  2:700095/1‑151 (MQ=255)
                                                                                                                                            cGATCCACTACCCCGCCGATAAAGCGCGCACTGGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc      >  2:1284348/1‑151 (MQ=255)
                                                                                                                                                 ccACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGt  >  2:219801/1‑151 (MQ=255)
                                                                                                                                                 ccACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGt  >  2:865366/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGT  >  NC_000913/2421295‑2421588

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