Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,388,285 A→G 5.2% intergenic (+24/+20) tyrR → / ← tpx aromatic amino acid biosynthesis and transport regulon transcriptional regulator; autorepressor; ATPase; phosphatase/lipid hydroperoxide peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,388,2850AG 5.2% 524.3 / 11.9 155intergenic (+24/+20)tyrR/tpxaromatic amino acid biosynthesis and transport regulon transcriptional regulator; autorepressor; ATPase; phosphatase/lipid hydroperoxide peroxidase
Reads supporting (aligned to +/- strand):  ref base A (75/72);  new base G (5/3);  total (80/75)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 166 total aligned reads displayed.

TTATGAGCTGCGTCCACAGGATATTTTGTTGCCGGATTATGACGCCGCAACGGTAGCCGTGGGCGAAGATGCGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTAAGCTTGCAGAAATTCAGCGTTACGGAAAGTGGAGAGGGTGATAACGTTGTTCAGACCTTCTGCGCCGCAGAAACGAGACTGGGCGAACGGCAGATCGG  >  NC_000913/1388015‑1388558
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                 
ttATGAGCTGCGTCCACAGGATATTTTGTTGCCGGATTATGACGCCGCAACGGTAGCCGTGGGCGAAGATGCGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAAGATGGTCTGAGTCAGAAGAAGAACGAAAGGTAAGCGCGGATATGCCTGATTGTGCAAcccc                                                                                                                                                                                                                                                                               >  4:94029/1‑275 (MQ=255)
     aGCTGCGTCCCAAGGATATTTTGTGGCCGGATTATGACGCCGCAACGGTAGCCGTGGGCGAAAATGCGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAg                                                                                                                                                                                                                                                                          <  4:287835/275‑1 (MQ=255)
                    atatTTTGTTGCCGGATTATGACGCCGCAACGGTAGCCGTGGGCGAAGATGCGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAAACAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATGTGCCTGATGGTGCAACACCATCAGGcctgtctgtgata                                                                                                                                                                                                                                                            >  2:309433/1‑262 (MQ=255)
                                                        gcGTGGGCGAATTTGCTATGGGAGGTTCCCTGAACGAACACCCCAACCTTTTTGAACGCTCGGTATTAACCCATCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGTATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCt                                                                                                                                                                                                                       <  4:321808/274‑1 (MQ=255)
                                                            gggCGAAGATGCGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTACCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATNCCTGATGGTGCAACACCAACAGGCAAATTAAATTATGCTTTAAGTACAGCAAGAGATGCTTCGAAATCAGGCTCgg                                                                                                                                                                                                                    >  3:67427/1‑274 (MQ=255)
                                                                                           gAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCaa                                                                                                                                                                                                                                                                                   <  4:271112/180‑1 (MQ=255)
                                                                                           gAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCa                                                                                                                                                                                                                                                                              >  2:116692/1‑185 (MQ=255)
                                                                                                          tttgaagagtCGGTATTAAACAAGCGTTATCGAAATAATCCAAGACAGAGCAAACGGGCAAAACGTCGCGGCGTTTAACATACCGCGATTGCCAATAAGTTGCGGGAAGATGGTCTGAGTCAGAAGAAGAAAGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAtttt                                                                                                                                                                       <  2:195296/264‑1 (MQ=255)
                                                                                                                   tCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAATAGTAAGCGCGAATATGCCTGCTGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCACTAAcc                                                                                                                                                            >  2:244887/1‑275 (MQ=255)
                                                                                                                           aaCCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAac                                                                                                                                                                                                                                                                                  <  7:1326407/149‑1 (MQ=255)
                                                                                                                              ccAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACAcc                                                                                                                                                                                                                                                                               >  7:912972/1‑149 (MQ=255)
                                                                                                                              ccAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACAcc                                                                                                                                                                                                                                                                               >  7:307122/1‑149 (MQ=255)
                                                                                                                              ccAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACAcc                                                                                                                                                                                                                                                                               >  6:723733/1‑149 (MQ=255)
                                                                                                                               cAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCa                                                                                                                                                                                                                                                                              >  7:1256078/1‑149 (MQ=255)
                                                                                                                                 gCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCa                                                                                                                                                                                                                                                                              >  4:108223/1‑147 (MQ=255)
                                                                                                                                  cTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCa                                                                                                                                                                                                                                                                           <  8:513562/149‑1 (MQ=255)
                                                                                                                                   tttATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAg                                                                                                                                                                                                                                                                          <  6:164686/149‑1 (MQ=255)
                                                                                                                                   tttATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAg                                                                                                                                                                                                                                                                          <  7:1433875/149‑1 (MQ=255)
                                                                                                                                   tttATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAg                                                                                                                                                                                                                                                                          <  5:1460010/149‑1 (MQ=255)
                                                                                                                                   tttATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGCCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAAAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAg                                                                                                                                                                                                                                                                          <  5:198187/149‑1 (MQ=255)
                                                                                                                                     tATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAg                                                                                                                                                                                                                                                                          <  2:108223/147‑1 (MQ=255)
                                                                                                                                       tCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCat                                                                                                                                                                                                                                                                      >  6:1385791/1‑149 (MQ=255)
                                                                                                                                        cGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCa                                                                                                                                                                                                                                                                           <  6:236840/143‑1 (MQ=255)
                                                                                                                                           aaTTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTa                                                                                                                                                                                                                                                                  >  6:438166/1‑149 (MQ=255)
                                                                                                                                             ttATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTaaa                                                                                                                                                                                                                                                                <  6:913484/149‑1 (MQ=255)
                                                                                                                                             ttATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTaaa                                                                                                                                                                                                                                                                <  8:1148166/149‑1 (MQ=255)
                                                                                                                                             ttATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTaaa                                                                                                                                                                                                                                                                <  5:77730/149‑1 (MQ=255)
                                                                                                                                             ttATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTaaa                                                                                                                                                                                                                                                                <  7:597071/149‑1 (MQ=255)
                                                                                                                                               aTCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAAtt                                                                                                                                                                                                                                                              >  5:301231/1‑149 (MQ=255)
                                                                                                                                                 cccAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGCATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTAt                                                                                                                                                                                                                                                            >  7:208539/1‑149 (MQ=255)
                                                                                                                                                   cAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCaa                                                                                                                                                                                                                                                                                   <  2:341701/124‑1 (MQ=255)
                                                                                                                                                    aGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCt                                                                                                                                                                                                                                                         <  6:270338/149‑1 (MQ=255)
                                                                                                                                                       aCGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                                                                                                                                                                                                      <  5:598961/149‑1 (MQ=255)
                                                                                                                                                        cgcgCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCa                                                                                                                                                                                                                                                     <  5:1007912/149‑1 (MQ=255)
                                                                                                                                                         gcgcAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTc                                                                                                                                                                                                                      <  2:367367/179‑1 (MQ=255)
                                                                                                                                                           gcAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTa                                                                                                                                                                                                                                                  <  5:738670/149‑1 (MQ=255)
                                                                                                                                                                   gCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCaga                                                                                                                                                                                                                                          >  7:137515/1‑149 (MQ=255)
                                                                                                                                                                   gCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCaga                                                                                                                                                                                                                                          <  8:964941/149‑1 (MQ=255)
                                                                                                                                                                     aaaaCGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgc                                                                                                                                                                                                                                        <  8:1256078/149‑1 (MQ=255)
                                                                                                                                                                        aCGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgctgc                                                                                                                                                                                                                                     >  8:294338/1‑149 (MQ=255)
                                                                                                                                                                           tctcGGCGTTTCACATACCGCGATTGGCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCTCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTc                                                                                                                                                                                                                                  <  8:1323200/149‑1 (MQ=255)
                                                                                                                                                                           tctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTc                                                                                                                                                                                                                                  <  7:836869/149‑1 (MQ=255)
                                                                                                                                                                           tctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTc                                                                                                                                                                                                                                  <  5:111586/149‑1 (MQ=255)
                                                                                                                                                                             tcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGt                                                                                                                                                                                                                                >  5:6307/1‑149 (MQ=255)
                                                                                                                                                                              cGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTa                                                                                                                                                                                                                               <  8:1015137/149‑1 (MQ=255)
                                                                                                                                                                               ggCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTaa                                                                                                                                                                                                                              <  5:458480/149‑1 (MQ=255)
                                                                                                                                                                                  gTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATcc                                                                                                                                                                                                                           <  7:917248/149‑1 (MQ=255)
                                                                                                                                                                                            ccGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtgg                                                                                                                                                                                                                 >  5:1095174/1‑149 (MQ=255)
                                                                                                                                                                                              gcgATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGgtggtg                                                                                                                                                                                                               >  8:792569/1‑149 (MQ=255)
                                                                                                                                                                                                gATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGAt                                                                                                                                                                                                             <  7:94969/149‑1 (MQ=255)
                                                                                                                                                                                                 aTTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGAtt                                                                                                                                                                                                            <  6:311022/149‑1 (MQ=255)
                                                                                                                                                                                                     ccAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCAt                                                                                                                                                                                                        <  6:632090/149‑1 (MQ=255)
                                                                                                                                                                                                      cAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATc                                                                                                                                                                                                       <  7:888149/149‑1 (MQ=255)
                                                                                                                                                                                                       aaTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtgg                                                                                                                                                                                                                 >  7:397666/1‑138 (MQ=255)
                                                                                                                                                                                                       aaTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCCGCCAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTAAGCTTGCAGAAATTCAGCGTTACGGaaa                                                                        <  2:94029/275‑1 (MQ=255)
                                                                                                                                                                                                                 gggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc                                                                                                                                                                                            <  7:792569/149‑1 (MQ=255)
                                                                                                                                                                                                                   gAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAGGc                                                                                                                                                                                                                                                                        <  7:75665/71‑1 (MQ=25)
                                                                                                                                                                                                                   gAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAGGc                                                                                                                                                                                                                                                                        >  8:75665/1‑71 (MQ=25)
                                                                                                                                                                                                                     atatGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaa                                                                                                                                                                                        >  5:81414/1‑149 (MQ=255)
                                                                                                                                                                                                                       atGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaaaa                                                                                                                                                                                      <  8:698081/149‑1 (MQ=255)
                                                                                                                                                                                                                        tGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAAt                                                                                                                                                                                     <  7:539165/149‑1 (MQ=255)
                                                                                                                                                                                                                         ggTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATc                                                                                                                                                                                    >  6:997678/1‑149 (MQ=255)
                                                                                                                                                                                                                          gTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATca                                                                                                                                                                                   >  5:66117/1‑149 (MQ=255)
                                                                                                                                                                                                                           tCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATcac                                                                                                                                                                                  >  8:769906/1‑149 (MQ=255)
                                                                                                                                                                                                                              gAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGATCGGTGGTGATTTCATCCACTAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGTCCATCAGCAGTTGCCACACCGTAATCTTGCAGAAATTCAGCGTTACGGAATGTGGAGAGGGTGATAACgttgt                                                  >  2:286407/1‑274 (MQ=255)
                                                                                                                                                                                                                                gTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGt                                                                                                                                                                             <  6:748352/149‑1 (MQ=255)
                                                                                                                                                                                                                                 tCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTc                                                                                                                                                                            >  5:99595/1‑149 (MQ=255)
                                                                                                                                                                                                                                 tCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTc                                                                                                                                                                            >  7:1032272/1‑149 (MQ=255)
                                                                                                                                                                                                                                 tCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTc                                                                                                                                                                            >  7:1005436/1‑149 (MQ=255)
                                                                                                                                                                                                                                     aagaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCAAATGCCACACCGTAAGCTTGCACAAATTCAGCGTTACGGAAAGTGGGAAGGGTGATACAGTTGTTCaagact                                          >  2:279013/1‑270 (MQ=255)
                                                                                                                                                                                                                                       gaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTc                                                                                                                                                                      <  7:1146581/149‑1 (MQ=255)
                                                                                                                                                                                                                                         agaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGt                                                                                                                                                                    <  5:1192725/149‑1 (MQ=255)
                                                                                                                                                                                                                                          gaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTc                                                                                                                                                                   >  6:90026/1‑149 (MQ=255)
                                                                                                                                                                                                                                          gaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTc                                                                                                                                                                   <  6:600161/149‑1 (MQ=255)
                                                                                                                                                                                                                                            aCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCaa                                                                                                                                                                 >  8:122913/1‑149 (MQ=255)
                                                                                                                                                                                                                                              gAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATa                                                                                                                                                               >  8:1059977/1‑149 (MQ=255)
                                                                                                                                                                                                                                               aaGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATaa                                                                                                                                                              >  6:49826/1‑149 (MQ=255)
                                                                                                                                                                                                                                               aaGAGGAAGGGCGAATAGGCCTGTGGGGGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAGTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTAAGCTTGCAGAAATTCAGCGTTACGGAAAGTGGAGAGGGTGATAACGTTGTTCAGACCTTCTGcgccgc                                <  4:128466/275‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgctgct                                                                                                                                                                                                                                    <  7:688375/73‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgctgct                                                                                                                                                                                                                                    >  8:688375/1‑73 (MQ=255)
                                                                                                                                                                                                                                                         gcgAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGc                                                                                                                                                    >  8:425642/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                tGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCt                                                                                                                                             <  6:570849/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGc                                                                                                                                           <  8:208539/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATa                                                                                                                                                               >  7:1206821/1‑128 (MQ=255)
                                                                                                                                                                                                                                                                   cTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATa                                                                                                                                                               <  8:1206821/128‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCa                                                                                                                                         >  5:990073/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    tGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCa                                                                                                                                         >  6:1222138/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCa                                                                                                                                         >  7:209202/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCa                                                                                                                                         <  6:1129103/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCa                                                                                                                                         >  8:534136/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCa                                                                                                                                         <  7:425642/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCt                                                                                                                                    >  8:251815/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         gTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCt                                                                                                                                    >  8:610838/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         gTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCt                                                                                                                                    <  6:81414/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggt                                                                                                                                                                                                                   >  5:912964/1‑68 (MQ=255)
                                                                                                                                                                                                                                                                           gCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCttt                                                                                                                                  >  6:1234236/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            cAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTc                                                                                                                                 >  8:1350589/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                             aaCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCa                                                                                                                                <  5:1385791/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                              acacCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCAcg                                                                                                                                                 <  4:375250/131‑1 (MQ=255)
                                                                                                                                                                                                                                                                              acacCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTAAGCTTGCAGAAATTCAGCGTTACGGAAAGTGGAGAGGGAGATAACGTTGTTCAGACCTTCTGCGCCGCAGAAACGAGACTGGGCGAACGGCAGATCgt  >  1:14615/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                 
TTATGAGCTGCGTCCACAGGATATTTTGTTGCCGGATTATGACGCCGCAACGGTAGCCGTGGGCGAAGATGCGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTAAGCTTGCAGAAATTCAGCGTTACGGAAAGTGGAGAGGGTGATAACGTTGTTCAGACCTTCTGCGCCGCAGAAACGAGACTGGGCGAACGGCAGATCGG  >  NC_000913/1388015‑1388558

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