Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,134,466 C→T A222V (GCC→GTC)  flgF → flagellar basal‑body rod protein FlgF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,134,4660CT96.9% 193.1 / ‑4.1 64A222V (GCC→GTC) flgFflagellar basal‑body rod protein FlgF
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (26/36);  total (27/37)
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.91e-01

GGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACT  >  NC_000913/1134222‑1134708
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
ggCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGGCAATGCCGTTGCGGCAATGGGCGACATGATTGTCAgcgc                                                                                                                                                                                                                                              >  2:424582/1‑251 (MQ=255)
     aaTCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCgcggg                                                                                                                                                                                                                                         >  2:480070/1‑250 (MQ=255)
      atcccGGCGCTCAAGCCGGGCGTTCCGGCAAAGACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTCCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCgcggcg                                                                                                                                                                                                                                        <  1:90992/247‑1 (MQ=255)
      aTCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCgcggcg                                                                                                                                                                                                                                        >  2:235416/1‑251 (MQ=255)
             cgcTCAATCCGGGCGATCCGGCAAATACGGATGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAACCCCAGGCCACGCGTGGGCCGGTACTGAAGGCAGATCCACCCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCACTGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGttttta                                                                                                                                                                                                                                  >  1:198146/1‑248 (MQ=255)
                 cAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATg                                                                                                                                                                                                                              <  1:312282/250‑1 (MQ=255)
                 cAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGAGCCAACCGTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGc                                                                                                                                                                                                                             <  1:237597/251‑1 (MQ=255)
                  aagccGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACACGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGAATGTCAGCGCGCGGCGTTTTGAAATGCa                                                                                                                                                                                                                            <  1:113677/248‑1 (MQ=255)
                  aaTCCGGGCGATCCGGCAAATACGGTTGCGCCAATAGGGCGTCTTACACTGGTGAAAGCCACGAGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGGGATGTCGTGGGTTCTGGAAGGCAGGAACGTCAATGCCGTTGCGGCAATGAGCGACAGGATTGTCAGCGCGCGGCGTTTTGAAATGc                                                                                                                                                                                                                             >  2:15964/1‑250 (MQ=255)
                     ccGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGc                                                                                                                                                                                                                             >  1:252223/1‑247 (MQ=255)
                        ggCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGAAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGGGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGGAATGCAGATgc                                                                                                                                                                                                                       >  2:281245/1‑249 (MQ=255)
                        ggCGATCCGGCAAATACGGTTGAGCCAGTAGGGCGTCTTAATCTGGTGAAAGCCACGAGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCAACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAGGGCAATAACGTCAAAGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGGGTTTTGGAATGCAGagga                                                                                                                                                                                                                       >  2:358922/1‑247 (MQ=255)
                              ccGGCAAATACGGTTGCTACAGTAGGGCGGCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTGTAAGCGCAGAAACCCAGGCCTCGCGTGGGCCGGTACTGCAGGCAGATCCAACCGTGCGTGAGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGAt                                                                                                                                                                                                                <  1:295874/251‑1 (MQ=255)
                                    aaTACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATcagc                                                                                                                                                                                                            >  1:379780/1‑249 (MQ=255)
                                      tACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGacgg                                                                                                                                                                                                                    >  2:419426/1‑236 (MQ=255)
                                      tACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAgg                                                                                                                                                                                                                    <  1:419426/239‑1 (MQ=255)
                                                ccAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGGCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATa                                                                                                                                                                                               >  2:339172/1‑250 (MQ=255)
                                                  aGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACg                                                                                                                                                                                            <  1:352311/251‑1 (MQ=255)
                                                      gggCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAgg                                                                                                                                                                                        <  2:300302/251‑1 (MQ=255)
                                                       ggCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGc                                                                                                                                                                                       <  2:377899/251‑1 (MQ=255)
                                                         cGTCTGAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGGGACGACGGCAGTTTGCGTTTAAGCGCAGAAACCCAGGCCACGCGAGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCg                                                                                                                                                                                     <  2:209321/251‑1 (MQ=255)
                                                          gTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCgg                                                                                                                                                                                    >  1:310094/1‑250 (MQ=255)
                                                                 aCTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGccaacc                                                                                                                                                                             <  1:410256/251‑1 (MQ=255)
                                                                    ggTGAAAGCCCCGGGCAGCGAAGTGCAGCGCGGTGACTACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATTATAACGCAGGCCGTGCcaaccaac                                                                                                                                                                          <  1:339294/251‑1 (MQ=255)
                                                                    ggTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCcaaccaac                                                                                                                                                                          <  1:480070/251‑1 (MQ=255)
                                                                    ggTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACACAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCCGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGAGTGCCAGCGCGCGGCGTTTTGAAATGCCGCAGAAGGTGATCAGCAGCGCCGATGATAACGCAGGCCGTGCcaaccaac                                                                                                                                                                          >  1:233012/1‑251 (MQ=255)
                                                                       gAAAGCCACGGGCAGCGAACTGCAGCGCGGGGACGACGGCAGTTGTCGTGTAATCGCAGAAACCCAGGCCACGCGTGGGCCGGTTCTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGTCCGTGCCAACCAActg                                                                                                                                                                        <  2:312822/250‑1 (MQ=255)
                                                                            ccACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGa                                                                                                                                                                                                                 <  2:168039/204‑1 (MQ=255)
                                                                            ccACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGa                                                                                                                                                                                                                 >  1:168039/1‑204 (MQ=255)
                                                                                 ggCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAg                                                                                                                                                             <  2:379780/251‑1 (MQ=255)
                                                                                  gCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGt                                                                                                                                                            <  1:134165/251‑1 (MQ=255)
                                                                                   cAGCGAAGGTCAGCGCGGTGACGACGGCATTTTTCGTTGAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGtt                                                                                                                                                           <  2:252223/251‑1 (MQ=255)
                                                                                     gCGAAGTTCAGCGCGTTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCAGTCGATGAGTTaa                                                                                                                                                         <  2:473271/251‑1 (MQ=255)
                                                                                      cGAAGTGCAGCGCGTTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAAt                                                                                                                                                        <  1:453485/251‑1 (MQ=255)
                                                                                          gTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGGCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGaa                                                                                                                                                    <  1:339172/251‑1 (MQ=255)
                                                                                            gCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGaaa                                                                                                                                                   <  1:415537/250‑1 (MQ=255)
                                                                                              aGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATACAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAgg                                                                                                                                                 >  1:367786/1‑250 (MQ=255)
                                                                                                cgcgGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATTATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATa                                                                                                                                              <  2:233531/251‑1 (MQ=255)
                                                                                                cgcgGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGACACGCGAGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGAt                                                                                                                                               >  1:426177/1‑250 (MQ=255)
                                                                                                     tGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGtt                                                                                                                                                                                                                                      >  2:29171/1‑158 (MQ=255)
                                                                                                     tGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGAGCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGtt                                                                                                                                                                                                                                      <  1:29171/158‑1 (MQ=255)
                                                                                                        cgacgGCATTTGTCGTTTAAGCGACGAAACCCAGGCCACGCGTGGGCTGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGGGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGAcga                                                                                                                                      <  2:121022/251‑3 (MQ=255)
                                                                                                            ggCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGtat                                                                                                                                  >  1:130790/1‑251 (MQ=255)
                                                                                                                 ttttCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGtt                                                                                                                             >  1:362763/1‑251 (MQ=255)
                                                                                                                  tttCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGtt                                                                                                                             <  2:362763/250‑1 (MQ=255)
                                                                                                                    tCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGc                                                                                                                           <  1:56514/250‑1 (MQ=255)
                                                                                                                      gTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGGTCTGGAAGGCAGTAAAGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGATTTGAAATGCAGCTGCGGGTGATCAGCAGCGTCGATGAAAACGCAGGCCGTGCCAACACACTGCTGACTATGAGTTAAATGACAGGATACAtacacaataaaagactcccg                                                                                                                        >  2:267599/1‑232 (MQ=255)
                                                                                                                          aaGCGCAGAAACCCAGGCCAAGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGAAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCCGCAGCGTCGATGATAACGCAGGCCGTGCCAAACAACTGCAGTCGGTGGGTTAATTGAAAGGATCAATGAAAAGTCTAAGTTGcacaatgg                                                                                                                    >  2:320101/1‑244 (MQ=255)
                                                                                                                             cgcAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGTTGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATgcgc                                                                                                                  >  1:184967/1‑250 (MQ=255)
                                                                                                                                 gAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGa                                                                                                                                                              <  2:2476/202‑1 (MQ=255)
                                                                                                                                 gAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGa                                                                                                                                                              >  1:2476/1‑202 (MQ=255)
                                                                                                                                      ccAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAAt                                                                                                                                                                                                                               <  2:97836/132‑1 (MQ=255)
                                                                                                                                      ccAGGACACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAAt                                                                                                                                                                                                                               >  1:97836/1‑132 (MQ=255)
                                                                                                                                                                aggcagATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGATGGTGATCAGCAGCGTCGATGAGCACGCAGGCCGCGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACACGTATAAGTTGCCCGCTGCGACAGTTTATAGATTCTTTGGGGACAATCGAAAttttt                                                                              >  2:105937/1‑249 (MQ=255)
                                                                                                                                                                                 cGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATTCCGTTGCGGCAATGAGCGTCATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACtt                                                             <  2:258759/251‑1 (MQ=255)
                                                                                                                                                                                  gtgtGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGta                                                                                                                                   >  1:363923/1‑180 (MQ=255)
                                                                                                                                                                                  gtgtGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTg                                                            <  1:281245/251‑1 (MQ=255)
                                                                                                                                                                                  gtgtGATGGCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGta                                                                                                                                   <  2:363923/180‑1 (MQ=255)
                                                                                                                                                                                                     ggAAGGCAGGAACGTCACTGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATTCAGATGAAGGTGATCAGCAGCGTCGATGATAACTCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAAGTTAGCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGtt                                          <  2:428709/250‑1 (MQ=255)
                                                                                                                                                                                                       aaGGCAGTATAGTCAATGGCGTGTCGGCAATGAGCGACATGATGGTCAGCGTGTGGCGTGTGGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCATGCCTAGCCAACCAACTGCTGTCGATGAGTTAAGTGAAATGTTACATGTCAAGTATAAGTTTTCATATGCGCAAGTTTATCGGGTCTATGGGGGCAAGCTCGTTTTATCTATTTTACGAGCCCTTTGAGGCAGGATaatggtcttcc                                       <  2:199113/251‑10 (MQ=255)
                                                                                                                                                                                                                                ggCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATg              >  1:36159/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gAGCGACAGGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACttt        <  2:250767/251‑1 (MQ=255)
                                                                                                                                                                                                                                         ggACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATTGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTc     <  1:317887/250‑1 (MQ=255)
                                                                                                                                                                                                                                            cATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGAGGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACt  <  2:133532/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
GGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACT  >  NC_000913/1134222‑1134708

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