Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,436,447 G→A 10.4% G414D (GGT→GAT)  rsmB → 16S rRNA m(5)C967 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009133,436,4470G→A10.4% 2.5 0.0 ‑0.144G414D (GGT→GAT) rsmB16S rRNA m(5)C967 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  new base (2/3):  ref base (19/20):  total (21/23)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.33e-01

TTATTAGATGCGCCTTGTTCAGCAACCGGTGTGATTCGTCGCCATCCAGATATTAAATGGTTACGTCGCGATCGCGATATCCCGGAACTCGCGCAATTGCAGTCTGAAATTCTCGACGCCATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGA  >  NC_000913/3436164‑3436747
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                            
ttattaGATGCGCCTTGTTCAGCAACCGGTGTGATTCGTCGCCATCCAGATATTAAATGGTTACGTCGCGATCGCGATATCCCGGCACTCGCGCAATTGCAGTCTGAAATTCTCGACGCCATTTGGTCGCATTTAACAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGACCTTTTCTAAACAGGCACACCAAGACAACCCGGCAAACACACTCTCCCTCGTCGCGAAGCgggagccg                                                                                                                                                                                                                                                                                             >  1:338681/1‑296 (MQ=255)
                   cggcaagcggtggtggtcgctgcaatcaatatatTAAAGGGGTACGGCGTCATAGGGAGATGCCGGAACTAGCGCAATTGCAGTCTGAAATTCTCGACGCCATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaagct                                                                                                                                                                                                                                                                          <  1:14935/271‑1 (MQ=255)
                         ccGGTGTGAGTCGTCGCCATCCAGAGATTAACGGGTTACGTCGCTATCGCGATACCCCGGAACGCGCGCAATTGCAGTCTGAAATTCTCGACGCCATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTTTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaa                                                                                                                                                                                                                                                                    <  2:277120/301‑1 (MQ=255)
                            gtgtGATTCGTCGCCATCCAGATATTAAATGGTTACGTCGCGATCGCGATATCCCGGAACTCGCGCAATTGCAGTCTGAAATTCTCGACGCCATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGCAGAGAATAGCCTGCAGCTTAAAGCCTTTTTGCAACGTACCGCTGATGCCGCACTTTGCGAACCAGCAACACCAGACCAACCAGGTCAACACAATCTACCTGGTGCCGCAAAGTGCGACGGCTTCTTTTACgccactatcaacagaca                                                                                                                                                                                                                                                                 >  1:172868/1‑286 (MQ=255)
                              gtgAGTCGTCGCCAGTCAGATATTACAGGGTGATGGTGCGATCGCGATATCCCGGAAGTAGCGGAAGTGCAGTCGGAAATTCGCGAAGCCATTTGGCCGCATGTAAAAACCGGTGGAACTCTGGTCTATGCTACCTGTTCGGTGTGACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTGTGCAACGTACCCCTGATGCCGAACTTGGCTAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAg                                                                                                                                                                                                                                                                <  2:234447/300‑1 (MQ=255)
                                           aTCCAGATATTAAATGGTTACGTCGCGATCGCGATATCCCGGAACTCGCGCAATTGCAGTCTGAAATTCTCGACGCCATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACACAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAAATAAAgg                                                                                                                                                                                                                                                  >  1:216540/1‑301 (MQ=255)
                                             ccaggaatgaaaGGGTTACGTCGCGATAGCGATATCCCGGAACTCGCGAAATTGCAGTCTGAAATTCTCGACGCCATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCCTGTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGt                                                                                                                                                                                                                                                <  1:276850/292‑1 (MQ=255)
                                                gATATTAAATGGTTACGTCGCGATCGCGATATCCCGGAACTCGCGCAATTGCAGTCTGAAATTCTCGACGCAATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCATTTTTGCAACGTACAGCTGATGCCGAACTTTGCGAAACAGGAACACAAGAAAAACCGGGTAAACACAAGCTACCGGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCCAACAATGGTGAGATAACGGGTcgg                                                                                                                                                                                                                                             >  2:390217/1‑300 (MQ=255)
                                                    ttAAATGGTTACGTCGCGATCGCGATATCCCGGAACTCGCGCAATTGCAGTCTGAAATTCTCGACGCCATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGc                                                                                                                                                                                                                                                                                           <  1:194690/251‑1 (MQ=255)
                                                      aaaTGGTTACGTCGCGATCGCGATATCCCGGAACTCGCGCAATTGCAGTCTGAAATTCTCGACGCCATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGACACAGGAACACCAGAGCAACCGGGTAAACACAATCTACCTGGTGCCGAAACGGGCGACGGCTTCTTTTACGCTACGCTAATCacaccgcgattcgctaacgggcccctacggg                                                                                                                                                                                                                                       >  1:246276/1‑271 (MQ=255)
                                                          ggTTACGTCGCGATCGCGATATCCCGGAACTCGCGCAATTGCAGTCTGAAATTCTCGACGCCATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCCACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGCCGGCTTCttt                                                                                                                                                                                                                                                                                     >  1:75112/1‑251 (MQ=255)
                                                              aCGTCGCGATCGGGCGATCCCGTAACTCGCTCAATTGCAGTCGGCAATGCTCGACGTCATTTGGCCGCATTTAAAAACCGGTGGCACTCTGGTCTATGTCACTTGTCCGTCGTTACCGGAATAGAATAGCCTGCAGATTAAAGTCTCGTCGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGGACACCAGAGCAAATGGGTAAACAAAATCTAACGGGTGCCGAAGAGGGCGACGGCTGCTGTTACGCTAAGCTAATCAACAAGTGATGAGATAACGGGTCGCGACTGATGAAAAt                                                                                                                                                                                                                                <  2:338689/300‑1 (MQ=255)
                                                                        cgcgatATCCCGGAACTCGCGCAATTGCAGTCTGAAATTCTCGACGCCATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCACCCGGGAAAACAACATCTAACTGGTGCCGACGAGGACGACGGCTTCTTTTACGCGAAGCTAATCAAAAAGTGCTTAGATAACGGGTCGCGCCTGAGAAAAATTTCCATTCttg                                                                                                                                                                                                                     >  1:318424/1‑299 (MQ=255)
                                                                                cccGGAACTCGCGCAATTGCAGTCTGAAATTCTCGACGCCATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATACCCTGCAGAGTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACCCCAGAGCAACCGGGTAAACAAAATCTAACTGGTGCCGAAGATGGCGACGGCTTCTTTTACGCTAAGCTACGCAAAAAGTCATGAGATAACGGCTCGCGACTGATGACAATTATCATTATGTGTGCAggg                                                                                                                                                                                                             >  2:280336/1‑300 (MQ=255)
                                                                                               aTTGCAGTCTGAAATTCTCGACGCCATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAAAAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGATGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCCAAAAGTGATGAGATAACGGGTCGCGACTGATGAATATTATCATTCTGGGTGCCGGCCAGGTTggacgcag                                                                                                                                                                                               >  2:151185/1‑294 (MQ=255)
                                                                                                           aaTTGTCGACGCCATTTGGCCGAATTTAAAAACCGGTGGAACTCTGGTGTCGGCCACCTGGTCGGTGTTACCGGAAGAGAATAGGCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAAc                                                                                                                                                                                  <  2:246281/301‑1 (MQ=255)
                                                                                                                  gaagcaaTTTGGCCGCATGTAAAAACCGGTGGAAGTCTGGTCTATGCCACCTGTGCGGTGTGATCGTAAGAGAATATCTTGCAGATTAAACTCTTTTTGCAACTTACCGCTGATGCAGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCGACCTGGTGACGAAGCGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCTACTGAGGAAAATTATCATTCTGGGTGCCAGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGtt                                                                                                                                                                            <  2:299099/294‑1 (MQ=255)
                                                                                                                               cGCATTTAAAAAACGGTGGAACTCTGGTATATGACACCTGTTTGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTGTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGa                                                                                                                                                              <  1:390208/301‑1 (MQ=255)
                                                                                                                                 cATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGCTGTTACCGGCAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCCAACTTTGCGAACCAGGAACACCAGAGCAAACGGGTAAACAAAATCTACATGATGCCGAAGAGGGCGAAGGAATCTTTTACGCTAAGCTCATCAAAAAGTGATGTGATAACGGGTCGCTCATGCTGAAAATTATCATCCTAGATGCCGTCCAGGTAGGAGGCACACTGGCGGAAAACCTAGTTGACGAGAACACCggt                                                                                                                                                             >  2:129101/1‑298 (MQ=255)
                                                                                                                                  aattAAAAACCGGGGGAACTCTGGTCTATGCCACCGGTTCGGTGTGACCGGAAGAGAATAGCCTGAAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGATGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGatat                                                                                                                                                           <  1:151184/299‑1 (MQ=255)
                                                                                                                                   tttAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGAAGATGAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTATGCGAAACAGGAACACCAGAGCAACCGGGTGAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTGTAGGCTAAGCTAAGCAAAAAGTGATGAGATAACGGGTCGCGACTGATTCAAATTATCAATCTGGGTGCGAGCCCGGTTGGCGGGCCACTGGCGGAAAACCTGGTTGGCGACAACACCGATCtt                                                                                                                                                          >  2:358936/1‑301 (MQ=255)
                                                                                                                                     tAAGAAGCGGGGGAAATCTGGTCCACTCCACCTTGTCGGTGTTACCGCAAGAGAATAGACTGCAGATTAAAACCTTTTTGCAACGTACCGCTGATGCCGAACTCTGCGAAACAGGAACAACAGAGCAACCGGGTAAACAAAATCTACCTGGGTCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTAc                                                                                                                                                        <  2:318432/301‑1 (MQ=255)
                                                                                                                                     tAAAACCCGTGGGAACTCTGTGATAGTCGATCTGTTCGTGGTTACCTGAACATAATATCCTGCATATTCAAGCCGTCTTTCAACGTACCGATTATGCCCAACTCGTCGAAACAGGAACATCAGAGCACACTGGTAAACAAAATCTACCTTGTGCCGAATTGGGCTGCTGCTTCTGTTACGCTAAGCTCATCAAAAAGTGATGAGATAACGGGTCGCGACTGCTGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCTAGAACAACGATATTAc                                                                                                                                                        <  2:364533/301‑1 (MQ=255)
                                                                                                                                     aaaaaaCCGGTGGAACTCTGGTCTATGTCACCGTTTCGGTGTTACCGGAAGAGAATATCCTGCAGATTAAAGGCTTTTGGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGGTCAACCGGGTAAACAAAAGCGACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTAc                                                                                                                                                        <  1:358927/300‑1 (MQ=255)
                                                                                                                                       aacaCCGGTGCAACTCTGGGATATGCCAGCTGTTCGGCGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCGTTTTGCAACGTACCGCTGATGCCGAACTGTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCAGAAGAGGGCGACGGCTTCTGTTACGCTAAGCTAATCAAAAAGTGATGAGATAACTGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACtg                                                                                                                                                      <  1:129100/298‑1 (MQ=255)
                                                                                                                                                       tGGTTTAGGACACCTTTGCGGGTTTACCGGAAAGGAATAGCATGCAGATTAAAGCCTTGTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAAGGGGGCGCGACTGAGGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGACCAACGATATTACTGTTGTCGATACCAACgg                                                                                                                                      <  1:280330/301‑1 (MQ=255)
                                                                                                                                                        ggTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGAAGGCTTCTTTTACGCTAAGCTAATCACACAGTGATTAGATAACGGGTCGAGACTGCTGAAAAATATCATTCTGGGTGCCGGCCCGGTTGGCGGCACCCGGGCGGCAAACCTGGTTGGCGAGCACAACCGTATTCGTGTTGTGGATCCCAACgtg                                                                                                                                     >  1:235010/1‑299 (MQ=255)
                                                                                                                                                              ggctgccTGTTAGGTGGGACCGGAAGGGAAGAGCCTGACGATGAAAGGAGTTTTGCGGCGTACCGCGTATGCCGAACGTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAGAATCTACCTGGTGCCGAAGAGGGCGACGGGTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCg                                                                                                                                <  1:47007/295‑1 (MQ=255)
                                                                                                                                                                    cTGTTCTGTGTTACCGGAAGAGAATAGACTGCAGATTAAAGCCTTTTTTCAACGTACCGCTGATGCCGAACGTTTCGAAACAGGAACACCAGAGCAACCGTGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAATAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCCTTCTGGGTGCCTGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACTATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCgg                                                                                                                         <  2:135426/301‑1 (MQ=255)
                                                                                                                                                                                 ccGGAAGAGAATAGCCTGCAGAGTAAAGACTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGACGGCACAGTGGCGTAAAACCTAGTTGGCGAGCACAACGATAGTACTGTTGTCGATACCCACGGGAAGCGGCTGCgtcacggccaggata                                                                                                            >  2:28174/1‑287 (MQ=255)
                                                                                                                                                                                        agaATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGaac                                                                                                                                                                   <  1:435661/239‑1 (MQ=255)
                                                                                                                                                                                        agaATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTAAGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCTACTGATGAAAATTATCATTATGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGTAAACATGTTTGGCGAGcac                                                                                                                                                                   >  2:435671/1‑239 (MQ=255)
                                                                                                                                                                                         gaATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGAAACACCAGAGCAACCGGGTACACCAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGCACAACGATATTACTGTTGTCGATACCAACGGTGATCGTCTGCGGACCTTACAGGATAAATTTgca                                                                                                    >  1:57522/1‑299 (MQ=255)
                                                                                                                                                                                          aaTAGCCTGCAGATGAAAGCCTTTTTGAAACTTACCGCTGATGCCGAACCTTGCGAAACAGGAACACCAGAGCAACCGGTTAAACAAAATCGACCTGATGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTACTCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAACCCTGGTTGGCTAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGAcc                                                                                                   <  2:235015/301‑1 (MQ=255)
                                                                                                                                                                                                     gATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATCTTACTGTTGTCGATACCAACGGTTAGCGTCTGCGGCCCTTACAGGATAAATTTGACCTGCGGGGCGTg                                                                                        >  1:123011/1‑301 (MQ=255)
                                                                                                                                                                                                            tcctttggtcAAAGTACCGATTAGGCCGAAATTTGCGAAAAAGGAAAACAAGAGCAACCGTGTAAACAAAATCTACATGGTGCCGAAGAGGGCGACGGCTTTTTTTACGATAAGATAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGTCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGc                                                                                 <  2:318205/292‑1 (MQ=255)
                                                                                                                                                                                                             catgtttGCAACGTACAGATGAGGACGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGGGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGGGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAgggg                                                                                  <  2:123012/295‑1 (MQ=255)
                                                                                                                                                                                                                         gTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAACGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGCTTGGCGAGAACAACGATATTACTGTTGTCGATCCCCACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTTGGGGTCGTGCCGGGGCAGGTCTCTCATc                                                                     >  1:255859/1‑300 (MQ=255)
                                                                                                                                                                                                                                  aTGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAGAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACCGGCGGAAAACCTTGTTGGCGAGTCCCACGATATTACTGTTATCTATACCAACCGTGAGCGTCTTCGGACCTTACAGGACAAATTTGACCTGTGGGTGGGGCGGGGGCGGGGGACTCATcacagggcttt                                                           >  1:191045/1‑291 (MQ=255)
                                                                                                                                                                                                                                  aTGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGTGCATGGGCCTCATCCACGCGGAtt                                                           >  2:89182/1‑301 (MQ=255)
                                                                                                                                                                                                                                        aaCTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGACCACCTGGTTGGCGAGACCAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGt                                                                                            >  2:364310/1‑262 (MQ=255)
                                                                                                                                                                                                                                                                  acatcGGGTAAACACAACCTACCTGATGGCGAAGAGCGAGACGGCTTCTTTTACGCTAAGCTAAGCAAAAAGTGATGAGAGAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCTGCCAGGTTGGCGGCACACTGGCTGAAAACCTGGTTGGCGAGAACAACGATATGACGGTTGGCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCAGGTGCCGACGACGCCGATATGc                           <  2:216545/297‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACtg                                                                                                                                                      <  2:11380/174‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACtg                                                                                                                                                      >  1:11380/1‑174 (MQ=255)
                                                                                                                                                                                                                                                                                           gTGCCGAAGAGGGCGACGGCTTATTTTACGCTAAGCTAATAAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGATTGGCGGCACACTGGAGGAACACCTGGTTGGCGAGCACATCGATCTAACTGTTTTAGCTACCAACGGTGAGCGTCTGCGGACCTTACAGGATTAATTTGATCTGTGGGTAGTGATGGGGCGTGGCTGTCTTCAACGAGTCTTGAGGGAGTCGGGTTCAGACGCCGCCGGTGTGCTGTTTGTCGGACCCAGTTCAGGTGa  >  2:103845/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                            
TTATTAGATGCGCCTTGTTCAGCAACCGGTGTGATTCGTCGCCATCCAGATATTAAATGGTTACGTCGCGATCGCGATATCCCGGAACTCGCGCAATTGCAGTCTGAAATTCTCGACGCCATTTGGCCGCATTTAAAAACCGGTGGAACTCTGGTCTATGCCACCTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGA  >  NC_000913/3436164‑3436747

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.