Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,239,610 C→G 5.7% intergenic (+65/‑334) alx → / → sstT putative membrane‑bound redox modulator/sodium:serine/threonine symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,239,6100CG 5.7% 903.3 / 39.3 263intergenic (+65/‑334)alx/sstTputative membrane‑bound redox modulator/sodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base C (160/87);  new base G (4/11);  total (165/98)
Fisher's exact test for biased strand distribution p-value = 4.79e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

Only 100 of 264 total aligned reads displayed.

GCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTT  >  NC_000913/3239463‑3239759
                                                                                                                                                   |                                                                                                                                                     
  gtggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                                                                                                                   >  2:222002/1‑150 (MQ=255)
  gtggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                                                                                                                   >  2:219535/1‑150 (MQ=255)
  gtggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                                                                                                                   <  1:2755985/150‑1 (MQ=255)
   tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCGTAAGCCGTGTACCCGGGCATTAACATGAGCACCGTTGTCTCCCTCTCCCTTCCTggg                                                                                                                                                  >  2:3035216/1‑150 (MQ=255)
   tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                  >  1:135082/1‑150 (MQ=255)
    ggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                                                                                                 >  1:1712916/1‑150 (MQ=255)
    ggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                                                                                                 <  1:1832992/150‑1 (MQ=255)
    ggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTATCTCCCTCTCCCTTCCAGGga                                                                                                                                                 >  2:2076352/1‑150 (MQ=255)
     gtgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                                                                                                >  1:4920452/1‑150 (MQ=255)
     gtgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                                                                                                >  2:5920851/1‑150 (MQ=255)
      tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAATCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                               >  1:3481040/1‑150 (MQ=255)
       gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGagag                                                                                                                                              <  1:3612044/150‑1 (MQ=255)
        tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                             >  2:1807523/1‑150 (MQ=255)
          tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGt                                                                                                                                           >  1:2139233/1‑150 (MQ=255)
          tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGt                                                                                                                                           >  1:3711318/1‑150 (MQ=255)
          tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAggga                                                                                                                                           <  1:4038401/150‑2 (MQ=255)
          tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAggga                                                                                                                                           <  1:1201959/150‑2 (MQ=255)
             cATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                                        >  2:1688871/1‑150 (MQ=255)
              tttCTGGTTATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCTGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCggg                                                                                                                                       <  1:4174460/149‑1 (MQ=255)
              aTTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCggg                                                                                                                                       >  2:4570524/1‑150 (MQ=255)
              aTTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATATATAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCggg                                                                                                                                       <  2:5183930/150‑1 (MQ=255)
               ttCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                                                                                                                      >  1:3172277/1‑150 (MQ=255)
                  tGGTGATGACGTTTATTATCAACGCCTGGGTTAATTATCGGCTTGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCTTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGa                                                                                                                                   <  1:3895424/150‑1 (MQ=255)
                  tGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                   <  2:864389/150‑1 (MQ=255)
                    gTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGTACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAgg                                                                                                                                 <  2:4802278/150‑1 (MQ=255)
                    gTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAgg                                                                                                                                 <  2:5233339/150‑1 (MQ=255)
                          aCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtt                                                                                                                           >  2:2759673/1‑150 (MQ=255)
                          aCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtt                                                                                                                           >  1:4346873/1‑150 (MQ=255)
                             tttATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttttt                                                                                                                        >  1:4757438/1‑150 (MQ=255)
                              ttattaTCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTc                                                                                                                       >  1:2325718/1‑150 (MQ=255)
                                  tatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAc                                                                                                                   >  1:299222/1‑150 (MQ=255)
                                  tatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAc                                                                                                                   >  1:761027/1‑150 (MQ=255)
                                   atCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAcc                                                                                                                  >  2:698518/1‑150 (MQ=255)
                                    tCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCg                                                                                                                 >  2:5544914/1‑150 (MQ=255)
                                     cAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGa                                                                                                                >  2:342196/1‑150 (MQ=255)
                                            tGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGc                                                                                                         >  1:472763/1‑150 (MQ=255)
                                            ggggTGAATTATCGGATTTAAACGCCGGGGGGTGTATATTTTTAAATCTGCCAACGCCGTTTACCCTGTCATTAACATGCCCACGGTTTTCGCCCTCCCCCTCCCTGGGAGCGGGTGGGGGTGAGGGAAATTTTTCGCACCGATGCTGGc                                                                                                         <  2:4867159/149‑1 (MQ=255)
                                              ggTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGTGGTGAGGGTAATTTTTCGCACCGATGCTGGCCt                                                                                                       >  1:1049952/1‑150 (MQ=255)
                                                 taataatCCGAATTATAACCAGCGGGTTTTATAATTTTTAATCTCAAAAAGACTTTTACCCTTTAATTAAAATTAGCACCTTTTACCCCAACCCCCTTCAAGTTAGAGGGTCGGGGTGAGGTTAAATTTTCGCACCGATGCTGGCCTGtt                                                                                                    <  2:2868925/149‑1 (MQ=255)
                                                  atttttCGGCATTATAAGCAGCGGGGTGCTAAATTTTTACTCTGCCTACGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCACCCTCTCCCTACCTGGTAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTc                                                                                                   <  2:3554141/145‑1 (MQ=255)
                                                  aaTTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTc                                                                                                   >  2:4299076/1‑150 (MQ=255)
                                                  aaTTAGCGGCATGATAAGCAGCGGGGTGGATAATTTTTATTCTGCCTAAGCCGTGTCCCCTGTCATTAACATGCGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTc                                                                                                   <  2:4989709/150‑1 (MQ=255)
                                                    ttATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTACTCCCTCTCCCTTCCCGGGTGAGGGTGGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTccc                                                                                                 <  2:3441359/150‑1 (MQ=255)
                                                    ttATCGGAATGATAAGCAGCGGGGTGTATAATTTTTATTCGGCCTAAGCCGTGTACCCTTACATTAACATGACCACCGTTTACTCCCTCTCCCTCCCCGGGAGTGGGTCGGGGTTAGGGTAATTTTTCGCACCGATGCTGGCCTGTTccc                                                                                                 <  2:472528/150‑1 (MQ=255)
                                                       tCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTc                                                                                              >  2:5824491/1‑150 (MQ=255)
                                                               aTAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAc                                                                                      >  2:3223981/1‑150 (MQ=255)
                                                               aTAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAc                                                                                      >  1:3654570/1‑150 (MQ=255)
                                                                tAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAcc                                                                                     >  2:1314075/1‑150 (MQ=255)
                                                                     cggggggTGTATATTTTTAAATCTGCATAAGCCGTTTACCCTGTCATTAAATTGCGCCCGTTTTTCTCCCTCTCCATCCCCGGTATAGGTTGGGGTTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                                <  2:2374006/147‑1 (MQ=255)
                                                                     agcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                                >  1:2238887/1‑150 (MQ=255)
                                                                     agcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCactc                                                                                >  1:2686764/1‑150 (MQ=255)
                                                                        ggggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTcccc                                                                             >  1:4395883/1‑150 (MQ=255)
                                                                        ggggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCGGTTCCCCTCACCCTAACCCACTcccc                                                                             >  1:4613665/1‑150 (MQ=255)
                                                                               atatTTTTTAATCTGCCTATGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGTGGGTCGGGGTGAGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCGccccgtttggg                                                                      <  1:5813048/150‑8 (MQ=37)
                                                                               atatTTTTTAATCTGCCTACGCGGTGTACCCTTTAATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCTGGTTGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACggg                                                                      <  1:4426688/150‑1 (MQ=25)
                                                                               aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCGCACCCTCACCCTCTCCCCAAACggg                                                                      >  2:654044/1‑150 (MQ=255)
                                                                               aTAATTTTTAATCGCCCTAACCGTTGTACCCTGTCATTAAATTGAGCACCGTTTTCTCCCCCCCCCTCCCTGGTCGTGGGTGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCGccccgtttggg                                                                      <  1:4300420/150‑8 (MQ=19)
                                                                                tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGGTGGCGTGTTCCCCTCACCCTAACCCTCTCCCCAAACgggg                                                                     >  1:2798402/1‑150 (MQ=255)
                                                                                    ttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAg                                                                 >  1:1449171/1‑150 (MQ=255)
                                                                                     tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAgg                                                                >  1:1304525/1‑150 (MQ=255)
                                                                                      ttAACCTGCCTACGCGGTGTCCCCTGTAATTAACTTGAGCACCTTTTTCTCCCCCTCCCTTCCGGGGCGTGGGTGGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAggg                                                               <  1:618158/150‑1 (MQ=25)
                                                                                          tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACt                                                           >  1:5908729/1‑150 (MQ=255)
                                                                                            tGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGa                                                         >  2:5650559/1‑150 (MQ=255)
                                                                                             gcccacgCCGTTTCCCCTGCCATAAACAGTACCACCTTTTCCTCCCCCCCCCTTCCGGGTTGAGGGCCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGAc                                                        <  1:5654839/144‑1 (MQ=17)
                                                                                              ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACGGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGTCTGAcc                                                       >  1:4586040/1‑150 (MQ=255)
                                                                                                 aaGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAg                                                    >  2:5435873/1‑150 (MQ=255)
                                                                                                    ccGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTcc                                                 >  2:2985408/1‑150 (MQ=255)
                                                                                                    ccGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTACAGGGAGAGGGTCGGGGTGAGGGTAATTTTTAGCACCGATGCTGGCCTGTTCCACACACCCTAAACCACTCCCCAAACGGGGCGAGTGGACTGACCGAGTca                                                 >  1:2904681/1‑149 (MQ=255)
                                                                                                         tACCGTGACATTACCATGACCACCGTTTTCCCCCTCCCCCTTCCGGGGAGAGGGTGGGGTTGTGGTTCATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCttttt                                            <  1:4512881/150‑1 (MQ=21)
                                                                                                         tACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCttttt                                            <  2:3484574/150‑1 (MQ=255)
                                                                                                         tACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCttttt                                            >  2:1836395/1‑150 (MQ=255)
                                                                                                          ccccTGTATTCAAAATGCGCACCTGTGTCCCCCTCCCCCTTCCCGGGTTTGGGGGGGGGTGTGTGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                           <  1:4208782/149‑1 (MQ=16)
                                                                                                          aCCCTGTCATTAACATGAGCACCGTTTTCTCCCCCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                           <  2:4393320/150‑1 (MQ=255)
                                                                                                          aCCCTGACATTAACATGAGCACCGTTTTCCCCCTCCCCCTCCCGGGGTGAGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                           <  1:5422279/150‑1 (MQ=25)
                                                                                                            ccgttCATTAAATTGAGCTCGGTTTTCACCCTCACCCTCCGGGGGTGTGGGGGGGGGTGTGGTCAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGa                                         <  1:1807523/146‑1 (MQ=21)
                                                                                                            ccTGTCATACACATGCGCACCTTTTCCCCCCCCCCCCTCCCGGGGTTTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGa                                         <  1:4824084/150‑1 (MQ=21)
                                                                                                             cTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAt                                        <  2:4481323/150‑1 (MQ=255)
                                                                                                             cTGACATTAACATGCGCCCCGTTTTCCCCCTCCCCCTACCCGGGAGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAt                                        <  2:2456558/150‑1 (MQ=21)
                                                                                                                tCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgtt                                     <  1:5511497/150‑1 (MQ=255)
                                                                                                                  aTTAACATGAGCACCGTTTTCTCCCTCGCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgttgt                                   <  2:3715526/150‑1 (MQ=255)
                                                                                                                    tAACAGGAGCACCTTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtca                                 <  1:5855940/150‑1 (MQ=255)
                                                                                                                        aTGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAg                             >  1:1679203/1‑150 (MQ=255)
                                                                                                                        aTGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAg                             >  2:1399940/1‑150 (MQ=255)
                                                                                                                        aTGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCACCCCAAACGGGGCGAGGGGACAGACCGAGTCCTTTTTTGATGTTGtcaacac                             >  2:1080380/1‑149 (MQ=255)
                                                                                                                          gAGCACCGTTTTCTCCCTCTCCCTTCCAGGGTGCGGGTCGGGGTTTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTc                           <  2:1010049/150‑1 (MQ=255)
                                                                                                                             cACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGGGTCCTTTTTTGATGTTGTCATCAGTCTgg                        >  1:4211821/1‑150 (MQ=255)
                                                                                                                              aCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGa                       >  2:5971773/1‑150 (MQ=255)
                                                                                                                              aCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGa                       >  2:917014/1‑150 (MQ=255)
                                                                                                                               ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                      >  2:4210633/1‑150 (MQ=255)
                                                                                                                               ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                      >  2:966746/1‑150 (MQ=255)
                                                                                                                                      ctccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGGTGTCATCAGTCTGGAAGCCGccc               >  2:765209/1‑148 (MQ=255)
                                                                                                                                       tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACg              >  2:4029876/1‑150 (MQ=255)
                                                                                                                                         cctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt            >  2:215019/1‑150 (MQ=255)
                                                                                                                                         cctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTACCCGCACCCAAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGTTGTTGTCATCAGTCTGGAAGCCGCACGtt            >  2:2561217/1‑150 (MQ=255)
                                                                                                                                         cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt            >  1:2314554/1‑150 (MQ=255)
                                                                                                                                           cccccctTCCGGGGTTTGGGGGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg          <  2:4241664/147‑1 (MQ=37)
                                                                                                                                            ctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCACCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGc         >  2:3538107/1‑150 (MQ=255)
                                                                                                                                                 ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGCTGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGCCTGACCGAGTCCTGTTTTGATGGTGTCATCAGTCTGGACGCCGCACGGTGGCGTTAt    >  1:3860641/1‑150 (MQ=25)
                                                                                                                                                 ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGTTGTTGTGATCAGTCTGGATGCCGCACGTTGGGTTTAt    >  1:5709643/1‑150 (MQ=255)
                                                                                                                                                 ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt    >  1:5315564/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTT  >  NC_000913/3239463‑3239759

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