Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,314,334 T→G 9.3% intergenic (+137/+10) 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,3340TG 9.3% 714.7 / 10.8 258intergenic (+137/+10)yjdP/phnPputative periplasmic protein/5‑phospho‑alpha‑D‑ribosyl 1,2‑cyclic phosphate phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (122/106);  new base G (3/27);  total (125/133)
Fisher's exact test for biased strand distribution p-value = 3.38e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

TGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGC  >  NC_000913/4314188‑4314479
                                                                                                                                                  |                                                                                                                                                 
  ggCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCct                                                                                                                                              <  2:1768961/150‑1 (MQ=255)
   ggtttgtggttggttATGGAGCAGGCATCGCCCCTCACCCCAACGGGGAGAGGGAGAAAACGAGCGCAATATTAAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGTGTGGTTTTGGGTGTGGGGGCCAGCGCGCACATTCCctc                                                                                                                                             <  1:2443830/143‑1 (MQ=255)
   gCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCCCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCctc                                                                                                                                             >  2:817511/1‑150 (MQ=255)
   gCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAAGATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCctc                                                                                                                                             >  2:3572835/1‑150 (MQ=255)
    cGATGCGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCctct                                                                                                                                            <  1:436257/150‑1 (MQ=255)
    cGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCctct                                                                                                                                            <  2:3419463/150‑1 (MQ=255)
    cGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCctct                                                                                                                                            >  2:3266867/1‑150 (MQ=255)
     gttGAGGTTGGGTTTGGAGCAGGCATCGCCCCCCACCCCAATGGGGAGAGGTAGAAAACGAGCGCAATATTCAATATCGAGCAAAATCGGTCCCCTCGCCCCTCTGGGGAGTGGTTTTGGGTGAGGGGGCCAGCGCGCACATTCCCtctc                                                                                                                                           <  1:1214028/148‑1 (MQ=255)
       tGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCtcat                                                                                                                                         <  2:586659/150‑1 (MQ=255)
                 tATGGAGCAGGCAGCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCAccc                                                                                                                               <  1:1612934/150‑1 (MQ=255)
                     gAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACAGTCCCTCTCATCACGCCACCCCact                                                                                                                           >  1:655773/1‑148 (MQ=255)
                      aGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAtc                                                                                                                          >  1:784379/1‑150 (MQ=255)
                        cAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAtctc                                                                                                                        >  2:2238469/1‑150 (MQ=255)
                         aGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAGTATTCAATATCGAGCACAATCGGACCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACAGTCCCGCTCATCACGCCCCCCCAGTCTcc                                                                                                                       >  1:2819165/1‑150 (MQ=38)
                         aGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTcc                                                                                                                       >  1:4274182/1‑150 (MQ=255)
                         aGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTcc                                                                                                                       <  2:2962014/150‑1 (MQ=255)
                         aGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTcc                                                                                                                       <  2:1476683/150‑1 (MQ=255)
                           gCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccat                                                                                                                     >  2:4009393/1‑150 (MQ=255)
                            cATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatc                                                                                                                    <  1:1235073/150‑1 (MQ=255)
                             aTCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatcc                                                                                                                   <  1:2990897/150‑1 (MQ=255)
                              tCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatccc                                                                                                                  >  2:3613738/1‑150 (MQ=255)
                               cGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatccca                                                                                                                 <  2:1596313/150‑1 (MQ=255)
                               cGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatccca                                                                                                                 <  1:1948864/150‑1 (MQ=255)
                                 tcccTCTCCCCATTGGGGTGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatcccatc                                                                                                               >  2:2187621/2‑150 (MQ=255)
                                 ccccTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatcccatc                                                                                                               >  2:3421103/1‑150 (MQ=255)
                                    cTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCaaa                                                                                                            >  2:3309423/1‑150 (MQ=255)
                                     tCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAc                                                                                                           <  2:112016/150‑1 (MQ=255)
                                     tCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAc                                                                                                           <  2:2414189/150‑1 (MQ=255)
                                         cccAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCa                                                                                                       >  1:382100/1‑150 (MQ=255)
                                           cAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCAcc                                                                                                     >  2:3500135/1‑150 (MQ=255)
                                           cAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCAcc                                                                                                     >  2:288139/1‑150 (MQ=255)
                                            aaTGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCt                                                                                                    >  2:3202064/1‑150 (MQ=255)
                                            aaTGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCt                                                                                                    <  1:2061975/150‑1 (MQ=255)
                                                gggagagggagaAAACGAGCGCACTATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCACTCTCCATCCCATCAAACCCCACCTCaaa                                                                                                <  2:634293/150‑1 (MQ=255)
                                                gggagagggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACAGTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCaaa                                                                                                <  2:2792011/150‑1 (MQ=255)
                                                 ggagagggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAAc                                                                                               <  1:798796/150‑1 (MQ=255)
                                                 ggagagggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAAc                                                                                               <  2:3893631/150‑1 (MQ=255)
                                                 ggagagggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAAc                                                                                               <  1:2532085/150‑1 (MQ=255)
                                                  gagagggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAAcc                                                                                              >  1:1008396/1‑150 (MQ=255)
                                                     agggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTg                                                                                           >  1:3268542/1‑150 (MQ=255)
                                                       ggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGAc                                                                                         >  1:2145282/1‑150 (MQ=255)
                                                       ggagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGAc                                                                                         >  1:4231837/1‑150 (MQ=255)
                                                        gagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACg                                                                                        <  1:3645706/150‑1 (MQ=255)
                                                                gAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCg                                                                                <  2:3035456/150‑1 (MQ=255)
                                                                   cgcAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGttt                                                                             >  2:4171041/1‑150 (MQ=255)
                                                                          ttCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCa                                                                      <  2:1242590/150‑1 (MQ=255)
                                                                          ttCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCa                                                                      <  2:3144192/150‑1 (MQ=255)
                                                                          ttCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACAGTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCa                                                                      <  1:3455787/150‑1 (MQ=255)
                                                                          ttCAATATCGAGCACAAACGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCa                                                                      >  1:2591865/1‑150 (MQ=255)
                                                                          tacAATATCCACAAAAAACGCTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCa                                                                      <  1:1262560/148‑1 (MQ=38)
                                                                             aaTATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGcc                                                                   >  1:1885449/1‑150 (MQ=255)
                                                                                atCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACg                                                                <  2:3807710/150‑1 (MQ=255)
                                                                                   gAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCAt                                                             >  2:3986333/1‑150 (MQ=255)
                                                                                    aGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATc                                                            >  1:4149326/1‑150 (MQ=255)
                                                                                    aGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATc                                                            <  1:1544511/150‑1 (MQ=255)
                                                                                     gCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCa                                                           >  2:311876/1‑150 (MQ=255)
                                                                                         aaTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACt                                                       >  2:3383885/1‑150 (MQ=255)
                                                                                            cGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggt                                                    <  2:4256493/150‑1 (MQ=255)
                                                                                             ggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACAGTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggtg                                                   <  1:887082/150‑1 (MQ=255)
                                                                                             ctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggtg                                                   >  1:3601907/5‑150 (MQ=255)
                                                                                              gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggtgg                                                  >  1:53910/1‑150 (MQ=255)
                                                                                               tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGc                                                 <  1:1857226/150‑1 (MQ=255)
                                                                                               tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGc                                                 <  1:2965176/150‑1 (MQ=255)
                                                                                                   cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGat                                             >  1:999947/1‑150 (MQ=255)
                                                                                                   cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGat                                             >  1:3841201/1‑150 (MQ=255)
                                                                                                    tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGata                                            >  1:3670454/1‑150 (MQ=255)
                                                                                                      gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATg                                          >  1:1163364/1‑150 (MQ=255)
                                                                                                        cccTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATggg                                        <  1:1742672/150‑1 (MQ=255)
                                                                                                         ccTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCACTCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGt                                       <  2:2541631/150‑1 (MQ=255)
                                                                                                           tctGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCa                                     <  2:870763/150‑1 (MQ=255)
                                                                                                              ggggAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGaa                                  <  1:3940589/150‑1 (MQ=255)
                                                                                                                 gagaGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCa                               <  2:2024926/150‑1 (MQ=255)
                                                                                                                    agGGTTAGGGTGAGGGGGCCAGCGCGCTCAGTCCCTCTCATCACGCCACCCCACTCTCCATCCCATCAAAGGGCACCTGTTTCGGTGACGGTTGTGCGTTTTCCATCAGCCTCGCATCTAACTGGTGGCTGTTATGGGTCAGAATCAccc                            <  2:3738383/150‑1 (MQ=255)
                                                                                                                    agGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCAccc                            <  1:3577590/150‑1 (MQ=255)
                                                                                                                      ggTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGt                          >  1:1058476/1‑150 (MQ=255)
                                                                                                                      ggTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGt                          >  2:1501179/1‑150 (MQ=255)
                                                                                                                       gTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTg                         >  2:1808862/1‑150 (MQ=255)
                                                                                                                         tAGGGTGAGGGGGCCAGCGCGCACAGTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGc                       <  1:1333793/150‑1 (MQ=255)
                                                                                                                            ggTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAg                    <  1:930621/150‑1 (MQ=255)
                                                                                                                              tGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCg                  <  2:2712169/150‑1 (MQ=255)
                                                                                                                              ggcgggggCCAGCGCGCACCGGCCCTCTCATCACGCCACCCCACTCTCCCTCCCAGCAAACCCCACCTCACCCCCTGACGGTAGTGCGTTTTCCCTCAGCCACGCATCAACCTGGTGGCGGCTATGGGTCAGAATCACCCGTGGCGAGCg                  <  1:509223/147‑1 (MQ=255)
                                                                                                                                 gggggCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGAt               >  1:3235162/1‑150 (MQ=255)
                                                                                                                                 gggggCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGAt               >  1:1348662/1‑150 (MQ=255)
                                                                                                                                 gggggCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGAt               >  1:2840642/1‑150 (MQ=255)
                                                                                                                                 gggggCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGAt               <  1:146945/150‑1 (MQ=255)
                                                                                                                                 gggggCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGGGCGTTTTCCATCAGCCACGCATCCAACTGGTGGCTGATAGGGGTCAGAATCACCCGTGGCGAGCGGAt               >  2:1714597/1‑150 (MQ=255)
                                                                                                                                  ggggCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATa              >  1:3026218/1‑150 (MQ=255)
                                                                                                                                  ggggCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATa              >  2:2368356/1‑150 (MQ=255)
                                                                                                                                   gggCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATaa             >  1:2288239/1‑150 (MQ=255)
                                                                                                                                   gggCCAGCGCGCACAGTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCACACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATaa             <  2:1191316/150‑1 (MQ=255)
                                                                                                                                    ggCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAAc            >  2:2622079/1‑150 (MQ=255)
                                                                                                                                     gCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAAcc           >  1:1683288/1‑150 (MQ=255)
                                                                                                                                       cAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTg         >  1:427648/1‑150 (MQ=255)
                                                                                                                                       cAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTg         >  1:2898145/1‑150 (MQ=255)
                                                                                                                                        aGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGa        >  1:2687458/1‑150 (MQ=255)
                                                                                                                                        aGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGa        <  2:500088/150‑1 (MQ=255)
                                                                                                                                        aGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGa        <  2:260310/150‑1 (MQ=255)
                                                                                                                                        aGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGa        <  2:75438/150‑1 (MQ=255)
                                                                                                                                        aGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGa        <  1:3370500/150‑1 (MQ=255)
                                                                                                                                              cacaGTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAgc  <  1:1106145/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGC  >  NC_000913/4314188‑4314479

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: