Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 508,795 A→G 7.1% intergenic (+236/+80) ybaQ → / ← copA putative DNA‑binding transcriptional regulator/copper transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913508,7950AG 7.1% 291.4 / 9.3 98intergenic (+236/+80)ybaQ/copAputative DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base A (54/37);  new base G (2/5);  total (56/42)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATT  >  NC_000913/508647‑508944
                                                                                                                                                    |                                                                                                                                                     
cccGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACATGTGCCTGAACCGTAGGCCGGATaa                                                                                                                                                      >  1:200680/1‑150 (MQ=255)
  cGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTTTCAGGCCTACAAACCTGTGCCTGCACCGTAGGCCGGATAAgg                                                                                                                                                    >  2:508714/1‑150 (MQ=255)
   gTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGc                                                                                                                                                   >  2:1420991/1‑150 (MQ=255)
   gTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGc                                                                                                                                                   >  2:238200/1‑150 (MQ=255)
    tCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCAGATGCGACGCAAACGAGTCTTATCAGGCCAAAAAACAACTGCCAGAACCCTAGGCCGGAAAAGCCg                                                                                                                                                  >  1:128446/1‑150 (MQ=18)
     cTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCCCCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGACCCGTAGGCCGGATAAGGCGt                                                                                                                                                 >  2:585047/1‑150 (MQ=255)
       tATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGttt                                                                                                                                               >  1:1011778/1‑150 (MQ=255)
        aTCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTa                                                                                                                                              >  2:1173152/1‑150 (MQ=255)
           aGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                           >  2:1553864/1‑150 (MQ=255)
             gCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                         <  2:721887/150‑1 (MQ=255)
                tACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                                                                                      >  2:1160525/1‑150 (MQ=255)
                 aCAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATGCGGCGTAAAcgccttatc                                                                                                                                     <  1:847042/150‑6 (MQ=17)
                  cAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                    <  1:426949/150‑1 (MQ=255)
                   aaaCTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                   <  2:409854/150‑1 (MQ=17)
                    aaCTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                  <  2:898968/150‑1 (MQ=17)
                    aaCTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                  >  1:1669940/1‑150 (MQ=255)
                    aaCTTGTGCCTGAACCGTACGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                  <  1:1394778/150‑1 (MQ=17)
                     ccTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                 <  1:51546/149‑1 (MQ=14)
                           gCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCAt                                                                                                                           >  1:1156227/1‑150 (MQ=255)
                            cccgacccGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTCTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATg                                                                                                                          <  2:1635604/144‑1 (MQ=255)
                            ccTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATg                                                                                                                          >  1:1478755/1‑150 (MQ=255)
                            ccTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTATTATCAGGCCTACAAACCTGTGCCTGAACCCTAGGCCGGATAAGGCGTTTACGCCGCAAACGGCAGTCATg                                                                                                                          >  1:270819/1‑150 (MQ=255)
                              tGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACCAACCTGTGCCTGAACCGTAGGCCGGAGAAGGCGTTTACGCCGCATCCGGCAGTCATGCg                                                                                                                        >  2:244345/1‑150 (MQ=255)
                              cgAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCCTACGCGTCTTATCAGGCCGACAAACCTGCGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCg                                                                                                                        <  1:1642289/149‑1 (MQ=255)
                                 aCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCAATGCCTGATGCGACGCTTACGCGACGTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGAAAAGCCGTTTACGCCGCATCCGGCAGTCATGCGTcc                                                                                                                     >  2:249452/1‑149 (MQ=14)
                                      aGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGcc                                                                                                                >  2:1681462/1‑150 (MQ=255)
                                       ggCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCa                                                                                                               <  1:1350021/150‑1 (MQ=255)
                                        gCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCaa                                                                                                              >  1:1030293/1‑150 (MQ=255)
                                        gCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCaa                                                                                                              >  1:1379030/1‑150 (MQ=255)
                                                 ggCGTTTACGCCTCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCAc                                                                                                     <  1:271439/150‑1 (MQ=255)
                                                 ggCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAACCCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAAAGCGCCAc                                                                                                     >  1:193535/1‑150 (MQ=255)
                                                 ggCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCAc                                                                                                     <  1:1681462/150‑1 (MQ=255)
                                                 ggCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCAc                                                                                                     >  1:1409537/1‑150 (MQ=255)
                                                  gCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCAcc                                                                                                    <  1:979685/150‑1 (MQ=255)
                                                      ttACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTaa                                                                                                >  2:523952/1‑150 (MQ=255)
                                                       tACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTaaa                                                                                               >  2:403727/1‑150 (MQ=255)
                                                       tACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTaaa                                                                                               >  2:1574125/1‑150 (MQ=255)
                                                        aCGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAag                                                                                              <  1:1138625/150‑1 (MQ=255)
                                                         cgccgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAACCCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCAACCCTAAagc                                                                                             >  1:1134866/1‑150 (MQ=255)
                                                         cgccgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAagc                                                                                             >  1:1338886/1‑150 (MQ=255)
                                                           ccgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGACGGTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCCCCCTAAagcag                                                                                           >  1:163765/1‑150 (MQ=255)
                                                             gcatccggcatccgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAgcg                                                                                         <  1:585047/150‑1 (MQ=255)
                                                                 ccggcatccgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATc                                                                                     <  2:1379030/150‑1 (MQ=255)
                                                                   ggcatccgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCg                                                                                   <  1:164254/150‑1 (MQ=255)
                                                                     catccgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCa                                                                                 >  1:310090/1‑150 (MQ=255)
                                                                         cgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAACCCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAatga                                                                             >  1:451756/1‑150 (MQ=255)
                                                                         cgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAatga                                                                             >  2:1124357/1‑150 (MQ=255)
                                                                            gCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGt                                                                          <  2:1514793/150‑1 (MQ=255)
                                                                             cACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTa                                                                         >  1:468962/1‑150 (MQ=255)
                                                                                    gCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTc                                                                  >  2:140213/1‑150 (MQ=255)
                                                                                    gCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTc                                                                  <  1:728135/150‑1 (MQ=255)
                                                                                    gCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTc                                                                  >  2:362282/1‑150 (MQ=255)
                                                                                     ccTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTcc                                                                 >  1:1576361/1‑150 (MQ=255)
                                                                                        gATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTc                                                              <  1:1420991/150‑1 (MQ=255)
                                                                                        gATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTc                                                              <  1:238200/150‑1 (MQ=255)
                                                                                            cGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGtt                                                          <  1:966616/150‑1 (MQ=255)
                                                                                            cGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGtt                                                          >  1:905215/1‑150 (MQ=255)
                                                                                            cGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGtt                                                          <  2:677319/150‑1 (MQ=255)
                                                                                            cGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGtt                                                          <  1:1629421/150‑1 (MQ=255)
                                                                                             gACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGttt                                                         >  1:529548/1‑150 (MQ=255)
                                                                                                 cTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAAc                                                     <  1:1173152/150‑1 (MQ=255)
                                                                                                          cTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAAc                                            <  1:617840/150‑1 (MQ=255)
                                                                                                           gtATCAGGCCTACACACCTGTGCCTGAACCGTAGGCCGGATAACGCGTTTACGCCGCATCCGGCAGTCATGCGTCGACGCCACATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAAcc                                           <  2:1018795/149‑1 (MQ=255)
                                                                                                            tATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATACGTCGATGCAAAATGCGACACCCTAAAGCAGCGCATCCGCAAATAACTACTTATTCCTTCGGTTTAAACCGCAGCACCCg                                          >  2:276744/1‑150 (MQ=25)
                                                                                                                aGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTg                                      <  1:1716828/150‑1 (MQ=255)
                                                                                                                aGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGAATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAACGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGGTg                                      >  2:515473/1‑150 (MQ=255)
                                                                                                                 ggCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTgg                                     >  2:1488341/1‑150 (MQ=255)
                                                                                                                  gCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGc                                    <  2:1669940/150‑1 (MQ=255)
                                                                                                                   ccTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCg                                   >  2:1478729/1‑150 (MQ=255)
                                                                                                                    cTACAAACCTGTGCCTGAACCGGAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGt                                  >  2:644747/1‑150 (MQ=255)
                                                                                                                     tACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGGTGGCGtt                                 <  2:625217/150‑1 (MQ=255)
                                                                                                                     caCAAACCCGCCCCAGAACCGTACGCCGGACAAGGCGTCTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGCTTAACCCGCAGCAACCGGTTGGCGtt                                 <  2:1122906/149‑1 (MQ=38)
                                                                                                                      aCAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTa                                >  2:1238817/1‑150 (MQ=255)
                                                                                                                      aCAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATACCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTa                                >  2:995948/1‑150 (MQ=255)
                                                                                                                       cAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTAc                               <  1:931873/150‑1 (MQ=255)
                                                                                                                         aaCCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTc                             >  1:360508/1‑150 (MQ=255)
                                                                                                                         aaCCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTc                             <  1:1709648/150‑1 (MQ=255)
                                                                                                                         aaCCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGAAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTAATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTc                             >  2:640573/1‑150 (MQ=255)
                                                                                                                              gtgCCTGAACCGTATGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGACGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTAc                        <  1:707389/150‑1 (MQ=255)
                                                                                                                              gtgCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTAc                        <  1:1160525/150‑1 (MQ=255)
                                                                                                                               tgCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACg                       >  1:527892/1‑150 (MQ=255)
                                                                                                                               ggCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCAGCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACg                       <  2:873365/149‑1 (MQ=255)
                                                                                                                                gCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACgg                      >  2:752475/1‑150 (MQ=255)
                                                                                                                                gCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACgg                      >  1:1566716/1‑150 (MQ=255)
                                                                                                                                 ccTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGt                     <  1:221748/150‑1 (MQ=255)
                                                                                                                                 ccTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGt                     >  1:32297/1‑150 (MQ=255)
                                                                                                                                 ccTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGt                     >  1:688264/1‑150 (MQ=255)
                                                                                                                                 ccTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGt                     >  2:1392828/1‑150 (MQ=255)
                                                                                                                                 ccTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGAGGTCCTTATTCCTTCGGTTTAAACCGCAGCACCCGGTTGGCGTTACTCACTACGGt                     >  2:130533/1‑150 (MQ=255)
                                                                                                                                  cTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTa                    <  2:128446/150‑1 (MQ=255)
                                                                                                                                   tGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTaa                   >  1:307611/1‑150 (MQ=255)
                                                                                                                                   tGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTaa                   >  1:706188/1‑150 (MQ=255)
                                                                                                                                         gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATcgacg             <  2:1086792/150‑1 (MQ=255)
                                                                                                                                         gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATcgacg             <  2:47949/150‑1 (MQ=255)
                                                                                                                                         gTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATcgacg             <  2:163765/150‑1 (MQ=255)
                                                                                                                                          actgCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATcgacga            >  2:991315/4‑150 (MQ=255)
                                                                                                                                           ctgCCGGATGCGGCGTAAACGCCTTATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACgag           >  2:1593435/3‑150 (MQ=255)
                                                                                                                                           ctgCCGGATGCGGCGTAAACGACTTATCCGGCAGTCATGCGTCGATGCCAAATGCGCCAACCTAAAGCAGCGCATCAGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGCTGGGCGATACACACTACGGTAATcgccgac           >  1:1290174/3‑149 (MQ=255)
                                                                                                                                                   tAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCAt   <  1:523952/150‑1 (MQ=255)
                                                                                                                                                    aaGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCAtt  <  1:276744/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATT  >  NC_000913/508647‑508944

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: