Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 777,208 A→C 10.2% K289N (AAA→AAC)  tolA → membrane anchored protein in TolA‑TolQ‑TolR complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913777,2080AC10.2% 602.8 / 17.4 216K289N (AAA→AAC) tolAmembrane anchored protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base A (131/56);  new base C (2/25);  total (133/83)
Fisher's exact test for biased strand distribution p-value = 3.45e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

GCAGAGAAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCAATAA  >  NC_000913/777059‑777357
                                                                                                                                                     |                                                                                                                                                     
 cagAGAAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCACAAg                                                                                                                                                      >  1:3557138/1‑150 (MQ=255)
 cagAGAAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAg                                                                                                                                                      >  1:1191379/1‑150 (MQ=255)
 cagAGAAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAg                                                                                                                                                      >  2:3383913/1‑150 (MQ=255)
 cagAGAAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAg                                                                                                                                                      >  1:3441914/1‑150 (MQ=255)
 cagAGAAAGATGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAATAAGCCGCTGCAGCAAAAg                                                                                                                                                      <  2:4309350/150‑1 (MQ=255)
  agagAAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGc                                                                                                                                                     >  2:1454089/1‑150 (MQ=255)
  agagAAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGc                                                                                                                                                     >  2:2898137/1‑150 (MQ=255)
   gagaAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCg                                                                                                                                                    >  1:740816/1‑150 (MQ=255)
      aaaGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGcc                                                                                                                                                 >  2:4466379/1‑150 (MQ=255)
         gCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCa                                                                                                                                              >  1:4414787/1‑150 (MQ=255)
         gCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCa                                                                                                                                              >  1:3016374/1‑150 (MQ=255)
           tGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCaga                                                                                                                                            >  2:1195615/1‑150 (MQ=255)
                ccGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAg                                                                                                                                       >  1:1930421/1‑150 (MQ=255)
                ccGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAg                                                                                                                                       >  2:3449056/1‑150 (MQ=255)
                 cGAAAAAAAAGCACAATAAAACGGGGCTGCTGAAAAGGCAGCCGCTTATAAGAACGCCGGGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCTGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAga                                                                                                                                      <  2:1016186/150‑1 (MQ=255)
                  gACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgat                                                                                                                                     >  2:2253579/1‑150 (MQ=255)
                     aaaaaaGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgatgat                                                                                                                                  >  2:2093561/1‑150 (MQ=255)
                     aaaaaaGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgatgat                                                                                                                                  >  1:3700373/1‑150 (MQ=255)
                        aaaGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATAtt                                                                                                                               <  2:3942442/150‑1 (MQ=255)
                             agcatAAAAAGCGGCTGCTGAAAAGGCTGCAGCTTATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCgg                                                                                                                          <  2:3330599/150‑1 (MQ=255)
                                           ctgctgAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTg                                                                                                            <  2:4356487/150‑1 (MQ=255)
                                           ctgctgAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTg                                                                                                            <  2:1917115/150‑1 (MQ=255)
                                               tgAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTaa                                                                                                        >  2:825892/1‑150 (MQ=255)
                                                gAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAg                                                                                                       >  1:2994944/1‑150 (MQ=255)
                                                gAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAg                                                                                                       >  1:568966/1‑150 (MQ=255)
                                                 aaaaGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGa                                                                                                      >  2:2186575/1‑150 (MQ=255)
                                                 aaaaGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGa                                                                                                      >  2:4195405/1‑150 (MQ=255)
                                                   aaGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAAt                                                                                                    >  1:2208172/1‑150 (MQ=255)
                                                    aGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATg                                                                                                   >  1:1783930/1‑150 (MQ=255)
                                                         gcagcTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACACAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCg                                                                                              >  2:3506430/1‑150 (MQ=255)
                                                           agcTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGaa                                                                                            >  2:4114311/1‑150 (MQ=255)
                                                                 tcagaccgcagcGGCAGAAAAAGCCGCCGCACACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGAAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACggg                                                                                      <  1:3103824/143‑1 (MQ=255)
                                                                 tAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACggg                                                                                      >  1:4205922/1‑150 (MQ=255)
                                                                  aaGAACGCAGCGGCAGAAAACGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAACGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACgggg                                                                                     <  2:545160/150‑1 (MQ=255)
                                                                  aaGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCCCACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACgggg                                                                                     >  2:4484176/1‑150 (MQ=255)
                                                                  aaGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACgggg                                                                                     >  1:3566867/1‑150 (MQ=255)
                                                                  aaGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACgggg                                                                                     >  2:4596474/1‑150 (MQ=255)
                                                                   aGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAGAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACggggg                                                                                    >  1:641176/1‑150 (MQ=255)
                                                                    gAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACgggggg                                                                                   >  1:4650378/1‑150 (MQ=255)
                                                                    gAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACgggggg                                                                                   >  1:4091604/1‑150 (MQ=255)
                                                                      atgcagcGTCTCAAAAAGCCGCCGCTGAAAATAAAGCGGCAGCGGAAAAAGCTGCAGCTTAAAAAGCCGCTTCAGAAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAg                                                                                 <  2:2865201/148‑1 (MQ=255)
                                                                             ggCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGaa                                                                          >  2:4663422/1‑150 (MQ=255)
                                                                                aGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAgg                                                                       >  1:4270381/1‑150 (MQ=255)
                                                                                       gccgccTCATAAAAAAAAGCGGCAGCGGAAAATGCTGCAGCTTAAAACGCCGGTGCAGCAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAAt                                                                <  2:4050905/150‑1 (MQ=255)
                                                                                        ccgccgCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATg                                                               >  2:2649596/1‑150 (MQ=255)
                                                                                         cgccgcAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGc                                                              >  1:3712548/1‑150 (MQ=255)
                                                                                         cgccgcAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGc                                                              <  2:4780607/150‑1 (MQ=255)
                                                                                          gccgcAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCt                                                             <  2:3217445/150‑1 (MQ=255)
                                                                                           ccgcATACAAAAACGCGGACGGGGCAATCGGTTCCGCTGAAAACGCCGCTGCAGCAACTGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGATCAATGCtt                                                            <  2:2839363/150‑1 (MQ=255)
                                                                                           ccgcAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCtt                                                            <  2:2993531/150‑1 (MQ=255)
                                                                                             gcataaaaaaatGCGGCAGCGGAAAAAGCTGCAGCTTAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCg                                                          <  2:3364325/144‑1 (MQ=255)
                                                                                             gcAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCg                                                          >  1:3969902/1‑150 (MQ=255)
                                                                                              cAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGc                                                         <  2:4728737/150‑1 (MQ=255)
                                                                                                gACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCt                                                       >  1:3338193/1‑150 (MQ=255)
                                                                                                   aaaaaaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGACCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGcc                                                    <  2:3566430/150‑1 (MQ=255)
                                                                                                       aaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCCAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGGACGAATGCACCGAAAACGGGGGGAGGGGCGACCGGGAACAATGCTTCGCCTGCCgggc                                                >  2:3824472/1‑149 (MQ=255)
                                                                                                       aaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGa                                                <  2:3948989/150‑1 (MQ=255)
                                                                                                       aaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGa                                                >  2:2669961/1‑150 (MQ=255)
                                                                                                       aaGCGGCAGCGGAAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGa                                                <  2:161683/150‑1 (MQ=255)
                                                                                                         gcggcagcggcaAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGt                                              >  1:4479810/1‑150 (MQ=255)
                                                                                                          cggcagcggcaAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTg                                             <  2:888571/150‑1 (MQ=255)
                                                                                                           ggctgcggaaAACGCTGCCGCTTAAAACGCCGCTGCACCAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCTAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTgg                                            <  1:3203638/150‑1 (MQ=255)
                                                                                                            gcagcggcaAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGt                                           >  1:3037430/1‑150 (MQ=255)
                                                                                                             cagcggcaAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTa                                          <  2:4797051/150‑1 (MQ=255)
                                                                                                             cagcggcaAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTa                                          >  2:135611/1‑150 (MQ=255)
                                                                                                             cagcggcaAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGAGGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGtc                                          >  2:4497254/1‑149 (MQ=255)
                                                                                                              agcggaaAACGCGGCAGCTTAAAAACCCGCTGCAGCAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTaa                                         <  2:23131/150‑1 (MQ=255)
                                                                                                               gcggcaAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAAt                                        >  1:792420/1‑150 (MQ=255)
                                                                                                                 ggcaAAAGCTCCTGCTGAAAACGCCGCTGCACCAAACGCGGCCGCAGCGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATAc                                      <  1:4661798/150‑1 (MQ=255)
                                                                                                                 ggcaAAAGCTCCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATAc                                      <  2:949652/150‑1 (MQ=255)
                                                                                                                 ggaaaacGCTGCAGCTGAAAACGCCGCGGCAGCAAACGCGGCCGCAGAGGCAGATGATCTTTTCGGTGAGCTACGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATAc                                      <  2:78337/147‑1 (MQ=255)
                                                                                                                  ccaAACGCTGCCTCTCAAACACCGGCGGCACCAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCTAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACt                                     <  1:90717/149‑1 (MQ=255)
                                                                                                                    aaaaGCTGCAGCTGAAAACGCCGCGGCAGAAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTaa                                   <  2:546514/150‑1 (MQ=255)
                                                                                                                    aaaaGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTaa                                   >  1:1122547/1‑150 (MQ=255)
                                                                                                                     aaaGCTGCAGCTTAAAAAGCCGCTGCCGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTaaa                                  <  2:626670/150‑1 (MQ=255)
                                                                                                                     aaaGCTGCAGCTGAAACAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTaaa                                  <  2:2508765/150‑1 (MQ=255)
                                                                                                                     aaaGCTCCAGCTTAAAAAGCCGCTGAACAAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAACAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTaaa                                  <  2:3706166/150‑1 (MQ=255)
                                                                                                                             aGCTTAAAAAGCCGCTTCAGAAAACGCGGCCGCAGAGGCAGATGATTTTTTCGGTGAGATAAGCTCTGTTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATgg                          <  1:3281751/150‑1 (MQ=255)
                                                                                                                              gCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGgc                         >  1:2369320/1‑150 (MQ=255)
                                                                                                                                  aaaaaGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCAt                     >  2:1287461/1‑150 (MQ=255)
                                                                                                                                     aaGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAg                  >  2:2243871/1‑150 (MQ=255)
                                                                                                                                     aaGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAg                  >  2:1974542/1‑150 (MQ=255)
                                                                                                                                       gCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAggg                >  1:1668749/1‑150 (MQ=255)
                                                                                                                                        ccGCTGCAGCAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAgggg               <  1:3845198/150‑1 (MQ=255)
                                                                                                                                        ccGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAgggg               >  2:1419851/1‑150 (MQ=255)
                                                                                                                                        ccGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAgggg               >  2:1146344/1‑150 (MQ=255)
                                                                                                                                        ccGCTGCACCAATAGCGGCCGCAGAGGCAGATGTTTTTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAgggg               <  1:3654744/150‑1 (MQ=255)
                                                                                                                                         cGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGc              >  2:3647865/1‑150 (MQ=255)
                                                                                                                                            tgctgcAAATGCGGGCGCTGTGGCTGTTGATATGTTGGGTGAGGTTCGCTCTGGTAAGAATGCACCGAATAGGGGGGGAGGGGGGCACGGGATCACTGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACTATGGCGCATCCGGGGCCGa           <  2:3834173/150‑1 (MQ=255)
                                                                                                                                               agcaAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGatat        >  2:1706639/1‑150 (MQ=255)
                                                                                                                                               agcaAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGatat        <  1:1426302/150‑1 (MQ=255)
                                                                                                                                               agaaaacGCGGCCGCAGAGGCATATTATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGatat        <  1:1665105/147‑1 (MQ=255)
                                                                                                                                               agaaaaaGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAATCGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGatat        <  2:545007/147‑1 (MQ=255)
                                                                                                                                                gcaAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATc       >  1:2780622/1‑150 (MQ=255)
                                                                                                                                                 caAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCa      >  1:303656/1‑150 (MQ=255)
                                                                                                                                                 caAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCa      <  1:2868431/150‑1 (MQ=255)
                                                                                                                                                 caAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCa      >  1:303721/1‑150 (MQ=255)
                                                                                                                                                 caAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCa      >  1:3633894/1‑150 (MQ=255)
                                                                                                                                                  aaaaGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCaa     >  1:3832717/1‑150 (MQ=255)
                                                                                                                                                     aGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCAATaa  <  1:4752394/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAGAGAAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCAATAA  >  NC_000913/777059‑777357

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑