Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,151,563 A→G 11.3% S140S (TCT→TCC exbB ← membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,151,5630AG11.3% 190.6 / 2.8 78S140S (TCT→TCCexbBmembrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base A (38/26);  new base G (8/2);  total (50/28)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTGCAACATCACCCAGCATCGCTTTAAAGCCGCCAATCTGGCGTGCAAATACGTTATAGATAACGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAAC  >  NC_000913/3151341‑3151786
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
cgcTGCAACATCACCCAGCATCGCTTTAAAGCCGCCAATCTGGCGTGCAAATACGTTATAGATAACGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTCGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAACCCGTCCCAAACACACCAACAAACGGAGAAAta                                                                                                                                                                                                                           >  2:674224/1‑228 (MQ=255)
     cAACATCACCCAGCATCGCTTTAAAGCCGCCAATCTGGCGTGCAAATACGTTATAGATAACGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATCCCCCGCGCAACGACTGCCAGGTTAGTGGTCTGCGTTTGTTTGCGCGATACCCATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAAACAAAAAAACGGGGaaaaac                                                                                                                                                                                                                          >  2:70159/1‑226 (MQ=255)
             cccAGCATCGCTTTAAAGCCGCCAATCTGGCGTGCAAATACGTTATAGATAACGACCGCAGGAATCGCTGCCACGCGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGTTGCCCCAGCCCGTACCAATCGAACCAACGAACTGGAAAATCCGGCCGACGGtt                                                                                                                                                                                                               >  1:719952/1‑228 (MQ=255)
                aGCATCGCTTTAAAGCCGCCAATCTGGCGTGCAAATACGTTATAGATAACGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCACTAAAGCTGTTCAGGATGCCCCAGACCGTACCAAACAAACAAACGAACGGGGAAACATAGCCGATGGtgtca                                                                                                                                                                                                            >  1:396535/1‑224 (MQ=255)
                       cTTTAAAGCCGCCAATCTGGCGTGCAAATACGTTATAGATAACGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGCTCCCAATAAAGCTGTTCATGATGCCCCCGACCGTACCCAACAGACCAACGAAAGGGGAAAACGCCCCCAAGGTTGCCaaaaaaa                                                                                                                                                                                                     >  2:312723/1‑222 (MQ=255)
                       cTTTAAAGCCGCCAATCTGGCGAGCAAATACGTTATAGATAACGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCaactaac                                                                                                                                                                                                     >  1:165652/1‑222 (MQ=255)
                        tttAAAGCCGCCAATCTGGCGTGCAAATACGTTATAGATAACGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACAAACGGCGAAATCGCGACGCTGGTTGCAAGaagacc                                                                                                                                                                                                    >  2:501793/1‑223 (MQ=255)
                            aaGCCGCCAATCTGGCGTGCAAATACGTTATAGATAACGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTCCCAAACAGACCAACGACCGGCGAAATTCGTCCGATGGTTGCCCGATAGCAGTCAc                                                                                                                                                                                               >  1:721297/1‑229 (MQ=255)
                                  ccAATCTGGCGTGCAAATACGTTATAGATAACGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGATAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGc                                                                                                                                                                                         >  1:620528/1‑229 (MQ=255)
                                     aTCTGGCGTGCAAATACGTTATAGATAACGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCa                                                                                                                                                                                      <  1:560473/229‑1 (MQ=255)
                                     aTCTGGCGGGCAAATACGTTATAGATAACGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCa                                                                                                                                                                                      <  1:659924/229‑1 (MQ=255)
                                         ggcagacAAATAGATTATAGATAAGGACCGAAGGAATCGCTGCCAGGAGGCCGATTGACGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAATGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCAtt                                                                                                                                                                                    <  1:460257/221‑1 (MQ=255)
                                                               aaCGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGCGGTCAGC‑‑‑‑GCCTGCGCGAAACCAATACAGCTGTGCATAATTCCCCCCACCGAACCAAACAGGCCAACGAACGGCGAAATGCCGCCGTTGGTTACAAGATATCCGTAACTGCGGCCCATTTGACGACCCACTGGGGGGAcaaaa                                                                                                                                                            >  2:71933/1‑225 (MQ=255)
                                                                 cGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACAACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTACGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGGGAAATCGCGCCGATGGATGCCAAAAACCCGTTACCGAGGCCCATATGACGACCCAAAGCGGCGACCCGACg                                                                                                                                                          >  2:323639/1‑229 (MQ=255)
                                                                  gACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCCAACAGACCAACGAACGGAGAAATCCCGCCGAGGGTTGCCAAATAGCCGTTCCCGAGGCCCATTTCTCAACCCACAgcggcccgccaacga                                                                                                                                                         >  2:735952/1‑219 (MQ=255)
                                                                       cAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCCATAAAGCTGTTCATGATGCCCCAGACCGAACCAAACAGACCAACGAACGGAAAAATCGCGCCCATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACCGCCCCGGCCCGAcccacacg                                                                                                                                                    >  2:191167/1‑222 (MQ=255)
                                                                       cAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAACGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGGTGCCAGATAGCCGGTCCCGCGGCCCACTTGACGACCCCCTGCGGCGACCCGACCCCCCac                                                                                                                                                    >  2:550163/1‑228 (MQ=255)
                                                                       cAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAAACCAACGAACGGGGAAATCGCGCCGCTGGTTGCCAGATAGCCGTTACCCCGGCCCATTTGACGACCCACTgcgggccccccaccccccaa                                                                                                                                                    >  2:723726/1‑213 (MQ=255)
                                                                            accGCTGCCCCGATGCCGATTGCCGTTGCTAACGGGACTTCTGGGATTCCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAAGAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGa                                                                                                                                               <  1:143989/227‑1 (MQ=255)
                                                                               gCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGCGTTTGTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACg                                                                                                                                                 <  1:70159/228‑1 (MQ=255)
                                                                                      cGAGGCCGGTTGCCCTTGCTAAAGAAGCTTCTGCGATCCCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGGGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGCCCCGACGCTCCAGACGGAAGCTGGTACg                                                                                                                                     <  2:702422/229‑1 (MQ=255)
                                                                                      cGAGGCCGATTGCCGTTGCTAACAGTGCCTCTGGGATACCCGGCGCAACGACTGCCAGGGTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTGACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACg                                                                                                                                     <  1:90268/229‑1 (MQ=255)
                                                                                       gagggcgtgtGCCGTTTCTAACAGGGCTTTTGCGATTCCCGGCGCAACGACTGCCAGGTTAGTGTTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGCCCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGTCCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGt                                                                                                                                    <  2:186677/220‑1 (MQ=255)
                                                                                         ggCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGCTACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCCATGGTTGTCAGATAGCCGTTCCCCCGGCCCATTTGCCGACCCACTACGGCTCACCCACCCTCCAGACGGAACACGGTACGtt                                                                                                                                   >  2:143264/1‑228 (MQ=255)
                                                                                                                gCTTCAGCGATCCCCGGCGCAACGACGTCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGc                                                                                                           <  1:314390/229‑1 (MQ=255)
                                                                                                                 cTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGCTACCAATACAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGCGAAATCGCGCCGCTGGTTGCCAGATAGCCGTTACCGCGCCCCATTTGACGCCCCCCCGCGCCGACCCGCCGCTCCATCCGGAACCTGGGACGCTCTTTAATACCTCCGTTATCGTctct                                                                                                          >  2:671157/1‑229 (MQ=255)
                                                                                                                 cTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGCTAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCt                                                                                                          >  1:410654/1‑229 (MQ=255)
                                                                                                                  ttCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGCACGGAGAAATCGCGCCGATGTTTGCCAGAAAGCCGGTACCGCGGCCACTTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGGGCGCTCTTTAAAACATCCGtatacgccgcgc                                                                                                         >  2:383809/1‑218 (MQ=255)
                                                                                                                         aTACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGAAGACCCACTGCGGCGACCCGACCCCCCAACCCGAAGCGGGTAAGTTCTTTATAACCTTCGTTATCGTCCCCGCCTTCt                                                                                                   >  2:538426/1‑228 (MQ=255)
                                                                                                                            cccGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCCTTTTACGCCCCACTGCGGCGACCCGCCCTCCCAGCCGGAAACTGGTACGTTTTTATATACCTCCGTTTATGTCGCTGCCtctcgacc                                                                                               >  2:314390/1‑222 (MQ=255)
                                                                                                                               ggCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAACCTGGTACGTTCTTTAATACCTTCGTTATCATACTGGCCTTCTGAccag                                                                                             >  2:211359/1‑225 (MQ=255)
                                                                                                                               ggCGAAACGACTGCCGGGTTAGTGGTATGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGGCCAACGAACGGAGAAAGCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGtt                                                                                            <  1:211359/229‑1 (MQ=255)
                                                                                                                                 cgcAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGGGAAATCGCGCCGATGGCTGCCAGATAGCCGTTACCGCGGCCAATTTGACGAACCACTGCGGCGCCCCGACGCTCCCGACGGAAGACGGTACGTTCTTTAATACTTTCGTTATCGTCGCTGCCTCTTGACAgg                                                                                             >  2:598866/1‑225 (MQ=255)
                                                                                                                                         cTGCCCGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGAGACCAATAAAGCTGTTCAGGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCAtt                                                                                  <  1:674225/229‑1 (MQ=255)
                                                                                                                                         cTGCCAGGTTAGTGGTCGGG‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCAtt                                                                                  <  1:549616/229‑1 (MQ=255)
                                                                                                                                              aGGTTAGTTGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCCGTTCATGCTGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGGCGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTgcg                                                                             <  1:109436/229‑1 (MQ=255)
                                                                                                                                                    gTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCAGGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCAt                                                                       <  1:501793/229‑1 (MQ=255)
                                                                                                                                                    gTGGCCTGC‑‑‑‑GTTTGCGCGATACCAAAAAAGGTGTTCAGGAGGCCCCAGACCGTCCCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCa                                                                        <  2:660206/228‑1 (MQ=255)
                                                                                                                                                     gggTCTAA‑‑‑‑GTTTGCGCGAGACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCCCTGCGGCGACCCGACGCTCCAGACGGAAGCTGGACCTTTCTTTAAGCCCTTCGTTATCGTCCCGGCCTTCTCACAGTTCCAGCTCATTcatgccttcagt                                                                      >  2:588241/2‑218 (MQ=255)
                                                                                                                                                       gTCTGC‑‑‑‑GTTTGCGCGATACCAATAATGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGCGAACCCACTGCGGCGACCAGACCCCCCCACCGGACGCTGGTACGTTCTTTAATCCCTCCGTTATCGTCGCTGCCTTTCGAAAGTTCACACTCATTCTCGCCTTTATTGa                                                                    >  2:180049/1‑229 (MQ=255)
                                                                                                                                                          tGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACAGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTCCGTTCTTTAATACACTCGTTCTAGGCACTGACTTCCGGCAGTTCCCTCACCTTCTGCGCTTCATTggtga                                                                 >  2:708884/1‑225 (MQ=255)
                                                                                                                                                           gC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGc                                                                  <  1:476887/227‑1 (MQ=255)
                                                                                                                                                                                 aaaaaccgTTCATGATGCCCCAGACCGTCCCAAACAGCCCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTCCCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCtttt                                              <  1:30452/222‑1 (MQ=255)
                                                                                                                                                                                   aaGCTGTTCATGATGCCCCAGGCCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTg                                             <  1:171735/228‑1 (MQ=255)
                                                                                                                                                                                     gCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGtt                                                                                            >  1:374505/1‑179 (MQ=255)
                                                                                                                                                                                     gCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGGGCAAATGCAGGCTTAAGCTTTTGCt                                           >  1:203389/1‑228 (MQ=255)
                                                                                                                                                                                        gTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGACGCTGGTACGTTCTTTAATACCTTCGTTATCGTCCCTGCCTTCTGCAAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATCACGGCTTAAGCTTTTGCTAcaa                                       >  2:155290/1‑227 (MQ=255)
                                                                                                                                                                                        gTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTTAGCAAATGCAGGCTTAAGCTTTTGCTACCa                                       >  1:34410/1‑229 (MQ=255)
                                                                                                                                                                                        gTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTAcc                                        >  1:168198/1‑228 (MQ=255)
                                                                                                                                                                                        gTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCCACCCGGCCCTCCAGAAGGCAGCTGGTAAGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACCGTTCCAGATCATTATTGGCTGCATTTGGCAAATGCAGGCTTACGCCTTTGCTACCa                                       >  2:410028/1‑229 (MQ=255)
                                                                                                                                                                                         ttCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGGCCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTTCCaa                                      >  1:17114/1‑229 (MQ=255)
                                                                                                                                                                                         ttCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCaa                                      >  1:20764/1‑229 (MQ=255)
                                                                                                                                                                                         ttCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCaa                                      >  1:11610/1‑229 (MQ=255)
                                                                                                                                                                                         ttCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGCCCCATTTGACGACCCACTGCGGCGACCCTCCGCTCCAGACGGAAGCTGGGACGTTCTTTTATACCTTCGTTATCGTCCCGCCCTTCTGACAGTACCAGCACATTCTGCCCTCCATTGACCAAAACCAGGGTTAATCGTTTGCACCCaa                                      >  2:417786/1‑229 (MQ=255)
                                                                                                                                                                                           cATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCaaaa                                    >  1:254821/1‑229 (MQ=255)
                                                                                                                                                                                              gatgCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGCCGACCCACTGCGGCGACCCGACGCTCCAGCCGGAACCTGGTACGTTCTTTACTACCTTCGTTTTCGTCGCGGCCTTCTGCCAGTTCCCGCTCATTTCGCGCTTCATTGCGCAAATGCAGCTTTAACCTTTTCTCCCaaaaacc                                 >  2:133426/1‑225 (MQ=255)
                                                                                                                                                                                                   cccAGACCGTACCAAACAGACCAACGAACGGGGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTGTCGGCCTTCTGACAGTTCCAGCTCATTCTGTGCTCAATTGAGACAATCCAGGCTTAAGCTTTTGCTACCAAAATCacgga                             >  2:2864/1‑224 (MQ=255)
                                                                                                                                                                                                   cccAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGCCTTAAGCTTTTGCTACCACAATCAgcggc                            >  1:380603/1‑229 (MQ=255)
                                                                                                                                                                                                   cccAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAACACCTTCGTTATCGTCGCTGCCTTCTGACCGTTCCAGCTCATTCTGCGCTTCATTCAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAccggc                            >  2:190340/1‑229 (MQ=255)
                                                                                                                                                                                                       gACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTaa                                                                                                                             <  2:230152/128‑1 (MQ=255)
                                                                                                                                                                                                       gACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTaa                                                                                                                             >  1:230152/1‑128 (MQ=255)
                                                                                                                                                                                                           gTACCAAACTGACCAACGAACGGAGAAACCGCGCCGATGGTTGCCAGGTAGCCGTTACCGCGGACCATTTGACGTCCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATCCCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGt                    <  1:323639/229‑1 (MQ=255)
                                                                                                                                                                                                            tACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCCCCTTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCCACCAAAATCCACGGGGATATCGtg                   >  2:116613/1‑228 (MQ=255)
                                                                                                                                                                                                            tACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCAGCCTTCTGACAGTTCGAGCTCATTCTCCGCTTCATTGGTCAAATTCAGGCTTAAGCTTTTGCTACCAAAATCAGCCGCGATTACGtt                   >  1:720176/1‑229 (MQ=255)
                                                                                                                                                                                                            tACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGATGCCAGATAGCCGTTACCGCGACCCATATGACGACCCACTGCGGCGACCCGACGCCCCCACCGGAAGCAGGTACGTGCTTAACTACCTTCGTTATCGTCCCTGCCTTCTACCAGGTCCAGCCCATTCTACACTTCATTGTGCAAAGGCAGGCTTAAGGTTTGCCTACCACAACAAGCGCGGATATAGtt                   >  2:154828/1‑229 (MQ=255)
                                                                                                                                                                                                             aCCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAACATCAGCGGCGATATCGTTg                  >  1:561918/1‑229 (MQ=255)
                                                                                                                                                                                                               caaccgcgccAACGCACGGAGAAAGCGCGCCGATTGTTGCCAGATAGCCGTTACCGCGGCCCATGTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATCCCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGc                <  2:56217/221‑1 (MQ=255)
                                                                                                                                                                                                               cAAACAGGACAACGAACGGAGAAATCGCGCCGATTGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGc                <  1:143264/229‑1 (MQ=255)
                                                                                                                                                                                                               cAAACAGACCAACGACCGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCACATTTGACGGCCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCGTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGc                <  1:417786/229‑1 (MQ=255)
                                                                                                                                                                                                               cAAACAGACCAACGAACGGAGAACCCGCCCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGc                <  2:203389/229‑1 (MQ=255)
                                                                                                                                                                                                               aaaaaagccaacggcacgGGGAATTCGCGCAGAGGGGTGCCGGATAGCGCTTACCTCGGCCAATTTGACGACCCACGGCGGCGACCCGCCGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCGTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGc                <  2:69100/214‑1 (MQ=255)
                                                                                                                                                                                                                 aaCAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCt              >  1:122778/1‑229 (MQ=255)
                                                                                                                                                                                                                      gccAACGAACGGGGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTGCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGttt         <  1:2864/228‑1 (MQ=255)
                                                                                                                                                                                                                      aCCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTaa                                                                                                                             <  2:657329/113‑1 (MQ=255)
                                                                                                                                                                                                                      aCCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTaa                                                                                                                             >  1:657330/1‑113 (MQ=255)
                                                                                                                                                                                                                      aCCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATGTTTGGCCTGGtt          >  1:588378/1‑228 (MQ=255)
                                                                                                                                                                                                                      aCCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCCTTACCGTCGCTGCCTTCTGACCATTCCAGCTCATTTTTCCATTCTTTGCGCCACTGCCACCTTAAGCTTTTGATACCCAAATCAccgcctgttccgttggcgcgggt          >  2:378872/1‑208 (MQ=255)
                                                                                                                                                                                                                             aacgGAGAATACCCCCCGGTGGTGCCCCGGTGGCCCTTCCCGCGGCCCATTTGACGCTCCACCGCGGCGACCCGACGCTCCAGGCGGACGCTGGGACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAAc  <  2:561917/229‑1 (MQ=255)
                                                                                                                                                                                                                                 gagaAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAAc  <  1:133426/225‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
CGCTGCAACATCACCCAGCATCGCTTTAAAGCCGCCAATCTGGCGTGCAAATACGTTATAGATAACGACCGCAGGAATCGCTGCCACGAGGCCGATTGCCGTTGCTAACAGAGCTTCTGCGATACCCGGCGCAACGACTGCCAGGTTAGTGGTCTGC‑‑‑‑GTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAAC  >  NC_000913/3151341‑3151786

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