Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,775 A→T 7.8% intergenic (+321/+16) codA → / ← cynR cytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,7750AT 7.8% 254.0 / 3.6 101intergenic (+321/+16)codA/cynRcytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (62/29);  new base T (3/6);  total (66/35)
Fisher's exact test for biased strand distribution p-value = 6.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

TCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG  >  NC_000913/357628‑357921
                                                                                                                                                   |                                                                                                                                                  
tccAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTCTATTTTATGCCCCAAACAGCCTGAACGCCATCAACTGCATACGTAGCAATATGAACTGGTAGGCCTGCTAAGAGTAGCGAATCAGTCAATTCAAGCCGCaca                                                                                                                                                  >  2:288289/1‑148 (MQ=12)
  cAACGCCATGTGCAAAAATGCCTTCGCCGCGGGTGTCGTCCTGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGAcc                                                                                                                                                <  2:881765/150‑1 (MQ=18)
    aCGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCtg                                                                                                                                              <  1:666695/150‑1 (MQ=255)
      gCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCtgtg                                                                                                                                            <  2:915368/150‑1 (MQ=255)
      gCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCCGAAGCGGCGTGACCGCCTTATCCGGCCGCCGTAGAGCACGGAACTCGTAGGCCAGAGAAGCGTAGCGCGTCAGTCAATGCCAGCCGCACACCggtg                                                                                                                                            >  1:648558/1‑150 (MQ=37)
      gCCATGTCCAAAAATGCCTCCCCCGCGGCGGGCTGCCAGCTGTAGTATATGCCTGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCtgtg                                                                                                                                            <  2:191975/150‑1 (MQ=17)
      gCCATGGGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCtgtg                                                                                                                                            <  1:599327/150‑1 (MQ=255)
            tgCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCg                                                                                                                                      <  1:482024/150‑1 (MQ=255)
            tgCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGACCTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCCATTCCAGCCGCAGACCTGTGTCAGCg                                                                                                                                      >  1:1932619/1‑150 (MQ=255)
            ggCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCg                                                                                                                                      <  2:1780260/149‑1 (MQ=255)
             gCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCTCAGACCTGTGTCAGCgg                                                                                                                                     >  1:953886/1‑150 (MQ=255)
               aaaaaTGCCTTCGCAGCGGCTGTCTGCCAGCTATATTTAATCCCCTAACCAACATCAACTCAATAACCGAAATACGTAGAGCTCTGAACACATAGGCATAATAAGCGTATCGCAACAGGAAATTACAGCAGCAGACCTGTGTCAGCGGCt                                                                                                                                   >  2:1020367/1‑150 (MQ=255)
                     gCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTg                                                                                                                             >  1:1792655/1‑150 (MQ=255)
                     gCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTg                                                                                                                             >  1:1656597/1‑150 (MQ=255)
                       cTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGAGCGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACATGTTTCAGCGGCGACCGTGAt                                                                                                                           >  1:1542595/1‑150 (MQ=17)
                        ctcccAGCGGCCGGCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAtt                                                                                                                          <  1:1128808/146‑1 (MQ=255)
                         ccGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCTTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTc                                                                                                                         <  1:1001959/149‑1 (MQ=255)
                            aaGCGGCGGTATGCCATCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCAtt                                                                                                                      <  2:1268085/149‑1 (MQ=255)
                             aGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCAttt                                                                                                                     >  2:649385/1‑150 (MQ=255)
                                ggCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCg                                                                                                                  >  2:897323/1‑150 (MQ=255)
                                 gCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAAGTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGc                                                                                                                 >  1:105374/1‑150 (MQ=255)
                                          cAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcg                                                                                                        >  1:662962/1‑150 (MQ=255)
                                           aGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAAAcgc                                                                                                       >  2:568479/1‑150 (MQ=255)
                                            gCTGTAGTTTATGCCGGATGCGGCGTGACCGCCTTATCCGGCCGACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCCGCCGCAGCCCTGTGTCAGCGGCTCCCGTGAGTCAGTTCCGCCAACCCCcgcg                                                                                                      >  1:782035/1‑150 (MQ=14)
                                            gCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGGTTCATTTCCGCCAACAACcgcg                                                                                                      >  1:774841/1‑150 (MQ=255)
                                             ttttAGTTTATGCCGGATGCGGCGTGAACGCATTATCAGGCGTATGTAGATCAGTGTACTCGTCTTCCTTATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCgcgc                                                                                                     <  1:664602/147‑1 (MQ=11)
                                                 aGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCAttt                                                                                                 >  1:527475/1‑150 (MQ=38)
                                                   tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTTCCGTGATTCATTTCCGCCAACAACCGCGCATGTAt                                                                                               >  1:1132168/1‑150 (MQ=255)
                                                   tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAt                                                                                               >  1:1572534/1‑150 (MQ=39)
                                                      aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATccc                                                                                            >  2:174558/3‑149 (MQ=14)
                                                       agCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGATCTCGTAGGCCTGATTAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCGACCGTGATTCATTTCCGCCAACACCCGCGCATTTATCCac                                                                                           >  1:362172/2‑149 (MQ=255)
                                                        gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCtg                                                                                                                                              <  1:996523/98‑1 (MQ=37)
                                                        gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCtg                                                                                                                                              >  2:485839/1‑98 (MQ=37)
                                                        gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCtg                                                                                                                                              <  1:485839/98‑1 (MQ=37)
                                                         ccGGATGCGGCGTGACCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGACGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCCACAACCGCGCATTTCTCAAACg                                                                                         >  2:1509604/1‑150 (MQ=11)
                                                         ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACg                                                                                         <  1:514936/150‑1 (MQ=255)
                                                          cGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGc                                                                                        >  2:845461/1‑150 (MQ=255)
                                                          cGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGc                                                                                        >  1:666412/1‑150 (MQ=255)
                                                              tGCGGCGTGAACGCCTTATCCGGCCTACTTAGAGCACTGAACTCGTCGGCCTGATAAGCGTAGCGCATAAGGCAATTCCAGCCGCAGACGTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGAGCGTTTATCCACCGCCAtg                                                                                    >  2:1167114/1‑150 (MQ=14)
                                                              tGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCATGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTCCCGTGAGTCATTTCCGCCGACAACCGCGCATTTATCCACGGCCAtg                                                                                    >  2:969483/1‑150 (MQ=14)
                                                                cggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTCGGCCTGATAAGCGTAGCGCATCAGCCAATTCCAGCCGCGGACCTGAGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCCTATATCCAACGCCAtgtg                                                                                  >  2:92273/1‑150 (MQ=14)
                                                                 ggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTGCCGTGATTCATTTCCGc                                                                                                                 >  1:1009317/1‑118 (MQ=12)
                                                                 ggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTGCCGTGATTCATTTCCGc                                                                                                                 <  2:1009316/118‑1 (MQ=255)
                                                                 ggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGc                                                                                 >  1:496465/1‑150 (MQ=255)
                                                                      gAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCCACGCCATGTGCaaaaa                                                                            >  1:816056/1‑150 (MQ=255)
                                                                      gAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaaaat                                                                            >  1:1579332/1‑149 (MQ=255)
                                                                      gAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaaaaa                                                                            >  2:681065/1‑150 (MQ=255)
                                                                        aCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATg                                                                          <  2:644355/150‑1 (MQ=255)
                                                                        aCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATg                                                                          >  2:333082/1‑150 (MQ=255)
                                                                         cGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCCGCGGCTCCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATgg                                                                         >  1:175620/1‑149 (MQ=255)
                                                                         cGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGc                                                                         >  1:88360/1‑150 (MQ=255)
                                                                         cGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGc                                                                         >  1:962419/1‑150 (MQ=255)
                                                                          gCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGcc                                                                        >  2:1398652/1‑150 (MQ=255)
                                                                          gCCTTATCCGGCCTACGTAGAGCACTGAAATCGTAGGCCTGATAAGCGTAGAGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGAAAAAATGcc                                                                        >  1:1569038/1‑150 (MQ=255)
                                                                           ccTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTg                                                                                                                             >  2:85760/1‑96 (MQ=37)
                                                                           ccTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTg                                                                                                                             <  1:85760/96‑1 (MQ=37)
                                                                           ccTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCt                                                                       >  1:1863307/1‑150 (MQ=255)
                                                                                 ccGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCgcgg                                                                 >  1:1059896/1‑150 (MQ=255)
                                                                                        aCGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGt                                                          >  1:1925836/1‑150 (MQ=255)
                                                                                         cGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTc                                                         >  2:1124306/1‑150 (MQ=255)
                                                                                           tAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTTCAAAAATGCCTTCGCGGCGGCTGTCTg                                                       <  2:1656596/150‑1 (MQ=255)
                                                                                           tAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTg                                                       <  2:1792654/150‑1 (MQ=255)
                                                                                            agagCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAtt                                                                                                                          <  2:1782962/82‑1 (MQ=38)
                                                                                            agagCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAtt                                                                                                                          >  1:1782963/1‑82 (MQ=38)
                                                                                             gagCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGcc                                                     >  1:1253875/1‑150 (MQ=255)
                                                                                              agCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCa                                                    <  2:496465/150‑1 (MQ=255)
                                                                                                cACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGGAAAAATGCCTTCGCGGCGGCTGTCTGCCAGc                                                  >  2:401528/1‑150 (MQ=255)
                                                                                                  cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaaaa                                                                             >  1:1438671/1‑121 (MQ=17)
                                                                                                  cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaaaa                                                                             <  2:1438670/121‑1 (MQ=17)
                                                                                                  cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                >  2:1198437/1‑150 (MQ=255)
                                                                                                    gAACTCGTAGGCCTGATAAGAGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCAGTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTaat                                              >  1:399457/1‑148 (MQ=255)
                                                                                                      aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGc                                                                                                                                       <  2:1364181/59‑1 (MQ=255)
                                                                                                      aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGc                                                                                                                                       >  1:1364182/1‑59 (MQ=255)
                                                                                                       cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                           >  1:222810/1‑150 (MQ=255)
                                                                                                        tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                          <  2:1883218/150‑1 (MQ=255)
                                                                                                         cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCCGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAGGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCGAGTTTTcc                                         >  1:1817316/1‑150 (MQ=255)
                                                                                                               ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACg                                                                                         >  1:398365/1‑96 (MQ=16)
                                                                                                               ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACg                                                                                         <  2:398365/96‑1 (MQ=255)
                                                                                                               ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                   >  2:1022684/1‑150 (MQ=255)
                                                                                                               ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                   >  2:1172402/1‑150 (MQ=255)
                                                                                                               ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAGTTTTCCGCCGca                                   >  1:1419113/1‑150 (MQ=255)
                                                                                                                      aGCGTAGCGCATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAACCGCGCATTTATCCAACGCCATGTGc                                                                                 <  1:1438656/97‑1 (MQ=255)
                                                                                                                      aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGc                                                                                                                                    <  1:223620/46‑1 (MQ=255)
                                                                                                                      aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                            >  2:935059/1‑150 (MQ=255)
                                                                                                                           aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                       >  1:1481658/1‑150 (MQ=255)
                                                                                                                            gcgcATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCATCAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                      <  1:1568588/150‑1 (MQ=255)
                                                                                                                             cgcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAATATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGtt                     <  2:368780/150‑1 (MQ=255)
                                                                                                                              gcATCAGGCAATTCCAGCCGCAGACCTGTGGCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                    <  1:632874/150‑1 (MQ=255)
                                                                                                                                   aGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTAATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc               <  2:924599/150‑1 (MQ=255)
                                                                                                                                    ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTGTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCGCGCCa              >  1:1881569/1‑150 (MQ=255)
                                                                                                                                    ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCa              >  1:1938104/1‑150 (MQ=255)
                                                                                                                                      cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAAACGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCaga            >  1:669912/1‑149 (MQ=255)
                                                                                                                                        aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGa                                                                                                                            <  1:381724/36‑1 (MQ=255)
                                                                                                                                        aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCAATGATCCAACGCCATGTGCAAAAATGACTTCGCGGCGGCTGTCTGACAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtag          >  2:342536/1‑150 (MQ=255)
                                                                                                                                           ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg       >  2:1749898/1‑150 (MQ=255)
                                                                                                                                            cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc      >  1:800366/1‑150 (MQ=255)
                                                                                                                                               ccGCAGACCTGTGTCATCGTCTACAGAGATTCATTTCCGCCAACAACCGCCCTTTTATCCAACTCCATGTTCAATAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg   <  2:105374/150‑1 (MQ=255)
                                                                                                                                               ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg   >  2:739863/1‑150 (MQ=255)
                                                                                                                                               cagCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCAAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGTCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg   <  2:378351/148‑1 (MQ=255)
                                                                                                                                                cGCAGACCTGAGTCAGCCGCGACCGGGTAACCTAGCCGACAAGACCCGCCCACTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg  <  1:1120132/150‑1 (MQ=25)
                                                                                                                                                   |                                                                                                                                                  
TCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG  >  NC_000913/357628‑357921

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