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,8310CT97.9% 142.2 / NA 47R152C (CGC→TGC) flgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/22);  minor base A (1/0);  total (25/22)
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.77e-01

CCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCA  >  NC_000913/1138592‑1139068
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
ccaaccaGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCGTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGAc                                                                                                                                                                                                                                    <  1:138899/251‑1 (MQ=255)
             tGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAat                                                                                                                                                                                                                        >  2:17285/1‑250 (MQ=255)
              gcggcgCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTTATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAatat                                                                                                                                                                                                                      <  1:592203/251‑1 (MQ=255)
                    cAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTAGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGCTGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGCAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACACACAGGTCAATATCGCGAt                                                                                                                                                                                                                >  2:618800/1‑251 (MQ=255)
                             aGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAg                                                                                                                                                                                                       <  2:492815/251‑1 (MQ=255)
                              ttAGCGGTCTGACTGCCCGCTATGCGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTCCCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCACACGCTGGAGAGTAACGCGGCAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCATCCGGATTGGTGAATCAGTTTCACACCACCGATCCATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGCTATGTGCCAGc                                                                                                                                                                                                      >  2:26703/2‑251 (MQ=255)
                                  cGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTg                                                                                                                                                                                                  >  2:200866/1‑251 (MQ=255)
                                   ggTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTg                                                                                                                                                                                                  >  1:624261/1‑250 (MQ=255)
                                    gTCTGACTGCCCGCTATGAGCAGATGTCGAAAAGCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTgat                                                                                                                                                                                                <  1:603709/251‑1 (MQ=255)
                                     tCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTgatc                                                                                                                                                                                               >  2:328620/1‑251 (MQ=255)
                                               cGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATcaacaa                                                                                                                                                                                     <  1:76381/251‑1 (MQ=255)
                                                   aTGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCACTATCGCGATAGGTGCCAGCGTTgatcagatc                                                                                                                                                                                          >  1:587495/1‑242 (MQ=255)
                                                   aTGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGGACCTCTTCGCTGGCAACACAGATGCAGGATTTCTGCACCAGCCTGCAGACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAAACAGTATAAAACCACCTATCAATATCTGTGCGACCAGGACAAACAGGTCAAGACCGCGATAGGGTCCAGCGTTgatcagatc                                                                                                                                                                                          <  2:587495/242‑1 (MQ=255)
                                                        cAGACGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATATCTTCACCAGCCTCCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGCTTGGGAAATCAGAAGGATTGGTGACTCAGTCTCAACCCACCGATCCATATCTGTGCGACCACGACAAACAGGTCAATACCGCGATAGGCTCCAGCGTTcatcagcaccaccactaccccac                                                                                                                                                                            >  2:343410/1‑234 (MQ=255)
                                                          gATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACt                                                                                                                                                                                   <  1:493923/242‑1 (MQ=255)
                                                          gATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACt                                                                                                                                                                                   >  2:493923/1‑242 (MQ=255)
                                                                       gACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCt                                                                                                                                                             >  1:439114/1‑251 (MQ=255)
                                                                        aCAATATGCCCTCCTCCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGTTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTg                                                                                                                                                            <  1:618800/251‑1 (MQ=255)
                                                                         cAATATGCTCTCCACCAGTACCTCTTCGCTGGCCACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGACAGGAGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTCAACAAATTGCCAGCCTg                                                                                                                                                            >  1:563936/1‑250 (MQ=255)
                                                                               gCTCTCCACCAGTGCCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATAGCTCGATAGGTGCCAGCGTTGAGCAGATCATCAACTACGCTAAACAAAGTTCCAGCCTTAACGAt                                                                                                                                                      <  2:234342/250‑1 (MQ=255)
                                                                                                          aaCACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTggggg                                                                                                                          <  2:186522/251‑1 (MQ=255)
                                                                                                             acaGATGCAGGATTTCTTCACCAGCCTGCACACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAACCAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGc                                                                                                                         >  2:597412/1‑249 (MQ=255)
                                                                                                             acaGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGCGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAAACAGTTTAAAACCACCGATCAATCTCTTTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTCAGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGCAGTGGGGGc                                                                                                                         >  2:200695/1‑249 (MQ=255)
                                                                                                                 aTGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGc                                                                                                                   <  2:437220/251‑1 (MQ=255)
                                                                                                                   gCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCg                                                                                                                  <  1:131100/250‑1 (MQ=255)
                                                                                                                    cAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGACGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTc                                                                                                                >  1:567793/1‑251 (MQ=255)
                                                                                                                                        gCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGAGAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGGCACCTAACAATCTGCTGGATc                                                                                            <  2:492697/251‑1 (MQ=255)
                                                                                                                                           aaCGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAAc                                                                                         <  1:529759/251‑1 (MQ=255)
                                                                                                                                              gCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATctgctg                                                                                                >  2:128795/1‑241 (MQ=255)
                                                                                                                                              gCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATctgctg                                                                                                <  1:128795/241‑1 (MQ=255)
                                                                                                                                                 ggTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATc                                                                                   <  2:351686/251‑1 (MQ=255)
                                                                                                                                                         aCGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTg                                                                           <  1:340908/251‑1 (MQ=255)
                                                                                                                                                          cgcgGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTg                                                                           >  1:395124/1‑250 (MQ=255)
                                                                                                                                                          cgcgGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCATGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGAGATCAACTGGTg                                                                           >  2:575293/1‑250 (MQ=255)
                                                                                                                                                              gAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAACCAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCg                                                                       >  1:456459/1‑250 (MQ=255)
                                                                                                                                                                 gACCCGGCAGCGCGCCAGGCGCTGATTTGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACt                                                                                                                                                                                   <  2:382459/139‑1 (MQ=255)
                                                                                                                                                                 gACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATCGGTGCCAGCGTTGATCAGATCAACAACt                                                                                                                                                                                   >  1:382459/1‑139 (MQ=255)
                                                                                                                                                                  aCCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAAGATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTa                                                                  <  2:361160/251‑1 (MQ=255)
                                                                                                                                                                                            gggAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACGGGTGAGCGAATTAAACCAGATTGTTGGTGGAGAAGTCAg                                        <  2:624261/251‑1 (MQ=255)
                                                                                                                                                                                             ggAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGAGCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGc                                       <  1:369073/251‑1 (MQ=255)
                                                                                                                                                                                                 aTCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTc                                   >  1:456166/1‑251 (MQ=255)
                                                                                                                                                                                                       aGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAAATACGCTAAACAAATTGCCACCCTGAACGATCCAATCTCACGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGCTCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCCGGATg                             >  2:629582/1‑251 (MQ=255)
                                                                                                                                                                                                           atGGTGAATCCGTTTAAAACCACCGATCAATCTCTGAGCGACCAGGACAAACAGGTCAATATCGCGCTAGGTGCCAGCGTTGATCAGATCAACAACCACGATAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGGGAGCGAATTAAACCAGATTGTTGGAGTAGAAGTCAGCGTTCAGGATggagg                         >  2:390680/2‑248 (MQ=255)
                                                                                                                                                                                                                                 ccGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTAAGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGGGAGCGAATTAAACCAGATTGGTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAAAATCACGATGGcc   >  1:445816/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCa  >  1:56416/1‑250 (MQ=255)
                                                                                                                                                                                                                                     tCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATgg                            >  2:58884/1‑222 (MQ=255)
                                                                                                                                                                                                                                     tCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATgg                            <  1:58884/222‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
CCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGCTAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCA  >  NC_000913/1138592‑1139068

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