Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,018,983 A→T 6.8% intergenic (+46/+30) matP → / ← ompA Ter macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,018,9830AT 6.8% 658.6 / 30.5 192intergenic (+46/+30)matP/ompATer macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base A (100/79);  new base T (5/8);  total (105/87)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-03

Only 100 of 200 total aligned reads displayed.

ATGCGGAACACCAGCATACGCGCAAAAAATCGATCGATCTGGAATTTATTGTCTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTA  >  NC_000913/1018732‑1019249
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                          
                             tcgatcgaTCTGGAATTTATTGTCTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACa                                                                                                                                                                                                                         <  2:237355/274‑1 (MQ=255)
                                               aTTGTCTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGttttt                                                                                                                                                                                                                                                              >  2:334554/1‑219 (MQ=255)
                                               aTTGTCTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTa                                                                                                                                                                                                       >  2:161222/1‑274 (MQ=255)
                                                                             cAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCa                                                                                                                                                                        <  2:53031/275‑1 (MQ=255)
                                                                                cttagcgtgaaaaCGCTTATGGAAAGGATAGGTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTATTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAgg                                                                                                                                                                      <  4:260144/265‑1 (MQ=255)
                                                                                       gTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGAt                                                                                                                                                              <  4:67564/275‑1 (MQ=255)
                                                                                                  tCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGc                                                                                                                                                                                                                                     >  4:290207/1‑193 (MQ=255)
                                                                                                    tGAAACGATTGTTGTGCTGATTGAAGAGGCGGAAAAAAAAGAGAAATACGCGAAAAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACg                                                                                                                                                  <  4:114905/274‑1 (MQ=255)
                                                                                                          gATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCgcagca                                                                                                                                                                                                                                                                         >  5:524317/1‑149 (MQ=255)
                                                                                                                  aGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGGtt                                                                                                                                                                                                                                                                 <  6:602710/149‑1 (MQ=255)
                                                                                                                  aGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGGtt                                                                                                                                                                                                                                                                 <  6:872243/149‑1 (MQ=255)
                                                                                                                  aGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGGtt                                                                                                                                                                                                                                                                 <  5:896311/149‑1 (MQ=255)
                                                                                                                  aGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGtt                                                                                                                                                                                                                                                                 <  6:1186674/149‑1 (MQ=255)
                                                                                                                  aGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGtt                                                                                                                                                                                                                                                                 >  7:953180/1‑149 (MQ=255)
                                                                                                                    cTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGtttt                                                                                                                                                                                                                                                               >  6:1427518/1‑149 (MQ=255)
                                                                                                                    cTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgcag                                                                                                                                                      >  2:29434/1‑254 (MQ=255)
                                                                                                                     tGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAGTTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGGttttt                                                                                                                                                                                                                                                              <  7:801073/149‑1 (MQ=255)
                                                                                                                      gATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGGtttttt                                                                                                                                                                                                                                                             <  5:349988/149‑2 (MQ=255)
                                                                                                                       aTTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCt                                                                                                                                                                                                                                                            >  6:154926/1‑149 (MQ=255)
                                                                                                                            aGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAg                                                                                                                                                                                                                                                       >  8:14522/1‑149 (MQ=255)
                                                                                                                             gATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGa                                                                                                                                                                                                                                                      >  5:1393298/1‑149 (MQ=255)
                                                                                                                                  ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTATTTCTACCAGACgaga                                                                                                                                                                                                                                                 >  5:1348482/1‑149 (MQ=255)
                                                                                                                                   gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGGAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                                                                                                                                                <  5:303827/149‑1 (MQ=255)
                                                                                                                                   gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                                                                                                                                                <  6:1348482/149‑1 (MQ=255)
                                                                                                                                   gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                                                                                                                                                >  7:538181/1‑149 (MQ=255)
                                                                                                                                       aCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTa                                                                                                                                                                                                                                            >  7:1138934/1‑149 (MQ=255)
                                                                                                                                          aaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCgg                                                                                                                                                                                                                                                                     >  5:815385/1‑121 (MQ=255)
                                                                                                                                          aaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                                                                                                                                                         <  6:636277/149‑1 (MQ=255)
                                                                                                                                           agagAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGcc                                                                                                                                                                                                                                        >  8:939361/1‑149 (MQ=255)
                                                                                                                                            gagaAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGttt                                                                                                                                                                                                                                                                <  7:588917/124‑1 (MQ=255)
                                                                                                                                                  tACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGTTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCt                                                                                                                                                                                                                                 <  5:154926/149‑1 (MQ=255)
                                                                                                                                                   aCGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTg                                                                                                                                                                                                                                >  7:5272/1‑149 (MQ=255)
                                                                                                                                                    cgcgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGa                                                                                                                                                                                                                               <  8:864414/149‑1 (MQ=255)
                                                                                                                                                       gAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAGCTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGtt                                                                                                                                                                                                                            >  6:106664/1‑149 (MQ=255)
                                                                                                                                                       gAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGtt                                                                                                                                                                                                                            >  8:809508/1‑149 (MQ=255)
                                                                                                                                                         aTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGGttttt                                                                                                                                                                                                                                                              <  7:852359/113‑1 (MQ=255)
                                                                                                                                                         aTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGGttttt                                                                                                                                                                                                                                                              >  8:852359/1‑113 (MQ=255)
                                                                                                                                                             aaTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCggg                                                                                                                                                                                                                                                                    >  5:208873/1‑103 (MQ=255)
                                                                                                                                                                   ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTc                                                                                                                                                                                                                                                             <  7:351496/104‑1 (MQ=255)
                                                                                                                                                                    tctcTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGa                                                                                                                                                                                                               <  7:456817/149‑1 (MQ=255)
                                                                                                                                                                     ctctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTAACGTTGTCACAGGTTTGGCCAGTAACCGGGTTGGATTCGCCCATACCAGGTCCGGAGAt                                                                                >  4:237813/1‑275 (MQ=255)
                                                                                                                                                                      tctGAAGCCGGATCTGCAGGCATTGCGGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAACGGCCAAAAAAACCCCGCAGCAGCGGGGAGTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATa                                                                                                                                                                                                             >  7:1287572/1‑149 (MQ=255)
                                                                                                                                                                      tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCt                                                                                                                                                                                                                                       >  6:386123/1‑123 (MQ=255)
                                                                                                                                                                      tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATa                                                                                                                                                                                                             >  6:1391586/1‑149 (MQ=255)
                                                                                                                                                                             cAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTaa                                                                                                                                                                                                      <  7:1368980/149‑1 (MQ=255)
                                                                                                                                                                               ggATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACt                                                                                                                                                                                                    >  8:785427/1‑149 (MQ=255)
                                                                                                                                                                               ggATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGcc                                                                      >  2:255798/1‑275 (MQ=255)
                                                                                                                                                                                gATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACtt                                                                                                                                                                                                   >  6:708534/1‑149 (MQ=255)
                                                                                                                                                                                gATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcgacga                                                                                                                                                                                  >  4:322485/1‑166 (MQ=255)
                                                                                                                                                                                    tGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCg                                                                                                                                                                              <  2:322485/166‑1 (MQ=255)
                                                                                                                                                                                        ggCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgcag                                                                                                                                                      >  2:2281/1‑186 (MQ=255)
                                                                                                                                                                                         gCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCgg                                                                                                                                                                                                                                                                     <  8:5301/74‑1 (MQ=255)
                                                                                                                                                                                            ttGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACg                                                                                                                                                  <  4:2281/186‑1 (MQ=255)
                                                                                                                                                                                                   gTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATc                                                                                                                                                                                <  7:785427/149‑1 (MQ=255)
                                                                                                                                                                                                     aaGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCg                                                                                                                                                                              <  8:813256/149‑1 (MQ=255)
                                                                                                                                                                                                       ggAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGa                                                                                                                                                                            <  6:524317/149‑1 (MQ=255)
                                                                                                                                                                                                       ggAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTaa                                              >  2:76422/1‑275 (MQ=255)
                                                                                                                                                                                                        gAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTaa                                                                                                                                                                                                      <  6:1250030/122‑1 (MQ=255)
                                                                                                                                                                                                        gAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAg                                                                                                                                                                           >  8:577290/1‑149 (MQ=255)
                                                                                                                                                                                                         aaTAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGc                                                                                                                                                                          >  6:41831/1‑149 (MQ=255)
                                                                                                                                                                                                            aaCTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCCCAGGTGTTGCCAGTAACCGGGtt                                                                                                           >  4:158599/1‑210 (MQ=255)
                                                                                                                                                                                                               tGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCa                                                                                                                                                                    >  8:711042/1‑149 (MQ=255)
                                                                                                                                                                                                                 aCGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGt                                                                                                                                                                  >  6:117279/1‑149 (MQ=255)
                                                                                                                                                                                                                   gAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATCCCACGTGCGGAGATCTTGTCTGcc                                                                      >  4:47935/1‑239 (MQ=255)
                                                                                                                                                                                                                   gAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGACAACTTAAGCATGCGGCTGAGTTACAACGTCTTTGATCCCTTTAACTTCGATCTCTACGCGACGATCCGGTGCCAGGCAGTcc                                                                                                                                                                >  8:12647/1‑148 (MQ=255)
                                                                                                                                                                                                                       gTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCa                                                                                                                                                            <  6:1262530/149‑1 (MQ=255)
                                                                                                                                                                                                                       gTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCa                                                                                                                                                            >  7:1438103/1‑149 (MQ=255)
                                                                                                                                                                                                                            ttCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                                                                                                                                                       >  6:996773/1‑149 (MQ=255)
                                                                                                                                                                                                                            ttCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGAGGCTGa                                                                                                                                                                                                                               >  6:536288/1‑77 (MQ=255)
                                                                                                                                                                                                                            ttCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGAGGCTGa                                                                                                                                                                                                                               <  5:536288/77‑1 (MQ=255)
                                                                                                                                                                                                                               aaTTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGt                                                                                                                                                          <  6:815740/143‑1 (MQ=255)
                                                                                                                                                                                                                               aaTTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgcagca                                                                                                                                                    >  8:495389/1‑149 (MQ=255)
                                                                                                                                                                                                                                aTTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCAc                                                                                                                                                   >  6:260995/1‑149 (MQ=255)
                                                                                                                                                                                                                                    aCTAAAGGCAAAAAAAGCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTAAAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGa                 <  2:128394/275‑1 (MQ=255)
                                                                                                                                                                                                                                      tAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGtt                                                                                                                                             >  7:734777/1‑149 (MQ=255)
                                                                                                                                                                                                                                      tAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGtt                                                                                                                                             >  5:1093730/1‑149 (MQ=255)
                                                                                                                                                                                                                                       aaaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTc                                                                                                                                                                 >  7:877887/1‑128 (MQ=255)
                                                                                                                                                                                                                                       aaaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGttt                                                                                                                                            >  8:767401/1‑149 (MQ=255)
                                                                                                                                                                                                                                        aaagCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTc                                                                                                                                           <  5:663725/146‑1 (MQ=255)
                                                                                                                                                                                                                                        aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTc                                                                                                                                           <  7:578193/149‑1 (MQ=255)
                                                                                                                                                                                                                                         aGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTa                                                                                                                                                                                          >  6:1435377/1‑101 (MQ=255)
                                                                                                                                                                                                                                          ggCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAg                                                                                                                                                                       >  7:254665/1‑119 (MQ=255)
                                                                                                                                                                                                                                          ggCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAg                                                                                                                                                                       <  8:254665/119‑1 (MQ=255)
                                                                                                                                                                                                                                          ggCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCAc                                                                                                                                         >  5:683121/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aaaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAGCTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGtt                                                                                                                                      >  7:1106812/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aaaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGtt                                                                                                                                      >  5:1072741/1‑149 (MQ=255)
                                                                                                                                                                                                                                                aaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt                                                                                                                                    <  7:1232366/149‑1 (MQ=255)
                                                                                                                                                                                                                                                aaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGGTCTTGTCTGcc                                                                      >  4:218884/1‑211 (MQ=255)
                                                                                                                                                                                                                                                aaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAACTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTc                                                                                                                                           >  1:25597/1‑142 (MQ=255)
                                                                                                                                                                                                                                                   aaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcac                                                                                                                                 >  5:420143/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   aaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcac                                                                                                                                 <  5:1237714/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcac                                                                                                                                 >  8:685106/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    aaCNCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGt                                                                                                                                       <  3:25597/142‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca                                                                                                                                >  5:581927/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    aaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca                                                                                                                                >  8:1128093/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    aaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGa                                                                  <  2:218884/211‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTa  <  4:76422/275‑1 (MQ=255)
                                                                                                                                                                                                                                                     aCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCAc                                                                                                                                         <  5:1261369/139‑1 (MQ=255)
                                                                                                                                                                                                                                                       cccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGgt                                                                                                                             >  8:392750/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        ccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTGTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGgtg                                                                                                                            >  7:842790/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                          
ATGCGGAACACCAGCATACGCGCAAAAAATCGATCGATCTGGAATTTATTGTCTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGCACGGCGCTCGGACAGACCCTGGTTGTA  >  NC_000913/1018732‑1019249

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