Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,314,233 A→T 6.1% intergenic (+36/+111) yjdP → / ← phnP putative periplasmic protein/5‑phospho‑alpha‑D‑ribosyl 1,2‑cyclic phosphate phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,314,2330AT 6.1% 467.3 / 16.3 132intergenic (+36/+111)yjdP/phnPputative periplasmic protein/5‑phospho‑alpha‑D‑ribosyl 1,2‑cyclic phosphate phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (57/67);  new base T (4/4);  total (61/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-02

Only 100 of 137 total aligned reads displayed.

CAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTG  >  NC_000913/4313967‑4314505
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                
  ggCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggag                                                                                                                                                                                                                                                                           <  2:119440/272‑1 (MQ=255)
               gCTTATCACCCTGACCGGTGGCGGGGCTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAgcg                                                                                                                                                                                                                                                           <  4:361952/275‑1 (MQ=255)
               gCTTATACCCCGTACCGTGGGCGGAACTATGAAGATTCGCGCGACGACGCCTGGCAGCGGGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAgc                                                                                                                                                                                                                                                            <  4:243070/274‑1 (MQ=255)
               gCTTATAACCCTGACCGGGGGGGAGACTATGAGGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAgcg                                                                                                                                                                                                                                                           <  4:135739/275‑1 (MQ=255)
                  tttagcccTGGGCGTGGGGGGGACGATGCAGATTCCGGGGACGACGGCTGGCGGCACGAAGTTAGCGACGATCGGCGCAGACAATATAACGGCCGCCGCCGCCGGTTTGAGGATGGCCGCGGGCAGCCGGACGATGGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTGGGGTATGGAGCAAGCAACGCCCCTCACCCCAATGGGGGGAGGGAGAAAACGAGCGCa                                                                                                                                                                                                                                                         <  4:359249/269‑1 (MQ=255)
                  tataTCCCGGGCCGTGGGCGGAACGATGAAGAGTCGGGGGAAGACGGGTGGCAGCGTGAAGGGAGTGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTGGAAGATGGCCGGAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCa                                                                                                                                                                                                                                                         <  4:410512/274‑1 (MQ=255)
                  tataACCCGGACCGTGGGCGAGACTATGAAGATTCCGGCGGCGACGGCTGGCAGTGTGAAGTGAGCGACGTTCGGCGCAGACAATATGACGACGGCCGCCGCCAGGTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCAGGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCa                                                                                                                                                                                                                                                         <  4:319496/274‑1 (MQ=255)
                  tagaaCCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGAGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGGGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCa                                                                                                                                                                                                                                                         <  2:176596/271‑1 (MQ=255)
                      aCCCTGACCGTAGGCGAAACTATGAAGATTCGGGCGACGACGGCTGGCAGCGGGAAATGAGCGACGATCTGCGCAGACAATAGGGCGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGAGGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATAtt                                                                                                                                                                                                                                                    <  2:337721/275‑1 (MQ=255)
                               gTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCaa                                                                                                                                                                                                                                                                                  <  2:347677/236‑1 (MQ=255)
                                        aCTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTAAGGAGCAGGCATCGCCCCTCACCCCAATGGGGGGAGGGAGAAAAACACCGCAATCTTCAAAATCGAGCAAAATc                                                                                                                                                                                                                                   >  4:312404/1‑274 (MQ=255)
                                                   tCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGGTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAACGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGTAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAACCGAGCGCAATATTCAATATAGAGCACAATCGGTCCCCTCGcc                                                                                                                                                                                                                       >  1:115859/1‑275 (MQ=255)
                                                             acgGCTGGCAGCGTGAAGTGGGGGACGATCGGCGCAGACAATATGCCGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                                                                             <  4:142331/275‑1 (MQ=255)
                                                                        cGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCTTCGCCCCTCACCCCAATGGGGAGAGGGGGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGGGAGGGTTAgg                                                                                                                                                                                                   >  2:190133/1‑274 (MQ=255)
                                                                            aaGAGAGCGACAGACGGCGCGGACAATAGGACGACCGCCGACGCCAGATTGAAGATAGCCGCAGGCAGCTGGACGAGCGCCAGCGCAAGCTTGATCAGGAGCGTCGGAAAATGGGGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                                                                               <  2:434201/274‑1 (MQ=255)
                                                                                         tCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGAGCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAgcgcg                                                                                                                                                                                  <  1:153499/274‑1 (MQ=255)
                                                                                                              ccgccgccgCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAatat                                                                                                                                                                                                                                              >  4:423462/1‑193 (MQ=255)
                                                                                                                  cgccgccAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAg                                                                                                                                                                                                                                          <  2:423463/193‑1 (MQ=255)
                                                                                                                                                cTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATgg                                                                                                                                                                                                                                                                               <  5:215188/126‑1 (MQ=255)
                                                                                                                                                cTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATgg                                                                                                                                                                                                                                                                               <  5:88493/126‑1 (MQ=255)
                                                                                                                                                cTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATgg                                                                                                                                                                                                                                                                               >  6:215188/1‑126 (MQ=255)
                                                                                                                                                       aTCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCa                                                                                                                                                                                                                                                  >  8:28535/1‑148 (MQ=255)
                                                                                                                                                            cagcgccagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACg                                                                                                                                                                                                                                                               >  7:735146/1‑130 (MQ=255)
                                                                                                                                                            cagcgccagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACg                                                                                                                                                                                                                                                               <  8:735146/130‑1 (MQ=255)
                                                                                                                                                                     cTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAatat                                                                                                                                                                                                                                              >  8:1069142/1‑138 (MQ=255)
                                                                                                                                                                     cTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAatat                                                                                                                                                                                                                                              <  7:1069142/138‑1 (MQ=255)
                                                                                                                                                                       tGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCTCCCCATTGgggtgag                                                                                                                                                                                                                                                                         >  7:1001593/1‑109 (MQ=255)
                                                                                                                                                                       tGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCg                                                                                                                                                                                                                                  >  5:729029/1‑148 (MQ=255)
                                                                                                                                                                         aTCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCa                                                                                                                                                                                                                                                  >  7:993197/1‑130 (MQ=255)
                                                                                                                                                                            aGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTccc                                                                                                                                                                                                                             >  6:474525/1‑148 (MQ=255)
                                                                                                                                                                             ggAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggagagg                                                                                                                                                                                                                                                                        <  8:132638/104‑1 (MQ=255)
                                                                                                                                                                             ggAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAat                                                                                                                                                                                                                                                       >  6:547831/1‑121 (MQ=255)
                                                                                                                                                                             ggAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAat                                                                                                                                                                                                                                                       <  5:547831/121‑1 (MQ=255)
                                                                                                                                                                               aGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACg                                                                                                                                                                                                                                                               >  7:1109813/1‑111 (MQ=255)
                                                                                                                                                                                 cGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACg                                                                                                                                                                                                                                                               <  8:244806/109‑1 (MQ=255)
                                                                                                                                                                                 cGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACg                                                                                                                                                                                                                                                               >  7:244806/1‑109 (MQ=255)
                                                                                                                                                                                   tCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAatat                                                                                                                                                                                                                                              >  4:368731/1‑124 (MQ=255)
                                                                                                                                                                                   tCGGCAACTGGAGGATGAAGAGCGGACAATGGAAGATGAGTCTGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCGCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCCAGTACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCCCATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCAt                                                                                       >  2:284700/1‑275 (MQ=255)
                                                                                                                                                                                       cAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAg                                                                                                                                                                                                                                          <  2:368732/124‑1 (MQ=255)
                                                                                                                                                                                        aaCTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTg                                                                                                                                                                                                                 <  6:893054/148‑1 (MQ=255)
                                                                                                                                                                                         aCTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATgg                                                                                                                                                                                                                                                                               >  6:971653/1‑85 (MQ=255)
                                                                                                                                                                                         aCTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATgg                                                                                                                                                                                                                                                                               <  5:971653/85‑1 (MQ=255)
                                                                                                                                                                                         aCTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCg                                                                                                                                                                                                                                            <  8:104908/120‑1 (MQ=255)
                                                                                                                                                                                             gaggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCTCCCCATTGgggtgagg                                                                                                                                                                                                                                                                        <  6:704228/88‑1 (MQ=255)
                                                                                                                                                                                             gaggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCTCCCCATTGgggtgagg                                                                                                                                                                                                                                                                        >  5:704228/1‑88 (MQ=255)
                                                                                                                                                                                             gaggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCgg                                                                                                                                                                                                                                 >  7:7341/1‑127 (MQ=255)
                                                                                                                                                                                             gaggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCgg                                                                                                                                                                                                                                 <  8:7341/127‑1 (MQ=255)
                                                                                                                                                                                             gaggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                                                                            <  6:729029/148‑1 (MQ=255)
                                                                                                                                                                                              aggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                                                                                           <  6:580485/148‑1 (MQ=255)
                                                                                                                                                                                               ggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCa                                                                                                                                                                                                                                                  >  7:208764/1‑108 (MQ=255)
                                                                                                                                                                                               ggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCa                                                                                                                                                                                                                                                  <  8:208764/108‑1 (MQ=255)
                                                                                                                                                                                                 aTGAAGAGCGGAGTATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCg                                                                                                                                                                                                                                            >  5:845182/1‑112 (MQ=255)
                                                                                                                                                                                                 aTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCg                                                                                                                                                                                                                                            <  6:845182/112‑1 (MQ=255)
                                                                                                                                                                                                    aaGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGa                                                                                                                                                                                                                                                              >  6:16518/1‑91 (MQ=255)
                                                                                                                                                                                                       agCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGc                                                                                                                                                                                                                        >  6:483216/1‑126 (MQ=255)
                                                                                                                                                                                                       agCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGc                                                                                                                                                                                                                        <  5:483216/126‑1 (MQ=255)
                                                                                                                                                                                                         cGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGt                                                                 <  2:205696/275‑1 (MQ=255)
                                                                                                                                                                                                          ggAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggagaggg                                                                                                                                                                                                                                                                       <  8:624278/76‑1 (MQ=255)
                                                                                                                                                                                                              aaTGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGCGCGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCATCCCAATCCCCATCCCATCAAACCCCACCTCAAACCCTGACCGTAGTCCGTTCTCCTCCACCAAACCCTCGAACTGTTGCCCGTTATGGGTCAGaa                                                            >  4:222149/1‑275 (MQ=255)
                                                                                                                                                                                                                 ggaggaTGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGGGAAAACGAGCGCAATATTCAATATc                                                                                                                                                                                                                                             >  5:956824/1‑95 (MQ=255)
                                                                                                                                                                                                                 ggAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGGGAAAACGAGCGCAATATTCAATATc                                                                                                                                                                                                                                             <  6:956824/95‑1 (MQ=255)
                                                                                                                                                                                                                  gAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCg                                                                                                                                                                                                                                            <  7:825337/95‑1 (MQ=255)
                                                                                                                                                                                                                  gAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                                                                                           >  6:880694/1‑128 (MQ=255)
                                                                                                                                                                                                                     gATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAatat                                                                                                                                                                                                                                              <  7:1083829/90‑1 (MQ=255)
                                                                                                                                                                                                                     gATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAatat                                                                                                                                                                                                                                              >  8:1083829/1‑90 (MQ=255)
                                                                                                                                                                                                                         aGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGagagt                                                                                                                                                                                                         <  6:893004/123‑2 (MQ=255)
                                                                                                                                                                                                                         aGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                                                                                         >  5:893004/1‑123 (MQ=255)
                                                                                                                                                                                                                            aTGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGc                                                                                                                                                              <  1:138279/163‑1 (MQ=255)
                                                                                                                                                                                                                                gCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCg                                                                                                                                                                                                                                  <  8:830897/91‑1 (MQ=255)
                                                                                                                                                                                                                                   aTGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                                                                                         >  8:83291/1‑113 (MQ=255)
                                                                                                                                                                                                                                   aTGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                                                                                         <  7:83291/113‑1 (MQ=255)
                                                                                                                                                                                                                                     gAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                                                                                 >  5:929606/1‑119 (MQ=255)
                                                                                                                                                                                                                                         ttGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTc                                                                                                                                                                                                                               <  5:812154/85‑1 (MQ=255)
                                                                                                                                                                                                                                               aTGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                                                                               >  6:739968/1‑111 (MQ=255)
                                                                                                                                                                                                                                                    gCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCaa                                                                                                                                                                                                                                                        <  5:1045825/49‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                                                                               <  7:357107/106‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                                                                               >  8:357107/1‑106 (MQ=255)
                                                                                                                                                                                                                                                     cAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                                                                              >  5:590565/1‑106 (MQ=255)
                                                                                                                                                                                                                                                      aGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggggg                                                                                                                                                                                          >  3:46251/1‑109 (MQ=255)
                                                                                                                                                                                                                                                       ggCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCCGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATCTGGGTTAGCATCACCCGTTCCGAGCGGCTTACCCGATTCCGCGCAAGcca                    >  2:239976/1‑272 (MQ=255)
                                                                                                                                                                                                                                                        gCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTg                                                                                                                                                                                                                 >  1:25362/1‑84 (MQ=255)
                                                                                                                                                                                                                                                        gCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                                                                                          >  5:445552/1‑91 (MQ=255)
                                                                                                                                                                                                                                                          aTCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAg                                                                                                                                                                                      <  1:46251/109‑1 (MQ=255)
                                                                                                                                                                                                                                                           tCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                                                                                     <  7:514054/141‑1 (MQ=255)
                                                                                                                                                                                                                                                           tCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                                                                                     >  8:514054/1‑141 (MQ=255)
                                                                                                                                                                                                                                                            cGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                                                                             <  3:25362/84‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ctctCCCCATTGGGGTGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                                                                            <  7:178392/80‑1 (MQ=38)
                                                                                                                                                                                                                                                                 ctctCCCCATTGGGGTGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                                                                            >  8:178392/1‑80 (MQ=38)
                                                                                                                                                                                                                                                                 ctctCCCCATTGGGGTGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAtct                                                                                                                                                   <  7:271630/137‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAatat                                                                                                                                                                                                                                              >  5:47112/1‑46 (MQ=255)
                                                                                                                                                                                                                                                                 cTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAatat                                                                                                                                                                                                                                              <  6:47112/46‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTcc                                                                                                                                                                          >  5:131470/1‑113 (MQ=255)
                                                                                                                                                                                                                                                                  tCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTcc                                                                                                                                                                          <  6:131470/113‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                                                                                     >  7:271055/1‑134 (MQ=255)
                                                                                                                                                                                                                                                                  tCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                                                                                     <  8:271055/134‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATc        >  2:422002/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                       ccAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggggg                                                                                                                                                                                          <  7:601617/92‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aaTGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCTTGGACAGCGGATAACCTGATTCAGCGCAAGCCCGGTATTTAAATCACAGTg  >  4:91218/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                          aTGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                                                                                          >  5:846424/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                          aTGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                                                                                          <  6:846424/73‑1 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                
CAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTG  >  NC_000913/4313967‑4314505

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