Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,314,295 T→A 5.0% intergenic (+98/+49) 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,2950TA 5.0% 563.5 / 14.8 160intergenic (+98/+49)yjdP/phnPputative periplasmic protein/5‑phospho‑alpha‑D‑ribosyl 1,2‑cyclic phosphate phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (62/89);  new base A (5/3);  total (68/92)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 168 total aligned reads displayed.

TCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCA  >  NC_000913/4314146‑4314444
                                                                                                                                                     |                                                                                                                                                     
  ggCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCctct                                                                                                                                                     >  2:1980410/1‑150 (MQ=255)
    cAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                   >  2:953971/1‑150 (MQ=255)
    cAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                   <  2:2061344/150‑1 (MQ=255)
    cAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                   <  1:1622802/150‑1 (MQ=255)
    cAACGGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGAGGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                   <  2:1169325/150‑1 (MQ=255)
      aCTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                 >  1:688124/1‑150 (MQ=255)
        tGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                               >  1:904075/1‑150 (MQ=255)
         ggaggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                              >  1:1545633/1‑150 (MQ=255)
          gaggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                             <  2:2415553/150‑1 (MQ=255)
          gaggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                             <  1:1125050/150‑1 (MQ=255)
                     gCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                  >  2:1630319/1‑150 (MQ=255)
                      cGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                 >  1:676742/1‑150 (MQ=255)
                       ggAGAATTGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                <  2:1517222/150‑1 (MQ=255)
                         agaATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                              >  1:1736207/1‑150 (MQ=255)
                          gaATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggggg                                                                                                                             <  2:637923/150‑1 (MQ=255)
                          gaATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggggg                                                                                                                             <  2:601309/150‑1 (MQ=255)
                             tGGAAGATGAGTAGGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCCCTGGGGAGAGGGAGAAAACGAGCGCAAGATTCAATATCGAGCACAATCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCa                                                                                                                          <  2:1153003/150‑1 (MQ=39)
                                          tGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTcc                                                                                                             >  2:2149342/1‑150 (MQ=255)
                                          tGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTcc                                                                                                             >  2:2149269/1‑150 (MQ=255)
                                           gggCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTccc                                                                                                            >  1:2340756/1‑150 (MQ=255)
                                            ggCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCct                                                                                                           >  2:2114536/1‑150 (MQ=255)
                                            ggCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCct                                                                                                           <  2:1110533/150‑1 (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:1090744/150‑1 (MQ=255)
                                                          gTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCAcc                                                                                             >  1:81568/1‑150 (MQ=255)
                                                              ggAGCCGGCATCGCCCCTCACCCCAAGGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGAACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCaa                                                                                         <  1:2536316/150‑1 (MQ=255)
                                                               gcgCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGTGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                        <  1:2396055/148‑1 (MQ=255)
                                                               gAGCAGGCATCGCCCCTCACCCCTTTGGGTTGTGGGAGAAATCGAGCTAAATATTCAATATCGAGCACTATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                        <  1:2149269/150‑1 (MQ=37)
                                                               gAGCAGGCATCGCCCCTCACCCCATTGGGGAGAGGGAGAAACCGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                        <  1:2218986/150‑1 (MQ=255)
                                                               gAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                        <  1:2149342/150‑1 (MQ=255)
                                                               gAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                        >  2:663350/1‑150 (MQ=255)
                                                               gAGCAGGCATCGCACCTCACCCCATTGGGTAGTGGGATAAAACGAGCGAAATATTAAATATCGAGCAAAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                        <  1:1249517/150‑1 (MQ=38)
                                                                  cAGGCATCGCCCCTCACCCCAAGGGGGAGAGGGAGACAACGAGCGCACTAGTCAATATCGAGCACAATCGGTCCCCTCGCCCCGCGGGGGAGCGGGGGAGGGTGAGGGGGCCCGCGCGCACATTCCCGCGCCTCACGCCACCCCaagcgc                                                                                     >  2:1395236/1‑146 (MQ=255)
                                                                     gCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccat                                                                                  >  1:2230423/1‑150 (MQ=255)
                                                                     gCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccat                                                                                  <  1:1753083/150‑1 (MQ=255)
                                                                        tCGCCCCTCACCCCACTGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGCGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatccc                                                                               <  2:2275228/150‑1 (MQ=255)
                                                                        gcGCCCCTCACCCCAATGGGGAGAGGGAGACCACGAGCGCAAGATTCAATCTCGAGCACACTCGGTCCCCCCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatccc                                                                               <  2:1244490/149‑1 (MQ=39)
                                                                             ccTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCaa                                                                          <  2:2120792/150‑1 (MQ=255)
                                                                               tCACCCCAATGGCGAGAGGGAGAAAACGAGCGCAATATTCCATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGATAGGGTGAGGGGGCCCGCCCGCACATTCCCTCTCAACACGCCACCCCAATCTCCATCCCATCAAAc                                                                        >  2:2653988/1‑150 (MQ=38)
                                                                                 aCCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCCGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAccc                                                                      >  2:744000/1‑150 (MQ=255)
                                                                                 aCCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAccc                                                                      <  2:2279369/150‑1 (MQ=255)
                                                                                 aCCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAccc                                                                      <  2:2652945/150‑1 (MQ=255)
                                                                                 aCCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAccc                                                                      <  1:1041571/150‑1 (MQ=255)
                                                                                   cccAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCa                                                                    <  1:369793/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)
                                                                                                 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:1526306/1‑150 (MQ=255)
                                                                                                  gagaAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACg                                                     >  2:1769069/1‑150 (MQ=255)
                                                                                                     aaaaCGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTa                                                  >  2:1448884/1‑150 (MQ=255)
                                                                                                      aaaCGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAg                                                 <  1:953971/150‑1 (MQ=255)
                                                                                                       aaCGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGt                                                <  1:2459674/150‑1 (MQ=255)
                                                                                                          gAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCg                                             >  1:1462210/1‑150 (MQ=255)
                                                                                                          gAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCg                                             >  1:2085995/1‑150 (MQ=255)
                                                                                                             cgcAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGttt                                          >  1:1164254/1‑150 (MQ=255)
                                                                                                                 ctagtCAATATCGCGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCa                                      <  1:1278632/146‑1 (MQ=255)
                                                                                                                 atatTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCa                                      <  1:789139/150‑1 (MQ=255)
                                                                                                                 atatTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACAGTCCCTCTCATCACGCCACCCCACTCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCa                                      <  1:314898/150‑1 (MQ=255)
                                                                                                                  tatTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCAt                                     <  1:1143196/150‑1 (MQ=255)
                                                                                                                    ttCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCa                                   <  1:2060890/150‑1 (MQ=255)
                                                                                                                    ttCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCa                                   <  1:712290/150‑1 (MQ=255)
                                                                                                                    ttCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCa                                   <  1:259743/150‑1 (MQ=255)
                                                                                                                    ttCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCa                                   <  1:259784/150‑1 (MQ=255)
                                                                                                                    ttCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACAGTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCa                                   <  1:645017/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:2631986/150‑1 (MQ=255)
                                                                                                                           tCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGc                            >  1:1598552/1‑150 (MQ=255)
                                                                                                                            cGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCa                           <  2:928094/150‑1 (MQ=255)
                                                                                                                               gCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCa                        >  2:723231/1‑150 (MQ=255)
                                                                                                                               gCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTACATCAGCCACGCATCa                        >  2:2482899/1‑150 (MQ=255)
                                                                                                                                     tCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtgg                  <  2:2096614/150‑1 (MQ=255)
                                                                                                                                      cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGc                                 >  1:679970/6‑134 (MQ=255)
                                                                                                                                      cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGc                                 <  2:679970/129‑1 (MQ=255)
                                                                                                                                      cGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACAGTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggt                 <  1:961487/150‑1 (MQ=255)
                                                                                                                                      cGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACAGTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggt                 <  1:477963/150‑1 (MQ=255)
                                                                                                                                        taacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGGCCCGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCACAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggtgg               >  1:759652/4‑150 (MQ=255)
                                                                                                                                        ggccccTCGCCCCTCTGGGGAGAGGGTTCGGGTGCGGGGGCCAGCGCGCACCTTCCCTCTCATCACGCCACCCCACTCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggtgg               <  1:2113224/148‑1 (MQ=255)
                                                                                                                                        gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtggtgg               <  1:2724013/150‑1 (MQ=255)
                                                                                                                                          acccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGGGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCt             >  2:1658288/2‑150 (MQ=255)
                                                                                                                                            ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGa           <  1:2091006/150‑1 (MQ=255)
                                                                                                                                            ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGa           <  1:2091086/150‑1 (MQ=255)
                                                                                                                                             cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGat          >  1:1970952/1‑150 (MQ=255)
                                                                                                                                                 ccccTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATgg      <  1:2232641/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)
                                                                                                                                                     |                                                                                                                                                     
TCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCA  >  NC_000913/4314146‑4314444

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: