Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,969,982 C→T A335T (GCG→ACG)  tap ← methyl‑accepting chemotaxis protein Tap

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,969,9820CT94.5% 165.3 / NA 57A335T (GCG→ACG) tapmethyl‑accepting chemotaxis protein Tap
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base T (32/20);  minor base G (1/2);  total (34/23)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGT  >  NC_000913/1969737‑1970225
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
cAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTgcc                                                                                                                                                                                                                                                 >  2:80373/1‑250 (MQ=255)
          aTCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTcgc                                                                                                                                                                                                                                       >  1:382175/1‑251 (MQ=255)
             tctTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGACACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGc                                                                                                                                                                                                                                    >  1:121571/1‑251 (MQ=255)
                 ttGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCAttt                                                                                                                                                                                                                                >  1:14103/1‑251 (MQ=255)
                                   ggCTGGCAAGATTGCGCACTTCACCCGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGAATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCCTGGTACTGACCTAAACACCGCCCGTCTGCGCCGTTGTCGCGGCTTTTTTTGCCAGTTCAGACGcc                                                                                                                                                                                                              >  2:265031/1‑251 (MQ=255)
                                    gCTGGCAAGCTTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCt                                                                                                                                                                                                             <  1:149103/251‑1 (MQ=255)
                                      tGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATCTTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGACGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGGTCCGCCGCCTgg                                                                                                                                                                                                           >  2:117442/1‑250 (MQ=255)
                                       ggCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTc                                                                                                                                                                                                          >  1:183013/1‑251 (MQ=255)
                                        gCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCg                                                                                                                                                                                                         >  1:324544/1‑251 (MQ=255)
                                         cAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGGCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGCCGCCTGTCgg                                                                                                                                                                                                        >  2:240143/1‑250 (MQ=255)
                                         aaaattggctcattcccccTGCCACTCCCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGt                                                                                                                                                                                                        <  1:488239/234‑1 (MQ=255)
                                                cccccctcccccTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTAt                                                                                                                                                                                                 <  1:76173/241‑1 (MQ=255)
                                                          ccTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCACGTCTGCGCCGTTGTTGCGGCATTTTTTGCCAGTTCCGCAGCCTGTCGTGCGTTATCGGCGTTTTg                                                                                                                                                                                       >  2:508261/1‑251 (MQ=255)
                                                           cTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGc                                                                                                                                                                                      >  2:476410/1‑251 (MQ=255)
                                                              ccACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCt                                                                                                                                                                                   >  1:176707/1‑251 (MQ=255)
                                                               cACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCtt                                                                                                                                                                                  >  1:454998/1‑250 (MQ=255)
                                                                   aCCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCgg                                                                                                                                                                              >  1:60002/1‑250 (MQ=255)
                                                                     cGCAAAACAACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATCTTGGGCTGGGAAGCAATTCCGTCAATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGAATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGCCGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTgg                                                                                                                                                                            >  2:489240/1‑251 (MQ=255)
                                                                           cccACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTa                                                                                                                                                                      <  1:132348/250‑1 (MQ=255)
                                                                            ccccGCCCCTGCCTCCCGGCGCGAGCCGCTTCAACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTGTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTaa                                                                                                                                                                     <  2:298476/248‑1 (MQ=255)
                                                                            ccACCCCCTTGCTCTCCGGGCCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTaa                                                                                                                                                                     <  2:176707/251‑1 (MQ=255)
                                                                                  cccTGCCCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTc                                                                                                                                                               <  1:240105/251‑1 (MQ=255)
                                                                                      gCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCAt                                                                                                                                                            <  1:236527/250‑1 (MQ=255)
                                                                                            cGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCCTACTggcg                                                                                                                                                     >  1:195185/1‑251 (MQ=255)
                                                                                             gggctcgAGCCCCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTggcg                                                                                                                                                     <  2:441208/246‑1 (MQ=255)
                                                                                                  cgAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCATTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTTATTGCTCCATATTGGCGTCCGt                                                                                                                                                >  2:254052/1‑251 (MQ=255)
                                                                                                  cgAGCCGCTTCCACTGCCGCATTCAGGGACAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGCTAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCATTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGt                                                                                                                                                >  2:36101/1‑251 (MQ=255)
                                                                                                    aggcgatgccACCGCCGCATTCAGGGCCAGTATAGTGGTCTGGAAAGCATTTCCGTAGATAACGTGGTTAATGTAGCAATTGTGTTGCGAGCAGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTAATTACTTGAACCGCTCCCGGCTTCGCCGTTGTCGTGGCATGGTGTGACAGGTCCTACTCCTGACGTGCGTTATCGGCGTTTTGCCCGACCGTGGCGGTTAATTTCTCCATACTGGCTGCCGttt                                                                                                                                              <  2:206884/243‑1 (MQ=255)
                                                                                                       cGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTtgtg                                                                                                                                           >  1:8495/1‑251 (MQ=255)
                                                                                                             cACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTACTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCGGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTTCCCCCTACTGGAGGGCGttttagccaacc                                                                                                                                     >  2:509576/1‑242 (MQ=255)
                                                                                                                      aTTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCct                                                                                                                            <  1:489240/251‑1 (MQ=255)
                                                                                                                                         ttCTGGAAAGCAAGTCCTTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTTGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGATCACCGCCCGTCTGCGCCGTTGTCGCGGCAtttttt                                                                                                                                                                                                                             <  2:366083/134‑1 (MQ=255)
                                                                                                                                         gTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGGTTGTGCCAACGATGCCGCCTGCTGTTCGGTCCGACTTgg                                                                                                         >  1:535194/1‑250 (MQ=255)
                                                                                                                                         gTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGAGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCAtttttt                                                                                                                                                                                                                             >  1:366083/1‑135 (MQ=255)
                                                                                                                                            tGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGACGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTTTGCCCGCGATGCCCCCTACTGTTccgttcggctttgg                                                                                                        >  1:61411/1‑236 (MQ=255)
                                                                                                                                                             gATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTTTGCCAGCGATGCCCCCTGCTGTTCGGTACGACTTGAGAGATCGTTTTTGCCTGCgc                                                                                     >  1:503356/1‑250 (MQ=255)
                                                                                                                                                             gATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGGTTTTGCCAGCGATTCCGCCTGCTTTTTGGGACGGCCTGTGAGATCGTTTTTGCCTGCgg                                                                                     >  1:194098/1‑250 (MQ=255)
                                                                                                                                                                atctcttaTAATTTCGACAATTTTTTGCGAGCTGGTGGCTATCTCCGGCATGGGGTGAGTCATGGTACTGACCTGAATACCGCCCGTCTGCGCCGGTGTCGCGGTATGTTTTGCCAGGTCCGACGCCTGTCGTGCGTTATCGGCGGTTTGCCCTACCTTGGCGGGTAATTGCTCAATACGGGCGTCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACTACTTGAGAGAGCGTTATTGCCTGAGACa                                                                                   <  2:410329/243‑1 (MQ=255)
                                                                                                                                                                 ccgctgataatggtgtacagttttaGCGAGATGGTGGCGATCTCCTGCATGGTGTGAGTCATGTGTATGTCCTGAACACCGCCTGGCTGCGCCGGTGTCGCGGGCTTTTTTGCCAGTTCCGATGCCTGTTGTGCGTTATCGGTGTTTTGCCCTCCCGTTGCGGTTAATTGCTCCATACTGGTGTCCGTTTGTGCCCGCGATGCCCCCTGCTGATCTGTACGACGTGTGAGATCGTTATTGCCTGGGACaac                                                                                 <  2:196066/231‑2 (MQ=255)
                                                                                                                                                                    cTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGCTAATTGCTCCCTACTTGCGGCCGTTTGGGCCAGCGGTGCCGCCTGCGGTTATGTACGACCTTAGAGCTCGTTTTTGccgtccgcaatctc                                                                              >  2:495184/1‑239 (MQ=255)
                                                                                                                                                                      gATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGCGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTTCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGGACGACTTGAGGGATCGTTATTGCCTGCGCCAATCTccc                                                                            >  1:362206/1‑250 (MQ=255)
                                                                                                                                                                      gATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTcc                                                                             <  1:29427/250‑1 (MQ=255)
                                                                                                                                                                        aaaTGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCGCcgcg                                                                          <  1:145659/250‑1 (MQ=255)
                                                                                                                                                                          atttCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCCTACTGGCGGCCGGTTGTGCCCGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTTCCTTCGGCAATCTcctcgtt                                                                        >  1:62523/4‑246 (MQ=255)
                                                                                                                                                                                ccAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCCGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAAACTCCGCGATAccca                                                                   >  1:53794/1‑248 (MQ=255)
                                                                                                                                                                                  aagttttcGCGAGCTGGTGGCGATCTCCTTCATGGTGTGATTCATGGTCTTCACCTGCGCCCGGCCCGCTGGCCCCGTTGTTCCGTCATGTTTTGCGAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATGGCTCCATACTGGCGGCCGTGTGTGCCAGCGATGCCGCCTGCGGTTCGGTACGATGTGAGAGATCGTTATTGCCTTCGACAATCTCGTCGATACCAAt                                                                  <  2:73022/246‑1 (MQ=255)
                                                                                                                                                                                             gAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTg                                                      <  1:79050/251‑1 (MQ=255)
                                                                                                                                                                                              aGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTgg                                                     <  1:560286/251‑1 (MQ=255)
                                                                                                                                                                                                 tggtggCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTTCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCtt                                                  >  1:110788/1‑251 (MQ=255)
                                                                                                                                                                                                  ggtggCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGGCAATCTCCGCGATACCAATGTGCATCTCCTGGCTTc                                                 >  1:540417/1‑251 (MQ=255)
                                                                                                                                                                                                  gggggCGATCTCCTGCATGGTGCGATTCATGGTACTTGCCTGAACGCCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTACGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGAGCGTTAGTGTCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTg                                                      <  2:382175/243‑1 (MQ=255)
                                                                                                                                                                                                      gCGATCTCCTGCATGTTGTTAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCtt                                             <  1:560099/251‑1 (MQ=255)
                                                                                                                                                                                                           ccccTGCATGGTGTGAGTCATGGTACTGCCCTGAACCCCGCCCGTCTGCGCCGTTGTCGCTGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGca                                       <  2:540417/249‑1 (MQ=255)
                                                                                                                                                                                                                   tGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGCGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACtt                               >  1:395011/1‑251 (MQ=255)
                                                                                                                                                                                                                            gtcaggggaCTGGCCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCATTTCCGACGCCTGGCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGTTTAATTGCTCCATATTTGCGGCCGTTTGTGCCAGCGATGCCTCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCTCATCACTTACCGTccc                       <  2:491369/243‑1 (MQ=255)
                                                                                                                                                                                                                                   tACTGACCTGAACACCGCCCGTCTACGCCGTAGTCGCGGCCTTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCATACCGTGGAGGTTCATAGCTCCATGCTGGCGGCCGTTTGTGCCCGCGCTGCCGCCTGCGTTTCTGCACGAATGTAGAGCTCGTTATTGCCTCCGGCAATCTCCCCGATACAACTGTACCTCCCCCGACTTCCCTTATCAAACACACTTACCGTCCCACGCaa                 >  2:279914/1‑250 (MQ=255)
                                                                                                                                                                                                                                                acCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCCCCTGCTGTTCGGTAAGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCCCATACCTTACCGTCCCACGCAAAGCCCGCTGCATgg   >  1:504084/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 ccgcccggCTGCGCCGTGGCCGCTGCCGTTTTTTCCAGTTCCTACTCCTGTCGTGCTTTATCGGCGTTTGGCCCTTCCGGGGCGGTTAATTGCTCCATCTTTGCGGCCTGTTGAGCCAGCGATGCCGCCTGCTGTTCGGGACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCCCATCACTTACCTGCCCACGCAAAGCCTGCTGCATGGt  <  2:216930/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
CAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGATACCAATGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGT  >  NC_000913/1969737‑1970225

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