Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,314,302 A→C 5.2% intergenic (+105/+42) 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,3020AC 5.2% 496.1 / 17.5 154intergenic (+105/+42)yjdP/phnPputative periplasmic protein/5‑phospho‑alpha‑D‑ribosyl 1,2‑cyclic phosphate phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (65/79);  new base C (5/3);  total (70/84)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

Only 100 of 165 total aligned reads displayed.

CTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCAC  >  NC_000913/4314153‑4314451
                                                                                                                                                     |                                                                                                                                                     
cTGGAGGAGGAAGAGCGGAGAAGGGAAGATGAGTCTGGGCGATGAGGTGGGGTATGGAGCAGGCATCGCCCCTCACCCCCATGGGGCGAGGGAGAAACCGCGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                       <  2:1028375/150‑1 (MQ=255)
 tGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                      >  1:904075/1‑150 (MQ=255)
  ggaggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                                     >  1:1545633/1‑150 (MQ=255)
  ggaggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGGATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                                     <  2:940060/150‑1 (MQ=255)
   gaggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                                    <  1:1125034/150‑1 (MQ=255)
              gCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                         >  2:1630319/1‑150 (MQ=255)
               cGGGGAATGGAAGATGTGTATGGGCGATGAGGTTGGGTATGGAGCTGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                        >  1:676538/1‑150 (MQ=255)
               cGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                        >  1:676742/1‑150 (MQ=255)
                ggAGAATTGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                       <  2:1517222/150‑1 (MQ=255)
                ggAGAATGGCAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                       <  1:1275887/150‑1 (MQ=255)
                ggAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                       <  1:2596161/150‑1 (MQ=255)
                  agaATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                                     >  1:2381986/1‑150 (MQ=255)
                  agaATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGTGAgggg                                                                                                                                     >  1:1056386/1‑150 (MQ=255)
                   gaATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAACCGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggggg                                                                                                                                    <  2:637927/150‑1 (MQ=255)
                   gaATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggggg                                                                                                                                    <  2:637923/150‑1 (MQ=255)
                      tGGAAGATGAGTAGGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCCCTGGGGAGAGGGAGAAAACGAGCGCAAGATTCAATATCGAGCACAATCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCa                                                                                                                                 <  2:1153003/150‑1 (MQ=39)
                                  aTGGGCGTTGTGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTAAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGGGCCAGCGCGCACAGTc                                                                                                                     <  1:322146/150‑1 (MQ=39)
                                  aTGGGCGTTGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTc                                                                                                                     <  1:693373/150‑1 (MQ=255)
                                   tGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTcc                                                                                                                    >  2:2149269/1‑150 (MQ=255)
                                   tGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTcc                                                                                                                    >  2:2149342/1‑150 (MQ=255)
                                    gggCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCACTATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCCGCGCGCACATTccc                                                                                                                   >  1:2219592/1‑150 (MQ=255)
                                     ggCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCACTATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCct                                                                                                                  >  2:2060890/1‑150 (MQ=255)
                                     ggCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCct                                                                                                                  >  2:2114536/1‑150 (MQ=255)
                                      gCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCctc                                                                                                                 <  2:1583220/150‑1 (MQ=255)
                                      gCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCctc                                                                                                                 >  2:2113224/1‑150 (MQ=255)
                                        gATGAGGTTGGTTAGGGAGCAGGCATCGCCCCTCACCCCAATGGGGTGTGGGAGAAAACGCGCGCAATATTCAATATCGAGCACAATCGGTCCCCGCGCCCCTCTGGGGAGTGGGTTAGGGTGTGGGGGCCAGCGCGCACATTCCCtctc                                                                                                               <  1:1262068/150‑1 (MQ=255)
                                        gATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCtctc                                                                                                               <  2:281782/150‑1 (MQ=255)
                                         aTGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCtca                                                                                                              >  2:2232641/1‑150 (MQ=255)
                                          tGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCtcat                                                                                                             <  1:1915274/150‑1 (MQ=255)
                                          tGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCtcat                                                                                                             <  1:2605372/150‑1 (MQ=255)
                                                        gcgCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGTGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                               <  1:2396055/148‑1 (MQ=255)
                                                        gAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                               <  1:2149342/150‑1 (MQ=255)
                                                           cAGGCATCGCCCCTCACCCCAAGGGGGAGAGGGAGACAACGAGCGCACTAGTCAATATCGAGCACAATCGGTCCCCTCGCCCCGCGGGGGAGCGGGGGAGGGTGAGGGGGCCCGCGCGCACATTCCCGCGCCTCACGCCACCCCaagcgc                                                                                            >  2:1395236/1‑146 (MQ=255)
                                                              gCATCGCCCCTCCCCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccat                                                                                         <  1:1753077/150‑1 (MQ=255)
                                                              gCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccat                                                                                         <  1:1753083/150‑1 (MQ=255)
                                                              gCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccat                                                                                         >  1:2230423/1‑150 (MQ=255)
                                                                aTCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatcc                                                                                       >  1:1192425/1‑150 (MQ=255)
                                                                 tCGCCCCTCACCCCACTGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGCGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatccc                                                                                      <  2:2275228/150‑1 (MQ=255)
                                                                 gcGCCCCTCACCCCAATGGGGAGAGGGAGACCACGAGCGCAAGATTCAATCTCGAGCACACTCGGTCCCCCCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatccc                                                                                      <  2:1244490/149‑1 (MQ=39)
                                                                  cGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatccca                                                                                     <  2:81568/150‑1 (MQ=255)
                                                                      ccTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCaa                                                                                 <  2:2120792/150‑1 (MQ=255)
                                                                      ccTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCaa                                                                                 >  1:726186/1‑150 (MQ=255)
                                                                        tCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAc                                                                               >  1:2314682/1‑150 (MQ=255)
                                                                          aCCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAccc                                                                             <  1:1041571/150‑1 (MQ=255)
                                                                          aCCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAccc                                                                             <  2:2652945/150‑1 (MQ=255)
                                                                          aCCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAccc                                                                             >  2:496180/1‑150 (MQ=255)
                                                                          aCCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAccc                                                                             <  2:2652949/150‑1 (MQ=255)
                                                                          aCCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAccc                                                                             <  2:2279369/150‑1 (MQ=255)
                                                                                    ggagagggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAAc                                                                   <  1:1665420/150‑1 (MQ=255)
                                                                                     gagagggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAAcc                                                                  <  1:1395181/150‑1 (MQ=255)
                                                                                      agagggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAAccc                                                                 <  1:1051953/150‑1 (MQ=255)
                                                                                      agagggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAAccc                                                                 >  1:1985569/1‑150 (MQ=255)
                                                                                          ggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGAc                                                             <  1:1728648/150‑1 (MQ=255)
                                                                                          ggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGAc                                                             >  1:585004/1‑150 (MQ=255)
                                                                                          ggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGAc                                                             >  2:1968472/1‑150 (MQ=255)
                                                                                          ggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCACCATCCCATCAAACCCCACCTCAAACCCTGAc                                                             >  2:615396/1‑150 (MQ=255)
                                                                                           gagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACg                                                            >  2:1769069/1‑150 (MQ=255)
                                                                                           gagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACg                                                            >  2:1526306/1‑150 (MQ=255)
                                                                                              aaaaCGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTa                                                         >  2:1448961/1‑150 (MQ=255)
                                                                                               aaaCGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAg                                                        <  1:953971/150‑1 (MQ=255)
                                                                                                aaCGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGt                                                       <  1:2514638/150‑1 (MQ=255)
                                                                                                   gAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCg                                                    >  1:1462210/1‑150 (MQ=255)
                                                                                                   gAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCg                                                    >  1:2085995/1‑150 (MQ=255)
                                                                                                          atatTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCa                                             <  1:789139/150‑1 (MQ=255)
                                                                                                           tatTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCAt                                            <  1:1143196/150‑1 (MQ=255)
                                                                                                           tatTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCAt                                            <  1:581683/150‑1 (MQ=255)
                                                                                                           tatTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCAt                                            <  1:1143200/150‑1 (MQ=255)
                                                                                                             ttCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCa                                          <  1:259743/150‑1 (MQ=255)
                                                                                                             ttCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCa                                          <  1:2060890/150‑1 (MQ=255)
                                                                                                               cAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTCGTGCGTTTTCCATCCGc                                        >  2:239282/1‑150 (MQ=255)
                                                                                                               cAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGc                                        >  2:1837850/1‑150 (MQ=255)
                                                                                                                  gctCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCAc                                     <  2:146067/148‑1 (MQ=255)
                                                                                                                    tCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGc                                   >  1:1598552/1‑150 (MQ=255)
                                                                                                                    tCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACAGTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGc                                   <  1:2632026/150‑1 (MQ=255)
                                                                                                                     cGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCa                                  <  2:688124/150‑1 (MQ=255)
                                                                                                                        gCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCa                               >  2:723299/1‑150 (MQ=255)
                                                                                                                        gCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTACATCAGCCACGCATCa                               >  2:2482899/1‑150 (MQ=255)
                                                                                                                            aaTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACt                           >  2:2397402/1‑150 (MQ=255)
                                                                                                                              tCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtgg                         <  2:2096614/150‑1 (MQ=255)
                                                                                                                               cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGc                                        <  2:679970/129‑1 (MQ=255)
                                                                                                                                 taacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGGCCCGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCACAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggtgg                      >  1:759652/4‑150 (MQ=255)
                                                                                                                                 taacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggtgg                      >  1:2134005/4‑150 (MQ=255)
                                                                                                                                 ggccccTCGCCCCTCTGGGGAGAGGGTTCGGGTGCGGGGGCCAGCGCGCACCTTCCCTCTCATCACGCCACCCCACTCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggtgg                      <  1:2113224/148‑1 (MQ=255)
                                                                                                                                 gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggtgg                      <  1:2724013/150‑1 (MQ=255)
                                                                                                                                 gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggtgg                      <  1:386867/150‑1 (MQ=255)
                                                                                                                                   acccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGGGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCt                    >  2:1658288/2‑150 (MQ=255)
                                                                                                                                    cccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTg                   >  1:87180/1‑150 (MQ=255)
                                                                                                                                     ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGa                  <  1:2091006/150‑1 (MQ=255)
                                                                                                                                     ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGa                  <  1:2091086/150‑1 (MQ=255)
                                                                                                                                          ccccTCTGGGGAGAGGGTGAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATGACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCATCCCCTGACGGTAGTGCGTTTGCCATCATCCACGCATCAAACTGGTTGCTGCTATgg             <  2:1800552/150‑1 (MQ=255)
                                                                                                                                           cccTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATggg            >  1:585974/1‑150 (MQ=255)
                                                                                                                                             ctctGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTc          >  2:1784826/1‑150 (MQ=255)
                                                                                                                                              tctTGGTAGAGGTTTAGGGAGACCGGGACATCGCTCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATAAAACCCCAACTCAAACCATTACCCTAGATCGTAATCCATCAGCCACGCATAAAACTTTTGTCAGATAAGTGTCa         >  2:2676672/1‑150 (MQ=255)
                                                                                                                                              tctGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCa         >  1:434372/1‑150 (MQ=255)
                                                                                                                                              tctGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCa         >  2:800536/1‑150 (MQ=255)
                                                                                                                                                  gggAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAAt     >  1:986642/1‑150 (MQ=255)
                                                                                                                                                  gggAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAAt     <  1:1299815/150‑1 (MQ=255)
                                                                                                                                                  gggAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAAt     <  1:1055623/150‑1 (MQ=255)
                                                                                                                                                    gagaGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCCAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCa   <  2:676538/150‑1 (MQ=255)
                                                                                                                                                     agagGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCAc  <  2:1342678/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCAC  >  NC_000913/4314153‑4314451

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