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

TCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGC  >  NC_000913/2421297‑2421584
                                                                                                                                                   |                                                                                                                                             
tCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAGTTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAcccc                                                                                                                                           >  2:614156/1‑151 (MQ=255)
  cTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgc                                                                                                                                         >  2:768937/1‑151 (MQ=255)
    aaGTAAAACAACGGAGACCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgccg                                                                                                                                       <  2:15446/151‑1 (MQ=255)
    aaGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACACCCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgccg                                                                                                                                       <  2:267582/151‑1 (MQ=255)
    aaGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgc                                                                                                                                         <  1:1238525/149‑1 (MQ=255)
             aacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgc                                                                                                                              >  2:1034894/1‑151 (MQ=255)
             aacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgc                                                                                                                              >  2:504333/1‑151 (MQ=255)
                   gaACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                         >  2:934578/1‑150 (MQ=255)
                    aaCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                         >  1:562316/1‑149 (MQ=255)
                           aaggccacacccgcccGTTATTATTCGAAAAAAGACCCTTCTTTTGCGTACTTTTATCGAAATTAAGA‑AAAAAAAAATAAACACCCGTCAGAATTTTTTTATCAGTGGGCCGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGat                                                                                                                <  1:1060106/137‑1 (MQ=255)
                           aaggccacaCCCTGCCGTTTTTATTCGTATAAAGACCCTTCTTTTGCGTTCTTTTATGGGAATAAAGG‑AAAAAAAAATAAACACCCGTCAGGATTTTTTTATCAATGTGACTACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                    <  1:1283315/142‑1 (MQ=255)
                             tGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                                                              >  2:225473/1‑151 (MQ=255)
                               acacaACCTCCCTTTATTATTCGAATAAAGAACTTTCGTTTGCGTTGTTTTATCGAATGAAAGA‑AAGAAAAAATCAACACCCGTAGGGATGTTTTTAACGATGGTACACCCCGCGCCCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                                                              <  2:350394/148‑1 (MQ=255)
                                 acacCCTGCCGTTATTATTGAAAAAAAGAACTTTCTTTTGCTTACTTTTATGGAAATAAAGA‑AAGAAAAAATTAACACCCGGAAGTATATTTTTACAAATGTGACTTCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGata                                                                                                               <  1:1213449/146‑1 (MQ=255)
                                  caACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                         >  2:208973/1‑151 (MQ=255)
                                    cccTCCCGTTATTTTTATTAAAAAGAACCTTCTTTTGCGTCCTTTTATCGGAATAAAGG‑AAGAAAAAATAACCACCCGTCAGGATATTTTTATCAATGTGACCATCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGa                                                                                                        <  1:693970/149‑1 (MQ=255)
                                     ccgcccgtttttattaGAAAAAAAACCCTTCTTTTGCGTCCTTTTATCGAAATAAAGA‑AAAAAAAAAAAAACACCCGAAAGGTTTTTTTTAAAAAGGTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                                                                            <  1:226060/135‑1 (MQ=255)
                                     ccTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATg                                                                                                      >  2:1096325/1‑151 (MQ=255)
                                       tcctgCTATTTTTCATAAAAAGAACCTTGTTTTGCTTACGTTTACAGAAATAAAGG‑AAAAAAAAATAAACACCCGTCAGTTTATTTTAACAAATGTGCCCCTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                                                     <  1:1348816/146‑1 (MQ=255)
                                       tGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                    >  2:851437/1‑151 (MQ=255)
                                         ccGTTATTTTTCGAAAAAAGACCCTTCTTTTGCGTACTTTTATCGAAATAAAGA‑AAAAAAAAATAAACACCCGTCAGGATATTTTTAACAATGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                  <  1:34993/151‑1 (MQ=255)
                                         ccGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                  >  2:233019/1‑151 (MQ=255)
                                          cGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                 >  2:1192544/1‑151 (MQ=255)
                                           ccTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGACTTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                 >  2:1024308/2‑150 (MQ=255)
                                               ttattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                             >  1:240095/1‑150 (MQ=255)
                                               ttattCGTAAAAAGACCTTTCTTTTGCGTACGTTTATC‑GAATTAAGGAAAAAAAAAATAAACACCCGTCGGGATTTTTTTATAATTGTGCCGACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGaa                                                                                               <  1:537462/148‑1 (MQ=255)
                                                tattgttaaaaaGACCCTTCTTTTGCGTTTGTTTATGGAAATTAAGA‑AAAAAAAAAAAACCACCCGTAAGGATATTTTTATAAAGGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                                                                            <  2:509869/141‑1 (MQ=255)
                                                tattagaaaaaaGACCCTTCGTTTGTTTTCTTGTATCGAAATAAAGG‑AAAAAAAAATAAACACCCGTAAGGATATTTTTAACAATGTGACACCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                                                                            <  2:1042574/141‑1 (MQ=255)
                                                tattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                                                           >  1:393779/1‑151 (MQ=255)
                                                  ttCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                         >  2:1042773/1‑151 (MQ=255)
                                                  ttCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                         >  1:365754/1‑151 (MQ=255)
                                                   tCGTAAAAAGACCTTTCGTTTGCGTACTTTTATCGGAATTAAGG‑AAAAAAAAATTAACACCCGTCAGGATATTTTTATCAATGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                         <  1:412516/150‑1 (MQ=255)
                                                   tCGTAAAAAGAACCTTTTTTTGCGTCCGTTTACGGAAATTAGGA‑AAAAAGAAATAAACAGCCGGCAGTATATTTTTAAAAATGTGCCGACCCCCGACGCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                                                          <  1:1114188/149‑1 (MQ=255)
                                                    cttaaaaaaacccTTCTTTTGCTTCCGTTTATCGAAATAAAGA‑AAAAAAAAATTAACACCCGTCAGGATATTTTTATAAAGGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                         <  1:834694/138‑1 (MQ=255)
                                                    cataaaaaGACCCTTCTTTTGCGTACTTTTATCGAAATAAAGA‑AAAAAAAAATAAACACCCGACAGGATATTTTTATAAATGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCa                                                                                       <  1:30478/146‑1 (MQ=255)
                                                    cGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                        >  1:1280924/1‑150 (MQ=255)
                                                    cGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                        >  2:1084039/1‑150 (MQ=255)
                                                         aaacacttttCGTTGGCGTACGTTTTTGGGAATAAAGG‑AAGAAGAAAATAACACCGGTAAGGATTTTGTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                  <  2:589193/143‑1 (MQ=255)
                                                             acccTTCGTTTGCTTTGTTTTATCGAAATAAAAA‑AAAAAAAAATAACCACCTGTCAGGTTTTTTTTACAAATGTGACCATCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                               <  2:997249/148‑1 (MQ=255)
                                                             aaCCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                               >  1:535139/1‑150 (MQ=255)
                                                             aaCCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCAACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                               >  2:253408/1‑150 (MQ=255)
                                                                cTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                            >  2:1302651/1‑150 (MQ=255)
                                                                   cGTTTGCGTACGTTTATCGGAATTAAGG‑AAGCAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                         >  2:11743/1‑150 (MQ=255)
                                                                   cGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                         >  1:997582/1‑150 (MQ=255)
                                                                   cGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                         >  1:834281/1‑150 (MQ=255)
                                                                     tttGTGTATTTTTAT‑GGAAATAAGGAAAAAAAAAATAAACCCCGGTCAGAATTTTTTTATCAGTGAGCCACTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                              <  1:1175010/143‑1 (MQ=255)
                                                                     tttGCGTACTTTTATAGAAATAAAGA‑AAAAAAAAATAAAAACCCGCAAGGATTTTGTTAACAATGTGACTACCCTCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaaa                                                                          <  1:731263/147‑1 (MQ=255)
                                                                     tttGCGTACTTTTAGCGAAATAAAGA‑AAAAAAAAATAAACACCGGTCAGAATTTTGTTATCAATGTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                           <  2:1144910/146‑1 (MQ=255)
                                                                     tttGCGTACGTTTATCGAAATAAAGA‑AAGAAAAAATTAACACCCGTCAGAATATTTTTATCAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                           <  1:501069/146‑1 (MQ=255)
                                                                     ttggCGTACGTTTATGGAAATTAAGG‑AAGAAAAAATTAACACCGGAGGGTATTTGTTTAAAAAGGTGACGCTCCACCACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACa                                                                        <  2:975363/146‑1 (MQ=255)
                                                                       tGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                    >  1:674059/1‑151 (MQ=255)
                                                                               tttATGAGAATAAAGA‑AAGAAAAAATAAACACCCGTCAGTATATTGTTAACAAGGTGACTACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCa                                                                 <  1:1088566/146‑1 (MQ=255)
                                                                                ttATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                           >  2:751130/1‑151 (MQ=255)
                                                                                    gggAATAAAGA‑AAAAAAAAATAACCACCCGAGAGAATTTTTTTATAAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag                                                       <  2:426147/150‑1 (MQ=255)
                                                                                    cgaaaTTAAAA‑AAGAAAAAATAAACACCGGTGAGGTTATTTTTATCATTGTGTCCTTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                             <  2:789936/142‑1 (MQ=255)
                                                                                    cGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                         >  1:989207/1‑149 (MQ=255)
                                                                                       aaTAAAGA‑AAAAAAAAATAAACACCCGTAAGGATATTTTTATAAATGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag                                                       <  1:769495/148‑1 (MQ=255)
                                                                                        ataaaGA‑AAAAAAAAATAAACACCCGAAAGGATTTTTTTATAAATGTGACGCCCCCCACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa                                                     <  1:67007/146‑1 (MQ=255)
                                                                                        aataaaaaaagaaaaaATAAACACCCGTCAGGATTTTGTAATAAATGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                                                          <  2:1238221/137‑1 (MQ=255)
                                                                                            aGG‑AAGAAAAAATTAACACCCGTCAGGATTTTTTTATCAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCg                                               <  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)
                                                                                                                                                   |                                                                                                                                             
TCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGC  >  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: ‑