Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,138,831 C→T R152C (CGC→TGC)  flgK → flagellar hook‑filament junction protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,138,8310CT96.1% 146.5 / ‑3.6 52R152C (CGC→TGC) flgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (29/20);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGT  >  NC_000913/1138587‑1139073
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
tattaCCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGAcc                                                                                                                                                                                                                                              >  1:504055/1‑251 (MQ=255)
    accaaccaGTTACGTGCGGCGCAGCCGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGCCAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGCAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACAGATCAATATCTGTGCGACCAgg                                                                                                                                                                                                                                           >  2:394178/1‑250 (MQ=255)
             ttACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGt                                                                                                                                                                                                                                 >  1:318981/1‑251 (MQ=255)
                    cggcgCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATc                                                                                                                                                                                                                          >  1:7736/1‑251 (MQ=255)
                         cAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGAt                                                                                                                                                                                                                     >  1:509533/1‑251 (MQ=255)
                                        ggTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTg                                                                                                                                                                                                       >  1:186044/1‑250 (MQ=255)
                                          tCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATCTGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGCTGCAGGATTTCTTCACCAGCCTGCAAAAGCTGGAGAGTAACGCGGAAGACCCGGCAGAGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGCATCAGTTTAAAACCACCGATCAATATCTGTTCGACCATGACAAACAGGTCAAAATCGCGATAGGTACCAACGTTgatc                                                                                                                                                                                                    >  1:586994/1‑251 (MQ=255)
                                                   ccGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGCTcaac                                                                                                                                                                                            >  1:573800/1‑250 (MQ=255)
                                                     gCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCaacaac                                                                                                                                                                                         >  1:644873/1‑251 (MQ=255)
                                                          gAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCCAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGTCAAACAGGTCCATATCGCGATAGGTGCCAGCGTTGATCAGATCAAAAACTACGc                                                                                                                                                                                    >  1:443109/1‑251 (MQ=255)
                                                             cAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACACGGCAGCGCGCCAGACGCTGATTGGGAAATCAGAAGGATTGGTGAATAAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGTCAAACAGGTCAATATCGCGATAGGTACCAGCGTAGATCAGACCACCAACTACGCTa                                                                                                                                                                                  >  2:180687/1‑250 (MQ=255)
                                                                aTGTCGAAAATCGACAATATGCTCTCCACCAGTACATCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGAGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAAc                                                                                                                                                                               >  1:492599/1‑250 (MQ=255)
                                                                             aCAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGACGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGACTCAGTTTAGAACCACCGCTCAATATCTGAGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTACACAAATTGCCAGCCTg                                                                                                                                                                 >  2:263769/1‑251 (MQ=255)
                                                                             aCAATATGCTCCCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTTCGCTAAACAAATTGCCAGCCTg                                                                                                                                                                 <  2:644873/251‑1 (MQ=255)
                                                                                      tctcCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCa                                                                                                                                                         >  1:568527/1‑250 (MQ=255)
                                                                                          caccaGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTaacca                                                                                                                                                                              >  2:79927/1‑222 (MQ=255)
                                                                                          caccaGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACa                                                                                                                                                                              <  1:79927/225‑1 (MQ=255)
                                                                                                 aCCTCTTCGCTGGCAACACAGAGGCAGGATTTCTTAACCAGCCTGCAAACGCTGGTGAGTACCGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGATATCAGAAGGATTGGTGACTCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACCGGTCAATATCGCGATCGGGTCCAGCGTTTATAAGATCAACAACTACGCTACACAAATTGCCAGCCTGAACGCTCAATTCTCgcgcat                                                                                                                                             >  2:299499/1‑249 (MQ=255)
                                                                                                            ggCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTgg                                                                                                                                  <  1:531339/251‑1 (MQ=255)
                                                                                                              cAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTgggg                                                                                                                                <  1:299499/251‑1 (MQ=255)
                                                                                                                acacAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGCCAGGCGTGGGGGc                                                                                                                              >  2:117828/1‑251 (MQ=255)
                                                                                                                         cAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGACGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGCTCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTc                                                                                                                     >  1:249371/1‑251 (MQ=255)
                                                                                                                         cAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTc                                                                                                                     <  1:664079/251‑1 (MQ=255)
                                                                                                                                        aGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATctgct                                                                                                      <  1:349957/251‑1 (MQ=255)
                                                                                                                                               aaaCGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCa                                                                                                >  1:413780/1‑250 (MQ=255)
                                                                                                                                                    cTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAAcgcg                                                                                           >  1:43728/1‑250 (MQ=255)
                                                                                                                                                               cgcgGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCATTTTAAAACCCCCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAg                                                                              <  1:531965/251‑1 (MQ=255)
                                                                                                                                                                gcgGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAg                                                                              <  1:347438/251‑1 (MQ=255)
                                                                                                                                                                 cgGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCGAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGc                                                                             <  1:62084/251‑1 (MQ=255)
                                                                                                                                                                     aGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAAt                                                                         <  1:117828/251‑1 (MQ=255)
                                                                                                                                                                           ggCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAAcc                                                                   <  1:525660/251‑1 (MQ=255)
                                                                                                                                                                                cgcgcCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACa                                                                                                                                                                              <  1:413508/139‑1 (MQ=255)
                                                                                                                                                                                cgcgcCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACa                                                                                                                                                                              >  2:413508/1‑139 (MQ=255)
                                                                                                                                                                                  cgcCAGACGCTGAGTGGGTAATCACAAGAATTGGTAAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAAGAAAAACAGGTCCCTATCTCGATAGGTGCCAGCGTTAATCAGATCAACAACTACGCTAAACAAATTTGCAGCCTGCAAGATCAAATCACGCGCCTGACAGGCGAGGAGGCAGGGGCGACACCTAACAATCTGCTGGATCAACACGCGCAACTGGAGAGCGAATAAATCCAGAttgt                                                            >  2:473936/1‑251 (MQ=255)
                                                                                                                                                                                         cgcTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGgtgt                                                      >  1:558626/1‑250 (MQ=255)
                                                                                                                                                                                              aTTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGt                                                >  1:330986/1‑251 (MQ=255)
                                                                                                                                                                                                      aTCAGAAGGATTTGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACGACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTc                                        <  2:321349/251‑1 (MQ=255)
                                                                                                                                                                                                                ttGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACACATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCCCCTAACAATCTGCTGGATCAACGCGATCAACTGGGGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGtgca                                       >  2:377881/1‑239 (MQ=255)
                                                                                                                                                                                                                ttGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCa                                       <  1:377881/242‑1 (MQ=255)
                                                                                                                                                                                                                     gactCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCATCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGCTTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACtt                         >  2:229421/4‑251 (MQ=255)
                                                                                                                                                                                                                       aTCAGATTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACCATCTGCTGGATCAACGCGGTCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGAGggcggaccttat                       >  2:389344/1‑244 (MQ=255)
                                                                                                                                                                                                                         cAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATaa                     <  1:373059/251‑1 (MQ=255)
                                                                                                                                                                                                                          aGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAg                                                 <  1:434595/222‑1 (MQ=255)
                                                                                                                                                                                                                          aGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAg                                                 >  2:434595/1‑222 (MQ=255)
                                                                                                                                                                                                                           gTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTTCGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACa                   <  2:568527/251‑1 (MQ=255)
                                                                                                                                                                                                                               aaaaCCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCAc               >  1:209851/1‑251 (MQ=255)
                                                                                                                                                                                                                               aaaaCCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCAc               >  1:501479/1‑251 (MQ=255)
                                                                                                                                                                                                                                    caccGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTTCCATCGTTGATCAGAGCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCTCTTATAACAGCTCGATgg          <  2:603310/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tccGATCACTATCTGTGCGACCCGAACAAACAGGTAAATATCGCGAGAGGAGCAAGCGAGGAGCAGATCAACAACTACGCTAAACAAATTGACAGCCTGCACGATCAACTCTCGCACCTGCCAGGCGAGAGGGCAGGGGCGGCACCTACAAATCTGCTGGATCAACGCGTTCAACTGGGGAGCGAATTAAACCAGATTGTTGGTGTAGAAGACAGCGTTCAGGAAGGCTGCCCTTTTAACAtaccgagagc         >  2:403321/2‑242 (MQ=255)
                                                                                                                                                                                                                                      ccGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGcc        <  1:602431/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTGGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATg    <  1:263769/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aaTATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAAATACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGTCCCGTATAACATCACGATGGCCAATGGt  <  1:211058/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
TATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGT  >  NC_000913/1138587‑1139073

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: