Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,778 C→T 6.7% intergenic (+324/+13) 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,7780CT 6.7% 280.6 / 2.3 101intergenic (+324/+13)codA/cynRcytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (64/30);  new base T (2/5);  total (66/35)
Fisher's exact test for biased strand distribution p-value = 4.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

CAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA  >  NC_000913/357630‑357926
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                                                   >  2:485839/1‑98 (MQ=37)
                                                      gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCtg                                                                                                                                                   <  1:485839/98‑1 (MQ=37)
                                                      gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCtg                                                                                                                                                   <  1:996523/98‑1 (MQ=37)
                                                       ccGGATGCGGCGTGACCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGACGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCCACAACCGCGCATTTCTCAAACg                                                                                              >  2:1509604/1‑150 (MQ=11)
                                                       ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACg                                                                                              <  1:514936/150‑1 (MQ=255)
                                                        cGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGc                                                                                             >  1:666412/1‑150 (MQ=255)
                                                        cGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGc                                                                                             >  2:845461/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:1172402/1‑150 (MQ=255)
                                                                                                             ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                        >  2:1022684/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)
                                                                                                                                                   aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga  <  1:825033/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA  >  NC_000913/357630‑357926

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: