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

TCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGC  >  NC_000913/2421297‑2421584
                                                                                                                                                   |                                                                                                                                             
tCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAGTTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAcccc                                                                                                                                           >  2:614156/1‑151 (MQ=255)
  cTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgc                                                                                                                                         >  2:768937/1‑151 (MQ=255)
    aaGTAAAACAACGGAGACCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgccg                                                                                                                                       <  2:15446/151‑1 (MQ=255)
    aaGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACACCCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgccg                                                                                                                                       <  2:267582/151‑1 (MQ=255)
    aaGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgc                                                                                                                                         <  1:1238525/149‑1 (MQ=255)
             aacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgc                                                                                                                              >  2:1034894/1‑151 (MQ=255)
             aacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgc                                                                                                                              >  2:504333/1‑151 (MQ=255)
                   gaACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                         >  2:934578/1‑150 (MQ=255)
                    aaCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                         >  1:562316/1‑149 (MQ=255)
                           aaggccacacccgcccGTTATTATTCGAAAAAAGACCCTTCTTTTGCGTACTTTTATCGAAATTAAGAAAAAAAAAATAAACACCCGTCAGAATTTTTTTATCAGTGGGCCGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGat                                                                                                                <  1:1060106/137‑1 (MQ=255)
                           aaggccacaCCCTGCCGTTTTTATTCGTATAAAGACCCTTCTTTTGCGTTCTTTTATGGGAATAAAGGAAAAAAAAATAAACACCCGTCAGGATTTTTTTATCAATGTGACTACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                    <  1:1283315/142‑1 (MQ=255)
                             tGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                                                              >  2:225473/1‑151 (MQ=255)
                               acacaACCTCCCTTTATTATTCGAATAAAGAACTTTCGTTTGCGTTGTTTTATCGAATGAAAGAAAGAAAAAATCAACACCCGTAGGGATGTTTTTAACGATGGTACACCCCGCGCCCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                                                              <  2:350394/148‑1 (MQ=255)
                                 acacCCTGCCGTTATTATTGAAAAAAAGAACTTTCTTTTGCTTACTTTTATGGAAATAAAGAAAGAAAAAATTAACACCCGGAAGTATATTTTTACAAATGTGACTTCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGata                                                                                                               <  1:1213449/146‑1 (MQ=255)
                                  caACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                         >  2:208973/1‑151 (MQ=255)
                                    cccTCCCGTTATTTTTATTAAAAAGAACCTTCTTTTGCGTCCTTTTATCGGAATAAAGGAAGAAAAAATAACCACCCGTCAGGATATTTTTATCAATGTGACCATCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGa                                                                                                        <  1:693970/149‑1 (MQ=255)
                                     ccgcccgtttttattaGAAAAAAAACCCTTCTTTTGCGTCCTTTTATCGAAATAAAGAAAAAAAAAAAAAACACCCGAAAGGTTTTTTTTAAAAAGGTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                                                                            <  1:226060/135‑1 (MQ=255)
                                     ccTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATg                                                                                                      >  2:1096325/1‑151 (MQ=255)
                                       tcctgCTATTTTTCATAAAAAGAACCTTGTTTTGCTTACGTTTACAGAAATAAAGGAAAAAAAAATAAACACCCGTCAGTTTATTTTAACAAATGTGCCCCTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                                                     <  1:1348816/146‑1 (MQ=255)
                                       tGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                    >  2:851437/1‑151 (MQ=255)
                                         ccGTTATTTTTCGAAAAAAGACCCTTCTTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTCAGGATATTTTTAACAATGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                  <  1:34993/151‑1 (MQ=255)
                                         ccGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                  >  2:233019/1‑151 (MQ=255)
                                          cGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                 >  2:1192544/1‑151 (MQ=255)
                                           ccTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGACTTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                 >  2:1024308/2‑150 (MQ=255)
                                               ttattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                             >  1:240095/1‑150 (MQ=255)
                                               ttattCGTAAAAAGACCTTTCTTTTGCGTACGTTTATCGAATTAAGGAAAAAAAAAATAAACACCCGTCGGGATTTTTTTATAATTGTGCCGACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGaa                                                                                               <  1:537462/148‑1 (MQ=255)
                                                tattgttaaaaaGACCCTTCTTTTGCGTTTGTTTATGGAAATTAAGAAAAAAAAAAAAACCACCCGTAAGGATATTTTTATAAAGGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                                                                            <  2:509869/141‑1 (MQ=255)
                                                tattagaaaaaaGACCCTTCGTTTGTTTTCTTGTATCGAAATAAAGGAAAAAAAAATAAACACCCGTAAGGATATTTTTAACAATGTGACACCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                                                                            <  2:1042574/141‑1 (MQ=255)
                                                tattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                                                           >  1:393779/1‑151 (MQ=255)
                                                  ttCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                         >  2:1042773/1‑151 (MQ=255)
                                                  ttCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                         >  1:365754/1‑151 (MQ=255)
                                                   tCGTAAAAAGACCTTTCGTTTGCGTACTTTTATCGGAATTAAGGAAAAAAAAATTAACACCCGTCAGGATATTTTTATCAATGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                         <  1:412516/150‑1 (MQ=255)
                                                   tCGTAAAAAGAACCTTTTTTTGCGTCCGTTTACGGAAATTAGGAAAAAAGAAATAAACAGCCGGCAGTATATTTTTAAAAATGTGCCGACCCCCGACGCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                                                          <  1:1114188/149‑1 (MQ=255)
                                                    cttaaaaaaacccTTCTTTTGCTTCCGTTTATCGAAATAAAGAAAAAAAAAATTAACACCCGTCAGGATATTTTTATAAAGGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                         <  1:834694/138‑1 (MQ=255)
                                                    cataaaaaGACCCTTCTTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGACAGGATATTTTTATAAATGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCa                                                                                       <  1:30478/146‑1 (MQ=255)
                                                    cGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                        >  1:1280924/1‑150 (MQ=255)
                                                    cGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                        >  2:1084039/1‑150 (MQ=255)
                                                         aaacacttttCGTTGGCGTACGTTTTTGGGAATAAAGGAAGAAGAAAATAACACCGGTAAGGATTTTGTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                  <  2:589193/143‑1 (MQ=255)
                                                             acccTTCGTTTGCTTTGTTTTATCGAAATAAAAAAAAAAAAAATAACCACCTGTCAGGTTTTTTTTACAAATGTGACCATCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                               <  2:997249/148‑1 (MQ=255)
                                                             aaCCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                               >  1:535139/1‑150 (MQ=255)
                                                             aaCCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCAACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                               >  2:253408/1‑150 (MQ=255)
                                                                cTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                            >  2:1302651/1‑150 (MQ=255)
                                                                   cGTTTGCGTACGTTTATCGGAATTAAGGAAGCAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                         >  2:11743/1‑150 (MQ=255)
                                                                   cGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                         >  1:997582/1‑150 (MQ=255)
                                                                   cGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                         >  1:834281/1‑150 (MQ=255)
                                                                     tttGTGTATTTTTATGGAAATAAGGAAAAAAAAAATAAACCCCGGTCAGAATTTTTTTATCAGTGAGCCACTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                              <  1:1175010/143‑1 (MQ=255)
                                                                     tttGCGTACTTTTATAGAAATAAAGAAAAAAAAAATAAAAACCCGCAAGGATTTTGTTAACAATGTGACTACCCTCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaaa                                                                          <  1:731263/147‑1 (MQ=255)
                                                                     tttGCGTACTTTTAGCGAAATAAAGAAAAAAAAAATAAACACCGGTCAGAATTTTGTTATCAATGTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                           <  2:1144910/146‑1 (MQ=255)
                                                                     tttGCGTACGTTTATCGAAATAAAGAAAGAAAAAATTAACACCCGTCAGAATATTTTTATCAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                           <  1:501069/146‑1 (MQ=255)
                                                                     ttggCGTACGTTTATGGAAATTAAGGAAGAAAAAATTAACACCGGAGGGTATTTGTTTAAAAAGGTGACGCTCCACCACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACa                                                                        <  2:975363/146‑1 (MQ=255)
                                                                       tGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                    >  1:674059/1‑151 (MQ=255)
                                                                               tttATGAGAATAAAGAAAGAAAAAATAAACACCCGTCAGTATATTGTTAACAAGGTGACTACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCa                                                                 <  1:1088566/146‑1 (MQ=255)
                                                                                ttATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                           >  2:751130/1‑151 (MQ=255)
                                                                                    gggAATAAAGAAAAAAAAAATAACCACCCGAGAGAATTTTTTTATAAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag                                                       <  2:426147/150‑1 (MQ=255)
                                                                                    cgaaaTTAAAAAAGAAAAAATAAACACCGGTGAGGTTATTTTTATCATTGTGTCCTTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                             <  2:789936/142‑1 (MQ=255)
                                                                                    cGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                         >  1:989207/1‑149 (MQ=255)
                                                                                       aaTAAAGAAAAAAAAAATAAACACCCGTAAGGATATTTTTATAAATGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag                                                       <  1:769495/148‑1 (MQ=255)
                                                                                        ataaaGAAAAAAAAAATAAACACCCGAAAGGATTTTTTTATAAATGTGACGCCCCCCACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa                                                     <  1:67007/146‑1 (MQ=255)
                                                                                        aataaaaaaagaaaaaATAAACACCCGTCAGGATTTTGTAATAAATGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                                                          <  2:1238221/137‑1 (MQ=255)
                                                                                            aGGAAGAAAAAATTAACACCCGTCAGGATTTTTTTATCAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCg                                               <  1:208706/151‑1 (MQ=255)
                                                                                                aaaaaaaaaacaacaccCGTCAAGATATTTTTATAAAGTTGACCATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                            <  2:314321/140‑1 (MQ=255)
                                                                                                 agaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtc                                           >  1:292122/1‑151 (MQ=255)
                                                                                                  aaaaaaaTAAACACCCGTCAGAATATTATTATCAATGTGACCCTCCACACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtct                                          <  2:284934/147‑1 (MQ=255)
                                                                                                   aaaaaaTTAACACCCGTCAGAATATTTTTATAAAGGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtct                                          <  1:504448/147‑1 (MQ=255)
                                                                                                     aaaaTAAACCCCCGTAAGATTATTTTTAACAATGTGACCCCCCCTCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                       <  2:572527/150‑1 (MQ=255)
                                                                                                      aaaTTAACACCCATAAGTATTTTTTTAACAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                      <  2:55994/151‑1 (MQ=255)
                                                                                                          tAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                  >  1:590555/1‑151 (MQ=255)
                                                                                                          tAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACACCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGCTGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                  >  2:870232/1‑151 (MQ=255)
                                                                                                           accacccGTCAGGATTTTTTTATAAATGTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                  <  1:935203/145‑1 (MQ=255)
                                                                                                            acaacCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                 >  1:1300845/1‑150 (MQ=255)
                                                                                                             cacccGTCGGTATTTTTTTACAAATGTGCCGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGa                                <  1:927797/147‑1 (MQ=255)
                                                                                                             caacCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGa                                >  2:279813/1‑150 (MQ=255)
                                                                                                             caacCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGa                                >  1:630632/1‑150 (MQ=255)
                                                                                                                ccGTCAGGATATTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaa                             <  1:948399/150‑1 (MQ=255)
                                                                                                                 cGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAt                           >  1:337349/1‑151 (MQ=255)
                                                                                                                  gacaagATATTGTAAACAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaaa                            <  2:838159/144‑1 (MQ=255)
                                                                                                                  gTCAGGATATTGTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAtt                          <  1:233190/151‑1 (MQ=255)
                                                                                                                    cAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTa                         >  1:50030/1‑150 (MQ=255)
                                                                                                                      ggATTTTTTTTTAAATGTGCCGTCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAtt                          <  1:109342/147‑1 (MQ=255)
                                                                                                                        atatTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCCCTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAaactaact                      >  2:1129870/1‑149 (MQ=255)
                                                                                                                        atatTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                     >  1:1173192/1‑150 (MQ=255)
                                                                                                                         ttttGTTATAAATGTTACGACCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCa                    <  2:1190640/148‑1 (MQ=255)
                                                                                                                         tatTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                     >  2:863325/1‑149 (MQ=255)
                                                                                                                          atTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAgg                  >  1:323876/1‑151 (MQ=255)
                                                                                                                            tGTTAGAAAGGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAgg                  <  2:505447/149‑1 (MQ=255)
                                                                                                                             gTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtg               >  1:383203/1‑151 (MQ=255)
                                                                                                                                   aaggtgACGACCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg             <  2:314782/144‑1 (MQ=255)
                                                                                                                                   aaTGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg         >  2:890290/1‑151 (MQ=255)
                                                                                                                                       tgACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGatat     >  1:939806/1‑151 (MQ=255)
                                                                                                                                          cgctCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc  <  2:1100586/148‑1 (MQ=255)
                                                                                                                                            acccACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACACTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGatat     <  2:562189/144‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGC  >  NC_000913/2421297‑2421584

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