Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 710,173 A→G 6.2% intergenic (+57/+27) chiQ → / ← fur chitosugar‑induced verified lipoprotein/ferric iron uptake regulon transcriptional repressor; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913710,1730AG 6.2% 867.8 / 30.1 258intergenic (+57/+27)chiQ/furchitosugar‑induced verified lipoprotein/ferric iron uptake regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (116/126);  new base G (9/7);  total (125/133)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 260 total aligned reads displayed.

AGTGTCTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGCTGACGCGC  >  NC_000913/710025‑710322
                                                                                                                                                    |                                                                                                                                                     
    tCTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTg                                                                                                                                                  >  2:337992/1‑150 (MQ=255)
     cTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTgg                                                                                                                                                 <  1:101032/150‑1 (MQ=255)
     cTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTgg                                                                                                                                                 <  2:3843958/150‑1 (MQ=255)
     cTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTgg                                                                                                                                                 >  1:223590/1‑150 (MQ=255)
      tGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGc                                                                                                                                                >  1:770761/1‑150 (MQ=255)
      tGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGc                                                                                                                                                >  1:3711642/1‑150 (MQ=255)
      tGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGc                                                                                                                                                >  1:1558753/1‑150 (MQ=255)
       gCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCt                                                                                                                                               <  1:2217999/150‑1 (MQ=255)
        cGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGGCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGTGCTGTAACTCACGCTTTTCTTATTTCCCTTGCATAAAAAAGCCTTCACGCAGGTTGGCtt                                                                                                                                              >  1:3967347/1‑150 (MQ=255)
          tGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCtttt                                                                                                                                            <  2:1114602/150‑1 (MQ=255)
          tGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCtttt                                                                                                                                            <  2:3414688/150‑1 (MQ=255)
           gCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtc                                                                                                                                           >  1:4059/1‑150 (MQ=255)
             aaCGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCATTAATAACAAATACGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtctc                                                                                                                                         <  1:2052138/150‑1 (MQ=255)
             aaCGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtctc                                                                                                                                         <  1:605345/150‑1 (MQ=255)
             aaCGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtctc                                                                                                                                         >  2:455637/1‑150 (MQ=255)
               cGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATATGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCGTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGt                                                                                                                                       <  1:3497676/150‑1 (MQ=255)
                 cAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTc                                                                                                                                     <  1:1046866/150‑1 (MQ=255)
                      cGGAAATTATCAGGCGGTGACAGCCGAGTGTGATAATGTCTGGCAGGAAATACGCAGTAAAAACAAATAACTGCGAGCTGTAACTCTCGCTTTTCTTAGTTCCCTGGCATAAAATAGCCTAGCCGCAGGTTGGCTTTTCTCGTTCAggct                                                                                                                                <  1:4227064/150‑1 (MQ=255)
                      cGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggct                                                                                                                                >  1:3633792/1‑150 (MQ=255)
                      cGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggct                                                                                                                                >  2:3539609/1‑150 (MQ=255)
                       ggAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggctg                                                                                                                               <  2:1186888/150‑1 (MQ=255)
                          aaTGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggctggct                                                                                                                            >  2:1424848/1‑150 (MQ=255)
                                cAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTg                                                                                                                      >  2:2269027/1‑150 (MQ=255)
                                 aGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGc                                                                                                                     >  2:782847/1‑150 (MQ=255)
                                       tGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGt                                                                                                               <  1:1789057/150‑1 (MQ=255)
                                           aGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTgcgc                                                                                                           >  1:841198/1‑150 (MQ=255)
                                             ccGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCAt                                                                                                         <  1:798666/150‑1 (MQ=255)
                                                 ttGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTc                                                                                                     <  1:321284/150‑1 (MQ=255)
                                                   gtgATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATACGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCAt                                                                                                   <  2:3382084/150‑1 (MQ=255)
                                                   gtgATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCAt                                                                                                   >  1:3074283/1‑150 (MQ=255)
                                                   gtgATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCAt                                                                                                   <  1:4329077/150‑1 (MQ=255)
                                                    tgATAAAGTCTGTCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGTTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATc                                                                                                  <  1:4256165/150‑1 (MQ=255)
                                                    tgATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATc                                                                                                  >  1:692331/1‑150 (MQ=255)
                                                    tgATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATc                                                                                                  >  1:1954136/1‑150 (MQ=255)
                                                     gATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCt                                                                                                 <  2:1140001/150‑1 (MQ=255)
                                                      aTAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTAGTTGCCTTCGTGCGCATGTTCATCtt                                                                                                >  2:3469208/1‑150 (MQ=255)
                                                       tAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTc                                                                                               <  1:3917354/150‑1 (MQ=255)
                                                          aGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTcgcg                                                                                            <  1:989494/150‑1 (MQ=255)
                                                          aGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTcgcg                                                                                            <  1:1982698/150‑1 (MQ=255)
                                                            tCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGc                                                                                          >  1:1709623/1‑150 (MQ=255)
                                                                   ggAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGcc                                                                                   <  1:1004603/150‑1 (MQ=255)
                                                                    gAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTGATTTCCCTGGCATAAAAAAGCCAACCCGCGGGGTGGCTTTTCTCGTTCAGGCTGGATTATGTGCGTTCGTGCGCATGTTCAACTTCGCGGCACTCTCCt                                                                                  >  1:1815403/1‑150 (MQ=255)
                                                                      aaTACGCAGTAATACCAAATAAGTGAGAGCTGTAACACTCGATTTTATTATTTCCCTTGCATAAAAAAGCCAACCCACAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTc                                                                                <  1:2601283/150‑1 (MQ=255)
                                                                        tACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCgg                                                                              <  2:657216/150‑1 (MQ=255)
                                                                          cGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGca                                                                            <  1:2366334/150‑1 (MQ=255)
                                                                           gCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGcac                                                                           >  1:3531886/1‑150 (MQ=255)
                                                                           gCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGcac                                                                           >  2:3936685/1‑150 (MQ=255)
                                                                                 ataaCAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGAc                                                                     <  1:920080/150‑1 (MQ=255)
                                                                                     cAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTa                                                                 >  2:21596/1‑150 (MQ=255)
                                                                                      aaaGAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTaa                                                                <  1:3982421/150‑1 (MQ=255)
                                                                                         tAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGa                                                             <  2:4251738/150‑1 (MQ=255)
                                                                                         tAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGa                                                             <  2:223571/150‑1 (MQ=255)
                                                                                               agagCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAgaga                                                       >  1:3153763/1‑150 (MQ=255)
                                                                                                   cTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACtgt                                                   >  2:1741693/1‑150 (MQ=255)
                                                                                                   cTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACtgt                                                   >  2:919060/1‑150 (MQ=255)
                                                                                                   cTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCACATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACtgt                                                   <  2:207631/150‑1 (MQ=255)
                                                                                                   cTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGACTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACtgt                                                   >  2:2638223/1‑150 (MQ=255)
                                                                                                     gTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTgg                                                 >  2:3349942/1‑150 (MQ=255)
                                                                                                      tAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGt                                                <  2:82373/150‑1 (MQ=255)
                                                                                                      tAACTCTCGCTTTTCTTAATTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCATATTTGCCTTAGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGt                                                >  1:3489451/1‑150 (MQ=255)
                                                                                                        aCTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGTCTGTGGTTa                                              <  2:3074283/150‑1 (MQ=255)
                                                                                                        aCTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTACGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTa                                              <  1:547876/150‑1 (MQ=255)
                                                                                                          tctcGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGt                                            >  1:4157189/1‑150 (MQ=255)
                                                                                                            tcGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCa                                          <  1:89862/150‑1 (MQ=255)
                                                                                                              gCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTATCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAgg                                        >  1:1339226/1‑150 (MQ=255)
                                                                                                               cTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGc                                       <  2:4211354/150‑1 (MQ=255)
                                                                                                                 tttCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGa                                     >  2:28395/1‑150 (MQ=255)
                                                                                                                   tCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAAt                                   <  2:2990146/150‑1 (MQ=255)
                                                                                                                   tCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAAt                                   <  1:663573/150‑1 (MQ=255)
                                                                                                                     ttATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGTTAGAGACTGTGGTTAGTCAGGCGAATGc                                 >  1:257547/1‑150 (MQ=255)
                                                                                                                     ttATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGc                                 >  2:3949767/1‑150 (MQ=255)
                                                                                                                     ttATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTACGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGc                                 >  1:1417477/1‑150 (MQ=255)
                                                                                                                      tATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGt                                                >  2:2559794/1‑134 (MQ=255)
                                                                                                                        tttCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTTGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCAt                              <  1:3998351/150‑1 (MQ=255)
                                                                                                                            ccTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGttt                          <  2:1719855/150‑1 (MQ=255)
                                                                                                                            ccTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGttt                          >  1:1186877/1‑150 (MQ=255)
                                                                                                                            ccTTGCATAAAAAAGCCAACCCGCAGGGTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGGCCAGATAGAGACGGTGGTTAGTCAGGCGACTGCCATGttt                          >  1:2328524/1‑150 (MQ=255)
                                                                                                                              ttGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTGTTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTg                        >  1:436297/1‑150 (MQ=255)
                                                                                                                              ttGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTg                        >  2:2749217/1‑150 (MQ=255)
                                                                                                                              ttGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTg                        >  2:1528563/1‑150 (MQ=255)
                                                                                                                                  aTAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGc                    >  2:610014/1‑150 (MQ=255)
                                                                                                                                    aaaaaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCaa                  <  1:3787318/150‑1 (MQ=255)
                                                                                                                                     gaaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCact                 >  2:568007/2‑148 (MQ=255)
                                                                                                                                     gaaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAAt                 >  1:2274038/2‑150 (MQ=255)
                                                                                                                                     gaaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAAt                 >  1:2916825/2‑150 (MQ=255)
                                                                                                                                     aaaaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAAt                 >  2:2703759/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAAtt                >  1:1140708/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAAtt                <  1:279918/150‑1 (MQ=255)
                                                                                                                                       aaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAAttt               <  2:1369088/150‑1 (MQ=255)
                                                                                                                                       aaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAAttt               <  2:1365774/150‑1 (MQ=255)
                                                                                                                                         aGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCa             >  2:4098134/1‑150 (MQ=255)
                                                                                                                                         aGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCa             <  2:432756/150‑1 (MQ=255)
                                                                                                                                           ccAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACg           >  1:2696457/1‑150 (MQ=255)
                                                                                                                                           ccAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGGAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACg           >  1:1241353/1‑150 (MQ=255)
                                                                                                                                            cAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGc          >  2:1820414/1‑150 (MQ=255)
                                                                                                                                            cAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGc          >  1:637597/1‑150 (MQ=255)
                                                                                                                                             aaCCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGCt         >  2:1299075/1‑150 (MQ=255)
                                                                                                                                               cccGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGCTGa       <  1:3052375/150‑1 (MQ=255)
                                                                                                                                                ccGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGCTGAc      >  2:1561743/1‑150 (MQ=255)
                                                                                                                                                ccGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGCTGAc      >  2:2512490/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGTGTCTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGCTGACGCGC  >  NC_000913/710025‑710322

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: