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

CATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGT  >  NC_000913/2421295‑2421588
                                                                                                                                                      |                                                                                                                                                 
cATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTAcc                                                                                                                                                 <  1:1467885/151‑1 (MQ=255)
   ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCccg                                                                                                                                              >  1:1017648/1‑151 (MQ=255)
      aaGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCccgccg                                                                                                                                           >  1:806961/1‑151 (MQ=255)
      aaGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCccgccg                                                                                                                                           >  1:179667/1‑151 (MQ=255)
         tAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGAt                                                                                                                                         >  2:1164785/1‑150 (MQ=255)
                    agaCCCTTCAAGGCCAAACCCCGCCGTTATTATTGTTAAAAAAACCCTTCTTTTGTATTCTTTTATCGTAATTAAGAAAAAAAAAATAAACACCCGTC‑AGGTTATTTTTATCATTGTGCCTA‑CCCACACCCCCGCCGATAAAGCGCGCACt                                                                                                                             <  1:767485/151‑1 (MQ=255)
                    aaaaCCTGCAAGGGCACACCCTGCCGTTATTATTCGAATAAAAAACCTTCTTTGGCGTACTTTTATCGAAATTAGGAAAAAAAAAATAAACACCCGCC‑AGGATATTTTTACCAATGTGACGA‑TCCACTCCCCCGCCGATAAAGCGCGCACt                                                                                                                             <  1:237779/149‑1 (MQ=255)
                      aaCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                             >  2:233019/1‑149 (MQ=255)
                       aCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGc                                                                                                                          >  2:1095255/1‑151 (MQ=255)
                        cctaaaaTGGCAAACCCTGCCGTTATTATTCGTATAAAGACCTTTCGTTTGCGTACTTTTATGGAAATAAAGGAAAAAAAAATAAACACCCGTA‑AGGATATTTTTACAAATGTGACCA‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                         <  2:208467/146‑1 (MQ=255)
                           gaaaggccacaCCCTGCCTTTTTTATTCGTATAAAAACCCTTCTTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTC‑AGGATATTTTTATCAGTGTGACGA‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAg                                                                                                                      <  1:841283/144‑1 (MQ=255)
                           gaaaTGCCACACCCTGCCTTTATTTTTCGTAAAAAGACCCTTCGTTTGCGTACGTTTATCGAAATAAAGAAAGAAAAAATAACCACCCGTC‑AGAATATTTTTATCAAGGGGACGA‑TCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAg                                                                                                                      <  2:747010/149‑1 (MQ=255)
                           gCAAGGGCACACCCTGCCATTATTATTCTTAAAAAAAACCTTTGTTGGCGTACGTTTTTCGGAATTAGGGAAAAAGAAATTAACAACCGTA‑AGGATATTTTTACTAAT‑TGGGGATTCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                        <  1:1399229/149‑1 (MQ=255)
                            aaaTGGCACACCCTGCAAATTTTATTTGAAAAAAAAACCTTGTTTTGCTTGCGTTTATAGAAATAAAGAAAAAAGAAAATAACAACCGCC‑AGTTTGTTTTTATAATGGTGACGA‑TCCGCAACCCCGCCGATAAAGCGCGCACTAGCGAAGa                                                                                                                     <  1:686188/150‑1 (MQ=255)
                             aagggCACACCCTGCCGTTTTTATTCGTAAAAAGAACTTTCTTTTGCGTCCTTTTACCGAAATAAAGAAAAAAAAAATAAACACCCGTC‑AGTATTTTTTTACAAATGTGACGA‑CCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGat                                                                                                                    <  1:1198556/148‑1 (MQ=255)
                             aaTGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                        >  1:690444/1‑147 (MQ=255)
                                gccacaCCCTCCCGTTATTATTCGTAAAAAAACCCTTCTTTTGCGTCCGTTTATCGAAATAAAGAAAAAAAAAATTAACACCCGTC‑AGGATATTGTTATAAATGGGACGA‑TCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATc                                                                                                                 <  2:987070/149‑1 (MQ=255)
                                      aCCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                           >  2:202606/1‑151 (MQ=255)
                                          gCCGTTATTATTCGTAAAAAAAACTTTCTTTTGCTTACTTTTATCGAAATAAAGAAAGAAAAAATAAACACCCGTA‑AGGATATTGTTATCAATGTGCCGA‑CCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                                                         <  1:829634/149‑1 (MQ=255)
                                              ttattatTCTAAAAAAGACCTTTCTTTTGCGTACTTTTAGCGAAATTAAGAAAGAAAAAATAAACACCCGTA‑AGGATTTTGTTAAAAATGTGACTA‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGa                                                                                                    <  2:583961/149‑1 (MQ=255)
                                               tattatTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                                                                                  >  1:720511/1‑151 (MQ=255)
                                                 ttattCAAAAAAAGACCTTTCTTTTGGTTACTTTTACCGAAATTAAGAAAAAAAAAATAAAAACCCGCA‑AGAATATTGTTATAAAGGAGACGA‑CCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                                 <  2:297407/150‑1 (MQ=255)
                                                  tattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                                                               >  2:145261/1‑151 (MQ=255)
                                                     tCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                            >  1:229976/1‑151 (MQ=255)
                                                        aaaaaaaacccTTCTTTTGCGTACTTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGTC‑AGGATTTTTTAATAAAGGTGAGGT‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAcc                                                                                         <  1:746260/142‑1 (MQ=255)
                                                         ataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                        >  1:617746/1‑151 (MQ=255)
                                                          taAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                       >  2:246318/1‑151 (MQ=255)
                                                           aaaGACCCTTCTTTGGCGTACGTTTATGGAAATAAAGAAAAAAAAAATTAACACCCGTA‑AGGATATTGTTATCAATGTGACGA‑TCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                      <  1:106815/151‑1 (MQ=255)
                                                           aaaGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                       >  2:550154/1‑150 (MQ=255)
                                                             aaaccttttttttGCGTACGTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGCA‑AGAATATTGTTATAAATGTGACGA‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                    <  1:778209/141‑1 (MQ=255)
                                                               acccTTCTTTTGCGTACTTTTATCGAATTTAAGAAAAAAAAAATAAACACCCGCC‑AGGATTTTTTTATCAATGTGCCGT‑TCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                                  <  1:892345/149‑1 (MQ=255)
                                                                 ccttgttttGCGTAATTTTCTCAAAATAAAGAAAGAAGAAATTAACAACCGTC‑AGTTTATTATTATAAATGTGCCGA‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                                <  1:736226/145‑1 (MQ=255)
                                                                        ttGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTC‑AGGATTTTTTTATAAAGGTGACGA‑TCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                          <  1:877512/150‑1 (MQ=255)
                                                                        ttGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                          >  2:461358/1‑150 (MQ=255)
                                                                         tgggtccttttaaaagAATTAAGAAAGAAGAAATCAACACCCAGA‑GGGTTATTGTTATAAATGGGACCC‑GCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                                                                         <  2:1458177/135‑1 (MQ=255)
                                                                         tGCTTACTTTTAGCGAAATAAAGAAAAAAAAAAAAAACACCCGTG‑AGAATATTGTTACAATTGTGACGA‑TCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                                                                         <  1:580849/150‑1 (MQ=255)
                                                                          gCGTTCTTTTATGGAAATAAAGGAAAAAAAAATAAAAACCCGTC‑AGGATATTTTTATCAATGTGACAA‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                        <  2:646864/150‑1 (MQ=255)
                                                                          gCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattatt                                                                       >  2:685845/1‑151 (MQ=255)
                                                                          gCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattatt                                                                       >  2:540330/1‑151 (MQ=255)
                                                                          gCGTACGTTTATCGAAATAAAGGAAAAAAAAATAAACAACCGTC‑AGAATTTTTTTATCAATGTGACAA‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                        <  1:1091217/150‑1 (MQ=255)
                                                                           cGTACTTTTATCGAAATAAAGGAAAAAAAAATTAACACCCGTC‑AGGATATTGTTATCAATGGGACGA‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                      <  1:681015/151‑1 (MQ=255)
                                                                            gTATGTTTATCGGAATTAAGAAAAAAAAAATTACCAACCGTC‑AGAATTTTTTTATCAAGGGCACGT‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCa                                                                     <  1:570597/151‑1 (MQ=255)
                                                                            gTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                      >  2:634621/1‑150 (MQ=255)
                                                                             tactttgatggaaaTTAAGGAAGAAAAAATAACCACCCGAA‑AGGATTTTTTTCAAAGGGTGACGT‑TCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                      <  1:1147496/137‑1 (MQ=255)
                                                                              ccttttATCGAAATAAAGAAAAAAAAAATAAACACCCTCC‑AGTATTTTTTAATCAATGTGACGA‑CCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg                                                                    <  1:572665/147‑1 (MQ=255)
                                                                              acttttATGGAAATAAAGAAAGAAAAAATAAACACCCGTC‑AGAATATTTTTATAAATGTGACGA‑CCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                   <  1:1303819/148‑1 (MQ=255)
                                                                               tttttATCGAAATAAAGAAAAAAAAAATAAACAACCGTA‑AGGTTTTTGTTATAAAGGTGACGA‑CCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                   <  2:330727/148‑1 (MQ=255)
                                                                                gTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                 >  1:1025037/1‑151 (MQ=255)
                                                                                gTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                 >  2:118820/1‑151 (MQ=255)
                                                                                 tttAGCGAAATAAGGAAAAAAAAAATAAACAACCGTA‑AGGATATTGTTATCAATGGGACGA‑TCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                 <  1:730656/150‑1 (MQ=255)
                                                                                  ttATGGAAATAAAGGAAAAAAAAATAAACACCCGCA‑AGGATTTTTTTATAATTGTGACGA‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg                                                                <  1:1086470/150‑1 (MQ=255)
                                                                                  ttATGGAAATAAAGAAAAAAAAAATAACCAACCGTA‑AGGATTTTGTTATAAATGTGACCA‑CCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg                                                                <  2:1279238/150‑1 (MQ=255)
                                                                                  ttATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                               >  1:1085137/1‑151 (MQ=255)
                                                                                   attCGTAATGAAGGAAAAAAAAATGAACACCCAGA‑AGTATATTTTTATAAATTTGATGC‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                                                              <  2:1375290/149‑1 (MQ=255)
                                                                                    atggaaaTAAAGAAAAAAAAAATAAACACCCGTC‑AGGATATTTTTACAAATGTGACAC‑CCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                                                              <  1:674634/145‑1 (MQ=255)
                                                                                    aTCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                               >  2:1332920/1‑149 (MQ=255)
                                                                                    aTCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                               >  2:798481/1‑149 (MQ=255)
                                                                                     tCGAAATAAAGAAAAAAAAAATTACCACCCGTA‑AGTATTTTTTTATAAATGTGACCA‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                             <  2:810148/150‑1 (MQ=255)
                                                                                     agcaaaatAAAGAAAAAAAAAATTACCACCGTAAAGGATATTTTTATAAATGTGACGT‑TCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                                                              <  2:1104558/143‑1 (MQ=255)
                                                                                      ggaaaTAAAGAAAAAAAAAATAAACAACCGTA‑AGAATTTTTTTATAAATGGGCCAT‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTaga                                                            <  2:1120905/147‑1 (MQ=255)
                                                                                      cGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                             >  1:479862/1‑149 (MQ=255)
                                                                                       ggAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                             >  1:1449603/1‑148 (MQ=255)
                                                                                         aaTAAAGAAAAAAAAAAAAAACACCCGTC‑AGTATTTTTTTAAAAATGTGACTC‑CCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAAt                                                        <  2:171619/151‑1 (MQ=255)
                                                                                          ataaaGAAAAAAAAAATTAACACCCGTA‑AGGTTTTTGTTATAAATGTGACGA‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa                                                       <  2:169714/148‑1 (MQ=255)
                                                                                           ttAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa                                                       >  2:761609/1‑150 (MQ=255)
                                                                                            tAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                     >  2:1202713/1‑151 (MQ=255)
                                                                                            aaaGAAAAAAAAAATAAACACCCGTC‑AGGATATTTTTACAAATTTGACAC‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATaa                                                      <  2:647230/149‑1 (MQ=255)
                                                                                             aaGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                    >  1:1127021/1‑151 (MQ=255)
                                                                                                aaagaaaaaATAAACACCCGGC‑AGAATATTTTTATAAATGGGACCA‑CCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGtt                                                 <  1:1276589/150‑1 (MQ=255)
                                                                                                  aaaaaaaaTTAACAACCTCG‑AGGTTGTTTTTACAAATGTGAGTA‑CCCCCAGCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtc                                               <  2:1265544/146‑1 (MQ=255)
                                                                                                   gaagaaATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtct                                              >  2:208632/1‑151 (MQ=255)
                                                                                                       aaaTTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                          >  1:858243/1‑151 (MQ=255)
                                                                                                       aaaTTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                          >  2:314687/1‑151 (MQ=255)
                                                                                                         aataaccaccggtaaGGATTTTTTTATCAATGTGACGC‑TCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                         <  2:424806/137‑1 (MQ=255)
                                                                                                            aacaacCGTG‑AGTTTATTTTAATAATTGTGCCAA‑TCCACCACGCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                     <  1:1316521/151‑1 (MQ=255)
                                                                                                                 ccGTC‑AGGATATTTTTATAAATGTGACGA‑TCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaaa                                <  1:1335147/151‑1 (MQ=255)
                                                                                                                       aGGATATTTTTAACAATGTGCCGA‑TCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAAc                           <  1:734284/151‑1 (MQ=255)
                                                                                                                       aGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAAc                           >  1:836806/1‑151 (MQ=255)
                                                                                                                         aatttttttAACAATTTGACGA‑CCCACGCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                         <  1:897561/144‑1 (MQ=255)
                                                                                                                              tGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAt                    >  1:597178/1‑151 (MQ=255)
                                                                                                                               tttATAAATGGGACAA‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAt                    <  1:300617/149‑1 (MQ=255)
                                                                                                                               gTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtg                   >  2:1186108/1‑151 (MQ=255)
                                                                                                                                 tATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg                 >  2:320989/1‑151 (MQ=255)
                                                                                                                                  aTCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGc                >  1:988072/1‑151 (MQ=255)
                                                                                                                                  aTCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGc                >  1:993568/1‑151 (MQ=255)
                                                                                                                                      aTGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGGACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg             >  2:366764/1‑150 (MQ=255)
                                                                                                                                      aTGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg             >  1:1071165/1‑150 (MQ=255)
                                                                                                                                          gaagA‑TCCACTCCCCCGCCGATACAGCGCTCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATACCCGTCTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGCTATc        <  2:789117/148‑1 (MQ=255)
                                                                                                                                           ccGA‑CCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCg       <  2:229980/150‑1 (MQ=255)
                                                                                                                                           aCGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCg       >  2:700095/1‑151 (MQ=255)
                                                                                                                                            cGA‑TCCACTACCCCGCCGATAAAGCGCGCACTGGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc      >  2:1284348/1‑151 (MQ=255)
                                                                                                                                                 ccACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGt  >  2:219801/1‑151 (MQ=255)
                                                                                                                                                 ccACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGt  >  2:865366/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTC‑AGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGT  >  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: ‑