Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,439,956 C→T S15N (AGC→AAC)  rplQ ← 50S ribosomal subunit protein L17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,9560CT97.0% 196.6 / NA 67S15N (AGC→AAC) rplQ50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (38/26);  minor base A (2/0);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATAT  >  NC_000913/3439714‑3440200
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
gAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCCGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGCCGACGATTAGCAACGCTATCAGTCTTGGCACGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCCGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTg                                                                                                                                                                                                                                              >  2:97264/1‑251 (MQ=255)
 aGTGTACCCACCGGCCCGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGTCAGAGTCGTCTTGATGATGTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGc                                                                                                                                                                                                                                             <  2:44880/251‑1 (MQ=255)
    gtAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGt                                                                                                                                                                                                                                          >  1:437739/1‑251 (MQ=255)
    gtAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAAGGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGt                                                                                                                                                                                                                                          <  1:225777/251‑1 (MQ=255)
         caccGGCACGGCTCGCGAAAAGCGGGCCCAGTTCGTTTAACAGTTTTGCCCCGGTCTCGTTATCCCGAGTACGGGCGAATGCCAGACGACGATTAGCAACGATATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGAAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGCACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGTTGTTGCGGTTCAGt                                                                                                                                                                                                                                     >  2:592671/1‑251 (MQ=255)
          accGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGtt                                                                                                                                                                                                                                    >  1:400809/1‑251 (MQ=255)
                  ggCTCGCGAAACGCTGGCCCAGTTCGTAAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGCCGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCa                                                                                                                                                                                                                            <  1:6845/251‑1 (MQ=255)
                      cgcgAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGc                                                                                                                                                                                                                                             <  1:123736/230‑1 (MQ=255)
                      cgcgAAACGCGGGACCAGCTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGc                                                                                                                                                                                                                                             >  2:123736/1‑230 (MQ=255)
                         gAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTAc                                                                                                                                                                                                                     >  1:574910/1‑251 (MQ=255)
                          taaCGCGGGCCCTGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGAGTAGCAACGCTATCAGTCTTGGCAAGAGTAATTAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTGGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTTACGTCCACTCTTa                                                                                                                                                                                                                      <  1:522002/248‑1 (MQ=255)
                           aaCGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCCGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGAACACTCTTACg                                                                                                                                                                                                                    >  2:105201/1‑250 (MQ=255)
                               ggggCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCTCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGgc                                                                                                                                                                                                               <  2:627434/250‑1 (MQ=255)
                                  gCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCa                                                                                                                                                                                                            <  1:620972/251‑1 (MQ=255)
                                                  aGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCATTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTtctc                                                                                                                                                                                            <  2:502010/251‑1 (MQ=255)
                                                   gTTTTGCCACGATCTCGTAAGCACGAGTACGGGCGAATGCCAGACGCCGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTCCGCGGCGAAGCTCTTTCGCTTAAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTTACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCa                                                                                                                                                                                           <  2:437739/251‑1 (MQ=255)
                                                      tgtccACGGTCTCGTTATCACGATAACGGCGGAGTGCCAGACGCCGATTAGCACCCGTATCAGTCTTGGAAAGATTAATCAGCGGCTCAACAACGCGGCCCAGCTCTTTCGCTTTAGGCCGAGTCGTCTTGATGATTTGATGACGAACCAGTGATCCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCCTGACCTTATCCTGCTCAgg                                                                                                                                                                                         <  2:353589/247‑2 (MQ=255)
                                                             gATCTCGTTATCACGAGTACGGGCGACTACCAGACGACGATTAGCAACGATATCAGTCTTGGCAAGCGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTCCGCTTTAGGCAGAGACGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGATGCGGTTCAGTTGACGACCACTCTTACGCTGGCGCATGACCTTATCCTTCTCAGTAAAACCtt                                                                                                                                                                                 >  2:596851/1‑251 (MQ=255)
                                                              agttctttaccACGAGTACGGGGGAATGCCAGACTACGATAAGCATCGCTATCAGTCTTGGCAAGGGTAAGCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGAGTTCATGACGAACAAGGGAACCTGCCATATTGCGGAACAGAGCTTGGCGATGGTTGCTGTTGCGGTTCAGTTTACCACCGCTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTa                                                                                                                                                                                <  1:358980/241‑1 (MQ=255)
                                                                  cGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCACCGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTGGATGATTTCATGACTAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCt                                                                                                                                                                                  <  2:228494/245‑1 (MQ=255)
                                                                     tATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGTCGCATGACCTTATCCTTCTCAGTACAACCTTAACCtgtg                                                                                                                                                                         >  2:532959/1‑251 (MQ=255)
                                                                     tATCACGAGTACGGGCGAAGGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGTCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGTAGAGGCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTTCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGGTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATTCTTCTCAGTAAAACCGTAACCtgt                                                                                                                                                                          <  2:296084/250‑1 (MQ=255)
                                                                        cACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATc                                                                                                                                                                      >  1:368294/1‑251 (MQ=255)
                                                                            aGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGACGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGTTCCGGt                                                                                                                                                                  >  2:511939/1‑251 (MQ=255)
                                                                                cGGGCGAGTGCCAGACGACGGTTAGCACCGCTCTCAGTCTTGGCAAGTGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCTTGACGAACCCGTGAACCTGCCCTATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCCGTTGACGACCACTCTTACGATGGCGCTTGACCTTATCCTTCTCAGTAAAACCTTCACCTGTGCTCCGATTACt                                                                                                                                                              >  2:35622/1‑251 (MQ=255)
                                                                                  ggCGAATGCCAGACGACGATTAGCAACGCTAACAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGGCTTGATGATTTCATGGCGAACCAGTGACCCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCg                                                                                                                                                            <  2:528759/251‑1 (MQ=255)
                                                                                           cAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATg                                                                                                                                                   <  1:592671/251‑1 (MQ=255)
                                                                                              acgacgATTAGCAACGCTAGCAGTCTTGGCAAGAGTAATCAGGGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCtt                                                                                                                                                <  2:528655/251‑1 (MQ=255)
                                                                                                 acgaTTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGACCCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCGTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGcc                                                                                                                                             <  2:400809/251‑1 (MQ=255)
                                                                                                 acgaTTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGcc                                                                                                                                             >  2:468463/1‑251 (MQ=255)
                                                                                                    aTTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCggg                                                                                                                                          >  2:292691/1‑250 (MQ=255)
                                                                                                        gCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGcc                                                                                                                                      >  1:462213/1‑251 (MQ=255)
                                                                                                         cAACGCTATCAGTCTTGGCAAGAGCAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGcc                                                                                                                                      >  2:135567/1‑250 (MQ=255)
                                                                                                           aCGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGCTGCTTGCCGGTGGCCAGt                                                                                                                                   >  2:471998/1‑251 (MQ=255)
                                                                                                                aTCATTCTTGGCAAGAGTAATCATCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGAGGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTcc                                                                                                                              <  2:624800/251‑1 (MQ=255)
                                                                                                                aTCAGTTTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTcc                                                                                                                              <  2:171974/251‑1 (MQ=255)
                                                                                                                   aGTCTTGGCAAGAGTAATCAGCGGCTCAACTCCGCGGCACAGCTCTTTCGCTTTAGGCAGAGTCTTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGACCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAgg                                                                                                                           <  2:588368/251‑1 (MQ=255)
                                                                                                                      cTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGAGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGGCCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGgcgc                                                                                                                        >  1:562202/1‑251 (MQ=255)
                                                                                                                      ataggaaaGAGTAGTCAACGGCTCAACTACGCGCCGCAGCTCTTGCGCTTTCGGCACAGTCGTCTTGATGATATCATGACGCACCCGTGAACCTGCCATATTGCGGAACCTAGCCTGGCGCTGGTTGCTGTTGCGGTTCAGTTGACGAAAACTCTTCCGCTGCCGCATGACCTTATCCTTCTCAGTAAACAATTAACCTGTGATCCGGTTACTCGTCAGCGCTGCTTGCCGACGGCCAGTGTTACAGgcgc                                                                                                                        >  2:47506/7‑251 (MQ=255)
                                                                                                                                   aTCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGAc                                                                                                            >  1:205632/1‑250 (MQ=255)
                                                                                                                                        cGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTCGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGCTCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCAAGGCGCATGCCCAGAGAccgtac                                                                                                       >  2:375781/1‑245 (MQ=255)
                                                                                                                                                      agcgcAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCATTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGACGGTGGCACGTTTTCCAGGCGCCTGCCAAGAGACAGTCCACGGGAAGCCAGctc                                                                                        >  2:498105/2‑249 (MQ=255)
                                                                                                                                                       gcgcAGCTCTTTCGTGTTAGCCAGAGTCGTCTTGATGATTCCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCATTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTGCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACg                                                                                       <  2:564583/251‑1 (MQ=255)
                                                                                                                                                         gcAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGCTGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCAGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCCGCACGTc                                                                                     >  1:381975/1‑251 (MQ=255)
                                                                                                                                                            gttctTTCGCTTTAGTCGGAGGCTGCTTGATATTTTCATGACGCCCCAGTAAACCTGCCATAGTGCGGAACATAGCCTGACGATTGTTGCTGTTTCTGTGCAGTGGACGACCACTCTTAAGAGGGCTCATGACCGTATCCGTCTCAGTAAAACCTTAACCTGTGATCCGGTTACGCCTCAGCGATGTATGCCGGCTGCCAGTATTCCAGGCGAATGCCCAGAGACAATACACGGGAAGCCAGCACGTCttt                                                                                  <  2:241637/249‑1 (MQ=255)
                                                                                                                                                                    gCTTTAGGCAGAGCCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAg                                                                          <  1:321290/251‑1 (MQ=255)
                                                                                                                                                                      tttAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGCTGGTTGCTGTTGCGGTTCAGTTGACAACCACTCTTACGATGGCGCATGACCTTCTACTTCTCAGTAAAACCTTAACATGTGATCCGGTTAATCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGACGCATGCCCAGAGACAGTCCACGGGAAACCAGCACGTCTTTAATCTCAga                                                                         >  2:331351/1‑249 (MQ=255)
                                                                                                                                                                            cAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGGCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAt                                                                  <  2:576666/251‑1 (MQ=255)
                                                                                                                                                                             agagTCGTCTTGATGATTTCATAACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACATGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAtt                                                                 >  1:115574/1‑251 (MQ=255)
                                                                                                                                                                                         atgatTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGt                                                     <  1:292691/251‑1 (MQ=255)
                                                                                                                                                                                         atgatTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGt                                                     >  2:138284/1‑251 (MQ=255)
                                                                                                                                                                                                cATGACGCACCAGTGAAACCGCCATATAGCGGAACATAGACTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATCGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTACCCGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTACGAGATTTTTTACCAAGGTTAGGCGt                                              >  2:32079/1‑251 (MQ=255)
                                                                                                                                                                                                 aTGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCgat                                             >  2:534119/1‑249 (MQ=255)
                                                                                                                                                                                                     cGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGAAGCTGATGCGGTTCAGTTGACGACCACTATTACGAAGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGCTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGCCAGTCCACGGGAAGCCAGAACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGCAGTTTTa                                          >  1:266528/1‑250 (MQ=255)
                                                                                                                                                                                                           agtgagtctatcattttgagggaCATAGCCTGGCGGTGGTTACTGATGCGGTTCAGTTGACAACCACTCATACGCTGGCGCATGAGCTTACCCTTCTCAGAGCAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTACCGACGCATGACCAGCGACCGTCCACGGGACGCCAGCACGTCTTTAATCTCAGTAAGCGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCt                                   >  2:591120/23‑251 (MQ=255)
                                                                                                                                                                                                                aCCTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGAGCATGCCCAGAGACAGTCCACGTGAAGCCAGCACGTCTTTAATCTCAGTCAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAAcc                              >  2:232387/1‑251 (MQ=255)
                                                                                                                                                                                                                 ccTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGGGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCt                             >  1:27445/1‑251 (MQ=255)
                                                                                                                                                                                                                 ccTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCt                             >  1:152244/1‑251 (MQ=255)
                                                                                                                                                                                                                 ccTGCCATATTGCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCt                             <  1:490951/251‑1 (MQ=255)
                                                                                                                                                                                                                 ccTGCCATATTGCGGAACATAGCCTGGCGATGGATGCTGTTGCGGTTCAGTTGACGACAACTCTTACGATGGCGCATGACATTCTCCTTCTCAGTAAAACCTTAACCTGTGATACGGTTACTCGTCAGCGATGCTTGACGGAGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCCCGGGAAGCCAGCACGTCTTTACTCTCAGTAAGAGCTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCt                             >  2:627233/1‑251 (MQ=255)
                                                                                                                                                                                                                        tatTGCGGAACATAGCCTGGCGATGGTTGCAGTTGCGGTTCAGTTGACGACAACTCTTACGATGGCGCCTGACCTTATCATTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCCGGCGCATGCCCAGAGAAAGTCCACGGGAAGCCAGCACGTCTTTCATCTAAGTAAGTGATTTTTTACACAGGTTAGGCGGTTTAAGGAGCTCAACCTCGGTACg                      >  2:77616/1‑251 (MQ=255)
                                                                                                                                                                                                                            gCGGAACATAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGt                  >  1:586947/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tAGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTCATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCCACCTCGGTACGCTGTACCATATc          >  2:277514/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aGCCTGGCGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCAATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCATAGACAGTCCACGGCAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGTTTAGGCGTTTTAAGGATCTCAACCTCGGTACTCTGTACCAGATc          >  2:444696/1‑250 (MQ=255)
                                                                                                                                                                                                                                          ggggATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGat    <  2:14710/248‑1 (MQ=255)
                                                                                                                                                                                                                                            cGATGGTTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGatat  >  2:509949/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cGATGGTTGCTGGTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCCTGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGttat  >  2:319181/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
GAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATAT  >  NC_000913/3439714‑3440200

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑