Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 508,803 T→A 8.3% intergenic (+244/+72) 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,8030TA 8.3% 264.4 / 11.5 97intergenic (+244/+72)ybaQ/copAputative DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base T (50/39);  new base A (3/5);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

Only 100 of 101 total aligned reads displayed.

TATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCC  >  NC_000913/508654‑508952
                                                                                                                                                     |                                                                                                                                                     
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:362282/1‑150 (MQ=255)
                                                                             gCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTc                                                                          >  2:140213/1‑150 (MQ=255)
                                                                             gCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTc                                                                          <  1:728135/150‑1 (MQ=255)
                                                                              ccTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTcc                                                                         >  1:1576361/1‑150 (MQ=255)
                                                                                 gATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTc                                                                      <  1:238200/150‑1 (MQ=255)
                                                                                 gATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTc                                                                      <  1:1420991/150‑1 (MQ=255)
                                                                                     cGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGtt                                                                  <  2:677319/150‑1 (MQ=255)
                                                                                     cGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGtt                                                                  <  1:1629421/150‑1 (MQ=255)
                                                                                     cGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGtt                                                                  <  1:966616/150‑1 (MQ=255)
                                                                                     cGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGtt                                                                  >  1:905215/1‑150 (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)
                                                                                                                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                             >  2:1392828/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)
                                                                                                                          ccTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGAGGTCCTTATTCCTTCGGTTTAAACCGCAGCACCCGGTTGGCGTTACTCACTACGGt                             >  2:130533/1‑150 (MQ=255)
                                                                                                                           cTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTa                            <  2:128446/150‑1 (MQ=255)
                                                                                                                            tGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTaa                           >  1:706188/1‑150 (MQ=255)
                                                                                                                            tGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTaa                           >  1:307611/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)
                                                                                                                                              cggCGTAAACGCCTTATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATCATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTg         >  1:628087/2‑150 (MQ=255)
                                                                                                                                              aGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTg         <  1:1187173/150‑1 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGc        >  2:85709/1‑150 (MQ=255)
                                                                                                                                                 cGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCg      >  2:1441194/1‑150 (MQ=255)
                                                                                                                                                    ttACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTc   >  1:1308627/1‑150 (MQ=255)
                                                                                                                                                     tACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGTCGCTcc  >  2:176317/1‑150 (MQ=255)
                                                                                                                                                     tACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTcc  <  1:539739/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCC  >  NC_000913/508654‑508952

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: