Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 131,353 A→G 5.2% intergenic (‑93/‑262) yacH ← / → acnB DUF3300 domain‑containing protein YacH/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913131,3530AG 5.2% 255.1 / 3.4 77intergenic (‑93/‑262)yacH/acnBDUF3300 domain‑containing protein YacH/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base G (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

AACGTAACCGCGGCAGCGTCAGACACAGCCAGTTGTGTCGATTGCGGTTCCACAGGCGCTTCCACTGTGCGGCTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGCCCGTTGCGGTACTGGGCATTTACCCTACAAACTGCTGTCTCACAGGAGCGTGAAGAGAATCGCCTGCCGCACTATGACAATG  >  NC_000913/131097‑131596
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
aaCGTAACCGGTGAAGCGTACGAAACAGTGTGTTGTGGCGATTGCGGTTCCACAGGCGCTTCCACTGGGCGGCTTTTTATATATAGCCCGGTAGAGGCTGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTc                                                                                                                                                                                                                                        <  2:162499/272‑1 (MQ=255)
                cGTCAAACACAGCCAGTTGTGTCGATTGCGGTTCCACAGGCGCTTCCACTGTGCGGCTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTgaga                                                                                                                                                                                                                                 <  1:108098/263‑1 (MQ=255)
                               gTTGTGTCGATTGCGGTTCCACAGGCGCTTCCACTGTGCGGCTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCCAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGACCTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAAtt                                                                                                                                                                                                         >  1:746821/1‑272 (MQ=255)
                                       gATTGCGGTTCCACAGGCGCTTCCACTGTGCGGCTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTTAATCCGAGTATAAAGAGGCGGGAGTTTAAATTTTGACTAATCTTGGGATTCGTTGTGAAAGGTGATTATCACCATTCGAATTAACGAAGt                                                                                                                                                                                                 >  1:29872/1‑272 (MQ=255)
                                                              cccggcgCGGCTTGTTATATCTGGCCAGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACAGGGGGTTTAAACGGCAAAGTCATTTTCATAAATCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATctct                                                                                                                                                                        <  2:746820/266‑1 (MQ=255)
                                                              cACTTTGCGGCTTTTTATATCTAGCCCGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATctct                                                                                                                                                                        <  1:968254/272‑1 (MQ=255)
                                                                    gctttttttttttTATCCCCCGGTGGGGGGGGAAAAAAGCCCCGGAACTGAAAATTAGCGCAAGAACATGGGGTTTAACCGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAg                                                                                                                                                                      <  2:29872/264‑1 (MQ=255)
                                                                      ggCTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAACTTTTGACTAATCTTGGGATTCGTTTAGAAAGGTTATTATCACCATGCGAATTTAAGAAGTTTTTACGGGAGGAAAACaaccctagaacct                                                                                                                                                                  >  2:1012784/1‑260 (MQ=255)
                                                                      ggCTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTg                                                                                                                                                                                                                                              >  1:298828/1‑196 (MQ=255)
                                                                      ggCTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGtt                                                                                                                                                                                                >  1:818542/1‑242 (MQ=255)
                                                                      ggCTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAggg                                                                                                                                                                                   >  1:276392/1‑253 (MQ=255)
                                                                      ggCTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCCAGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAACGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGGTGAGAAAGGTGATTAAACACAGGCGAATAAAAGAAGATTTTACGGAGGGGAAACaaccatagcacaaa                                                                                                                                                                 >  2:53761/1‑260 (MQ=255)
                                                                       gCTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTgg                                                                                                                                                                                                                                             <  2:298827/196‑1 (MQ=255)
                                                                       gCTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGttt                                                                                                                                                                                               <  2:818541/242‑1 (MQ=255)
                                                                       gCTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGa                                                                                                                                                                                  <  2:276391/253‑1 (MQ=255)
                                                                        cTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCACGGCACTGCAAATTAGCGCCAGCACATGGGGTTTACACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCAACCCGCAGGTGCTGTGAATCCGAGTATACAGAGGCGGTAGTTTAAATTTTGACTAATATTGTGTTTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTGTACGGAGGGAAACAAACTATAAACCAGCCt                                                                                                                                                              >  1:557315/1‑272 (MQ=255)
                                                                                         ggTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGAGTAATCTTGGGAGTCGTTGAGAAAGGTGAATATCACCATGGGCATTAACGAAGTTTTTACTGATGGAAACAACTTCAAGACGAACCTTAACGaataaagcactt                                                                                                                                             >  2:315918/1‑261 (MQ=255)
                                                                                                      caAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCAGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTgag                                                                                                                                                                                                                                  >  1:551205/1‑176 (MQ=255)
                                                                                                       aaaGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTgaga                                                                                                                                                                                                                                 <  2:551204/176‑1 (MQ=255)
                                                                                                             ccGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTgag                                                                                                                                                                                                                                  >  2:444645/1‑170 (MQ=255)
                                                                                                              cGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTgaga                                                                                                                                                                                                                                 <  1:444646/170‑1 (MQ=255)
                                                                                                                  aCTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTgag                                                                                                                                                                                                                                  >  2:125938/1‑165 (MQ=255)
                                                                                                                   cTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTgaga                                                                                                                                                                                                                                 <  1:125938/165‑1 (MQ=255)
                                                                                                                       aaaTTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCtt                                                                                                                                                                                                                                               >  1:381210/1‑147 (MQ=255)
                                                                                                                       aaaTTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGCAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTa                                                                                                                >  1:792426/1‑272 (MQ=255)
                                                                                                                        aaTTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTg                                                                                                                                                                                                                                              <  2:381209/147‑1 (MQ=255)
                                                                                                                         aTTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAggg                                                                                                                                                                                   >  1:699749/1‑203 (MQ=255)
                                                                                                                         aTTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAAACCATCCTTAACGATTCAGCCACTTTTTTATGTTTCTTTTTTTTAAACAAATTaa                                                                                                               >  1:815884/1‑271 (MQ=255)
                                                                                                                          ttAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGa                                                                                                                                                                                  <  2:699748/203‑1 (MQ=255)
                                                                                                                            aGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGAtt                                                                                                                                                                                                                                         >  1:283649/1‑148 (MQ=255)
                                                                                                                             gcgcCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTc                                                                                                                                                                                                                                        <  2:283648/148‑1 (MQ=255)
                                                                                                                                 cAGCACATGGGGTTTAAACGGTAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTa                                                                                                                                                                                            >  1:385448/1‑188 (MQ=255)
                                                                                                                                 cAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTACCGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGt                                                                                                       >  1:280186/1‑271 (MQ=255)
                                                                                                                                  aGCACATGGGGTTTAAACGGTAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTAc                                                                                                                                                                                           <  2:385447/188‑1 (MQ=255)
                                                                                                                                       attgggTTTAAACGGCAAAGTCATTTTCATAATTCTGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATc                                                                                                <  1:1012785/269‑1 (MQ=255)
                                                                                                                                       aTGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATc                                                                                                <  1:315919/272‑1 (MQ=255)
                                                                                                                                                 aaagaaaagtgaaTTTTAATAATTCGAATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCt                                                                                      <  2:557314/260‑1 (MQ=255)
                                                                                                                                                 aaCGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGGCTAATCTTGGGATTCGtt                                                                                                                                                                                                                                     >  1:927421/1‑131 (MQ=255)
                                                                                                                                                 aaCGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAggg                                                                                                                                                                                   >  1:580290/1‑179 (MQ=255)
                                                                                                                                                  aCGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGGCTAATCTTGGGATTCGTTg                                                                                                                                                                                                                                    <  2:927420/131‑1 (MQ=255)
                                                                                                                                                  aCGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGa                                                                                                                                                                                  <  2:580289/179‑1 (MQ=255)
                                                                                                                                                  aCGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACACGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGAtt                                                                                                                                                                                                                                         >  1:835554/1‑126 (MQ=255)
                                                                                                                                                   cGGCAAATTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCTGCGAACTGCCACCCGCAGGTGCTGTGAAGCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCt                                                                                      <  1:717420/270‑1 (MQ=255)
                                                                                                                                                   cGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTc                                                                                                                                                                                                                                        <  2:835553/126‑1 (MQ=255)
                                                                                                                                                        aaGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGTATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGa                                                                                                                                                                                  >  1:583455/1‑173 (MQ=255)
                                                                                                                                                        aaGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACATCGTCAAAATCCTGCTCTTCTGCCCGtt                                                                               >  1:575209/1‑272 (MQ=255)
                                                                                                                                                         aGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGaa                                                                                                                                                                                 <  2:583454/173‑1 (MQ=255)
                                                                                                                                                         aGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTg                                                                                                                           >  1:666713/1‑227 (MQ=255)
                                                                                                                                                          gTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGt                                                                                                                          <  2:666712/227‑1 (MQ=255)
                                                                                                                                                          gTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCa                                                                                                     <  1:563550/248‑1 (MQ=255)
                                                                                                                                                                        tCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGa                                                                                                                                                                                  >  2:747052/1‑157 (MQ=255)
                                                                                                                                                                         cGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGaa                                                                                                                                                                                 <  1:747053/157‑1 (MQ=255)
                                                                                                                                                                          ggATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGCCCGTTGCGGTACTGGGCAttt                                                               <  1:53761/270‑1 (MQ=255)
                                                                                                                                                                                   gaaatgcGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGTTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGCCCGTTGCGGTACTGGGCATTTACCCTACAAAc                                                    <  2:815883/266‑1 (MQ=255)
                                                                                                                                                                                    gAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTc                                                                                                         >  1:264143/1‑218 (MQ=255)
                                                                                                                                                                                    gAAATCGCAATGGACGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAggg                                                                                                                                                                                   >  1:808006/1‑144 (MQ=255)
                                                                                                                                                                                     aaaTCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGa                                                                                                                                                                                  <  2:808005/144‑1 (MQ=255)
                                                                                                                                                                                     aaaTCGAAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCg                                                                                                        <  2:264143/218‑1 (MQ=255)
                                                                                                                                                                                           cAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTgag                                                                                                                                                                                                                                  >  2:689791/1‑92 (MQ=255)
                                                                                                                                                                                            aaTGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTgaga                                                                                                                                                                                                                                 <  1:689792/92‑1 (MQ=255)
                                                                                                                                                                                                                     ggggCTGTGAATCCGAGTTTAAAAGGGCGTTAGTTTAAATTTTGACGAATCTTGGGGTTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGCCCGTTGCGGTACTGGGCATTTACCCTACAAACTGCTGTCTCACAGGAGCGTGAAGAGAATCGCCTg                  <  2:280185/269‑1 (MQ=255)
                                                                                                                                                                                                                       tGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGGCTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTg                                                                                                                           >  1:861625/1‑165 (MQ=255)
                                                                                                                                                                                                                        gCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGGCTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGt                                                                                                                          <  2:861624/165‑1 (MQ=255)
                                                                                                                                                                                                                         cTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGc                                                                                            >  1:181408/1‑194 (MQ=255)
                                                                                                                                                                                                                          tgtgAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCt                                                                                           <  2:181408/194‑1 (MQ=255)
                                                                                                                                                                                                                           gggAATCCGGGTTAAAAGAGGAGGTATTTTAAATTTTGACTAATCTGGGGTTTCTGTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGTAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGCCCGTTGCGGTACTGGGCATTTACCCTACAAACTGCTGTCTCACAGGAGCGTGAAGAGAATCGCCTGCCGCAc            <  2:642723/270‑1 (MQ=255)
                                                                                                                                                                                                                                      tataAAGGTGGGGTTGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAAGGGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGCCCGTTGCGGTACTGGGCATTTACCCTACAAACTGCTGTCTCACAGGAGCGTGAAGAGAATCGCCTGCCGCACTATGACAATg  <  2:792425/271‑1 (MQ=255)
                                                                                                                                                                                                                                       ataAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTgag                                                                                                                                                                                                                                  >  2:995944/1‑48 (MQ=255)
                                                                                                                                                                                                                                       ataAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTaa                                                                                                               >  1:326156/1‑161 (MQ=255)
                                                                                                                                                                                                                                       ataAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGCCCGTTGCgg                                                                           >  2:387258/1‑197 (MQ=255)
                                                                                                                                                                                                                                        taAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTgaga                                                                                                                                                                                                                                 <  1:995945/48‑1 (MQ=255)
                                                                                                                                                                                                                                        taAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAac                                                                                                              <  2:326155/161‑1 (MQ=255)
                                                                                                                                                                                                                                        taAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGCCCGTTGCGGt                                                                          <  1:387259/197‑1 (MQ=255)
                                                                                                                                                                                                                                                cGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGCCCGTTGCGGTACTggg                                                                    >  1:612948/1‑195 (MQ=255)
                                                                                                                                                                                                                                                cGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGCCCGTTGCGGTACTggg                                                                    <  2:612947/195‑1 (MQ=255)
                                                                                                                                                                                                                                                   tAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCTTCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTg                                                                                                                           >  1:301365/1‑137 (MQ=255)
                                                                                                                                                                                                                                                    aGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCTTCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGt                                                                                                                          <  2:301364/137‑1 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
AACGTAACCGCGGCAGCGTCAGACACAGCCAGTTGTGTCGATTGCGGTTCCACAGGCGCTTCCACTGTGCGGCTTTTTATATATAGCCCGGTAGAGGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAGGTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAGGTGATTATCACCATGCGAATTAACGAAGTTTTTACGGAGGGAAACAATCTCTAGACCATCCTTAACGATTCAGCCACTTTTTTATGTTGCTTTTTTGTAAACAGATTAACACCTCGTCAAAATCCTGCTATTCTGCCCGTTGCGGTACTGGGCATTTACCCTACAAACTGCTGTCTCACAGGAGCGTGAAGAGAATCGCCTGCCGCACTATGACAATG  >  NC_000913/131097‑131596

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: