Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,566,544 C→A S13I (AGC→ATC)  glgP ← glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,566,5440CA99.6% 1054.1 / NA 280S13I (AGC→ATC) glgPglycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (147/132);  minor base T (0/1);  total (147/133)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 280 total aligned reads displayed.

GGTTTGAACGTAACCAGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGC  >  NC_000913/3566398‑3566690
                                                                                                                                                  |                                                                                                                                                  
ggTTTGAACGTAACCAGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATaa                                                                                                                                                 >  1:2440203/1‑150 (MQ=255)
ggTTTGAACGTAACCAGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATaa                                                                                                                                                 <  2:2527541/150‑1 (MQ=255)
  tttGAACGTAACCAGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGc                                                                                                                                               >  1:1067680/1‑150 (MQ=255)
     gAACGTAACCAGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTg                                                                                                                                            >  2:2317161/1‑150 (MQ=255)
       aCGTAACCAGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTggg                                                                                                                                          >  2:1183346/1‑150 (MQ=255)
        cGTAACCAGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGc                                                                                                                                         >  2:1511340/1‑150 (MQ=255)
              cAGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGaa                                                                                                                                   >  1:1295126/1‑150 (MQ=255)
               aGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAat                                                                                                                                  <  1:2329694/150‑1 (MQ=255)
               aGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAat                                                                                                                                  <  1:15479/150‑1 (MQ=255)
               aGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAat                                                                                                                                  <  2:111340/150‑1 (MQ=255)
                 cgcTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAatat                                                                                                                                <  2:2796911/150‑1 (MQ=255)
                  gcTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCATTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATg                                                                                                                               >  2:1099505/1‑150 (MQ=255)
                   cTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCTTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGt                                                                                                                              <  2:2382548/150‑1 (MQ=255)
                   cTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGt                                                                                                                              >  1:2112939/1‑150 (MQ=255)
                   cTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGt                                                                                                                              <  2:1276986/150‑1 (MQ=255)
                       aCGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAAc                                                                                                                          >  1:2900705/1‑150 (MQ=255)
                       aCGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAAc                                                                                                                          <  2:678796/150‑1 (MQ=255)
                              gATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTATGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCAt                                                                                                                   >  2:2620513/1‑150 (MQ=255)
                              gATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCAt                                                                                                                   >  1:2572531/1‑150 (MQ=255)
                              gATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCAt                                                                                                                   >  1:1658546/1‑150 (MQ=255)
                              gATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCAt                                                                                                                   >  2:704738/1‑150 (MQ=255)
                               aTCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACgg                                                                                                                        <  1:2381321/144‑1 (MQ=255)
                               aTCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCAtt                                                                                                                  <  1:61459/150‑1 (MQ=255)
                               aTCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCAtt                                                                                                                  >  1:2853539/1‑150 (MQ=255)
                                tCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTc                                                                                                                 <  2:2774510/150‑1 (MQ=255)
                                    gcACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCatat                                                                                                             <  2:185217/150‑1 (MQ=255)
                                      acaGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAg                                                                                                           >  1:965404/1‑150 (MQ=255)
                                        aGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGc                                                                                                         <  1:294653/150‑1 (MQ=255)
                                          cAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGt                                                                                                       >  1:2836706/1‑150 (MQ=255)
                                           aaaTAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGtt                                                                                                      >  1:2167955/1‑150 (MQ=255)
                                           aaaTAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGtt                                                                                                      <  1:2657601/150‑1 (MQ=255)
                                            aaTAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGttt                                                                                                     >  1:1554441/1‑150 (MQ=255)
                                             aTAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTc                                                                                                    >  2:958782/1‑150 (MQ=255)
                                             aTAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTc                                                                                                    >  2:2360030/1‑150 (MQ=255)
                                             aTAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTc                                                                                                    >  2:2284149/1‑150 (MQ=255)
                                              tAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTcc                                                                                                   <  2:1562444/150‑1 (MQ=255)
                                               aaCGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCt                                                                                                  >  2:493210/1‑150 (MQ=255)
                                                  gTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                                                                               <  1:2775201/150‑1 (MQ=255)
                                                  gTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                                                                               <  2:2767581/150‑1 (MQ=255)
                                                     gCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAAc                                                                                            <  2:1658546/150‑1 (MQ=255)
                                                       gTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTa                                                                                          >  1:1299047/1‑150 (MQ=255)
                                                          cAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAttt                                                                                       <  2:1757799/150‑1 (MQ=255)
                                                           aGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTc                                                                                      <  1:1418242/150‑1 (MQ=255)
                                                            gCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCa                                                                                     >  1:1197476/1‑150 (MQ=255)
                                                            gCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCa                                                                                     <  1:1130473/150‑1 (MQ=255)
                                                                 tCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGa                                                                                >  1:1744957/1‑150 (MQ=255)
                                                                 tCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGa                                                                                <  1:162864/150‑1 (MQ=255)
                                                                   aTGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATa                                                                              >  2:91375/1‑150 (MQ=255)
                                                                    tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAg                                                                             >  1:1961756/1‑150 (MQ=255)
                                                                     gTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGt                                                                            >  2:1781909/1‑150 (MQ=255)
                                                                        tATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaaa                                                                         <  2:1889212/150‑1 (MQ=255)
                                                                           tGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCt                                                                      <  1:1243637/150‑1 (MQ=255)
                                                                             gCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCa                                                                    >  1:404294/1‑150 (MQ=255)
                                                                             gCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCa                                                                    >  2:2239240/1‑150 (MQ=255)
                                                                               gacgacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACg                                                                  >  2:322517/1‑150 (MQ=255)
                                                                                acgacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACgg                                                                 <  2:1295126/150‑1 (MQ=255)
                                                                                acgacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACgg                                                                 <  2:1818047/150‑1 (MQ=255)
                                                                                  gacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTGTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTa                                                               >  1:768887/1‑150 (MQ=255)
                                                                                   acCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTAc                                                              <  1:520901/150‑1 (MQ=255)
                                                                                     cGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGa                                                            >  1:1459170/1‑150 (MQ=255)
                                                                                       ggTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACt                                                          >  1:1018125/1‑150 (MQ=255)
                                                                                         tCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTc                                                        <  2:1698182/150‑1 (MQ=255)
                                                                                         tCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTc                                                        <  2:1937264/150‑1 (MQ=255)
                                                                                            tttCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcgcc                                                     >  1:1594690/1‑150 (MQ=255)
                                                                                            tttCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcgcc                                                     >  2:1359272/1‑150 (MQ=255)
                                                                                            tttCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcgcc                                                     >  1:646673/1‑150 (MQ=255)
                                                                                              tCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcgccgc                                                   <  2:592770/150‑1 (MQ=255)
                                                                                               ccAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCcgcg                                                  >  1:2422278/1‑150 (MQ=255)
                                                                                                 aaTCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGAc                                                <  2:2853539/150‑1 (MQ=255)
                                                                                                   tCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCt                                              <  2:2435080/150‑1 (MQ=255)
                                                                                                   tCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCt                                              <  1:720102/150‑1 (MQ=255)
                                                                                                   tCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCt                                              >  1:483777/1‑150 (MQ=255)
                                                                                                     gTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGc                                            <  1:1856463/150‑1 (MQ=255)
                                                                                                     gTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGc                                            <  2:1038333/150‑1 (MQ=255)
                                                                                                      tAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGcc                                           >  1:1721376/1‑150 (MQ=255)
                                                                                                      tAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGcc                                           >  1:2688843/1‑150 (MQ=255)
                                                                                                              aGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTaaaa                                   <  1:2065588/150‑1 (MQ=255)
                                                                                                                cTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATc                                 >  2:2088781/1‑150 (MQ=255)
                                                                                                                cTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATc                                 >  2:2739267/1‑150 (MQ=255)
                                                                                                                cTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATc                                 >  1:130755/1‑150 (MQ=255)
                                                                                                                     aaGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTg                            >  2:2962324/1‑150 (MQ=255)
                                                                                                                     aaGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTg                            <  2:1603146/150‑1 (MQ=255)
                                                                                                                     aaGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTg                            <  2:2058568/150‑1 (MQ=255)
                                                                                                                      aGCGATAGTGTTCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAATCGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGc                           <  2:2893451/150‑1 (MQ=255)
                                                                                                                      aGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGc                           <  2:1471442/150‑1 (MQ=255)
                                                                                                                          aTAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATa                       >  2:627461/1‑150 (MQ=255)
                                                                                                                            agagTGCTTAAGAGCTTCTACGTTAATCGTGGGCGATGAATATGTAAACGGAGAATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCGGCCAGCTAAAATCCATTGCCATAGc                     <  2:79365/150‑1 (MQ=255)
                                                                                                                            agagTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGc                     <  1:2280520/150‑1 (MQ=255)
                                                                                                                             gagTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGcc                    >  2:2235764/1‑150 (MQ=255)
                                                                                                                              agTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCt                   >  1:2008456/1‑150 (MQ=255)
                                                                                                                              agTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCt                   >  1:1816392/1‑150 (MQ=255)
                                                                                                                                  cTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACg               >  1:1772862/1‑150 (MQ=255)
                                                                                                                                      agagCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGc           >  2:1017289/1‑150 (MQ=255)
                                                                                                                                       gagCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGca          >  2:1306998/1‑150 (MQ=255)
                                                                                                                                       gagCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGca          <  2:1948343/150‑1 (MQ=255)
                                                                                                                                       gagCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGca          <  2:2139986/150‑1 (MQ=255)
                                                                                                                                            tCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAAc     >  1:853476/1‑150 (MQ=255)
                                                                                                                                             cTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAAcc    >  1:219982/1‑150 (MQ=255)
                                                                                                                                               aCGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGc  <  2:1299047/150‑1 (MQ=255)
                                                                                                                                               aCGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGc  >  1:1838851/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGTTTGAACGTAACCAGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGC  >  NC_000913/3566398‑3566690

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