Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,580,643 C→A 6.0% Q943H (CAG→CAT treSB ← bifunctional trehalose synthase B/maltokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,580,6430CA 6.0% 417.9 / 9.4 130Q943H (CAG→CATtreSBbifunctional trehalose synthase B/maltokinase
Reads supporting (aligned to +/- strand):  ref base C (65/57);  new base A (2/6);  total (67/63)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 132 total aligned reads displayed.

CCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACA  >  NC_002947/4580495‑4580792
                                                                                                                                                    |                                                                                                                                                     
ccTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCt                                                                                                                                                      <  1:1820604/150‑1 (MQ=255)
ccTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCt                                                                                                                                                      >  2:1210378/1‑150 (MQ=255)
 cTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtg                                                                                                                                                     >  2:1649433/1‑150 (MQ=255)
    gTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtgctt                                                                                                                                                  <  2:2396168/150‑1 (MQ=255)
          gCTGTGTTTTGCCCGGCGCTCTTCCAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAg                                                                                                                                            >  1:1518185/1‑150 (MQ=255)
           cTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAgg                                                                                                                                           <  1:1490302/150‑1 (MQ=255)
              tgtTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGgtg                                                                                                                                        <  1:1190092/150‑1 (MQ=255)
               gtTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGgtgt                                                                                                                                       <  1:155530/150‑1 (MQ=255)
               gtTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGgtgt                                                                                                                                       <  2:1200598/150‑1 (MQ=255)
                tgttGGCCGGCGCTTCTCTAGGGGGCGGGCGGGTGCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCCTGTACACGCATAAGCAGGCCACCCTGTGCCTGCTTGCTGAGGGTGTc                                                                                                                                      <  2:92273/148‑1 (MQ=255)
                 tttGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAGGgtggcc                                                                                                                                     <  1:1815826/150‑4 (MQ=255)
                  ttGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAACCGATGATTTATGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGa                                                                                                                                    <  2:666000/150‑1 (MQ=255)
                    gCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAGGgtggcctgc                                                                                                                                  <  2:1082856/150‑7 (MQ=255)
                          cgcTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGCGCCTGCTTGCTCAGGGTGTCGATGTGCTgg                                                                                                                            <  1:2277212/150‑1 (MQ=255)
                          cgcTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTgg                                                                                                                            <  1:1657593/150‑1 (MQ=255)
                          cgcTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGATCAGGGTGTCGTTGTGCTgg                                                                                                                            <  2:2368086/150‑1 (MQ=255)
                               ttCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAg                                                                                                                       >  2:1815982/1‑150 (MQ=255)
                                       gCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCCCCCGGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCAc                                                                                                               >  1:723382/1‑150 (MQ=255)
                                             cgggTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGtt                                                                                                         <  2:930370/150‑1 (MQ=255)
                                                gTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTtgct                                                                                                      >  1:508544/1‑150 (MQ=255)
                                                gTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTtgct                                                                                                      >  1:1273384/1‑150 (MQ=255)
                                                  tcgccTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTtgctgc                                                                                                    <  2:198883/147‑1 (MQ=255)
                                                    aCCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTtgctgc                                                                                                    >  2:575687/1‑148 (MQ=255)
                                                    aCCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAg                                                                                                  >  1:730027/1‑150 (MQ=255)
                                                      cTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgt                                                                                                >  2:943799/1‑150 (MQ=255)
                                                       ttCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgtc                                                                                               >  2:1479002/1‑150 (MQ=255)
                                                       ttCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgtc                                                                                               >  1:1933157/1‑150 (MQ=255)
                                                       ttCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgtc                                                                                               >  1:1905677/1‑150 (MQ=255)
                                                       ttCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgtc                                                                                               >  2:535184/1‑150 (MQ=255)
                                                       ttCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgtc                                                                                               >  1:39966/1‑150 (MQ=255)
                                                        tCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgtcg                                                                                              >  1:408420/1‑150 (MQ=255)
                                                           aaaTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAGGgtg                                                                                                                                        >  2:1276010/1‑105 (MQ=255)
                                                           aaaTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAGGgtg                                                                                                                                        <  1:1276010/105‑1 (MQ=255)
                                                            aaTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGa                                                                                          <  1:1975773/150‑1 (MQ=255)
                                                            aaTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGa                                                                                          <  1:451633/150‑1 (MQ=255)
                                                                gATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCCGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCaa                                                                                      >  1:577640/1‑150 (MQ=255)
                                                                gATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCaa                                                                                      <  1:1648203/150‑1 (MQ=255)
                                                                gATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCAACCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCaa                                                                                      >  1:755837/1‑150 (MQ=255)
                                                                     ggTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACt                                                                                 <  1:1210378/150‑1 (MQ=255)
                                                                     ggTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACt                                                                                 <  1:982358/150‑1 (MQ=255)
                                                                      gTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTg                                                                                >  1:1352357/1‑150 (MQ=255)
                                                                      gTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCATGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTg                                                                                >  2:2440514/1‑150 (MQ=255)
                                                                        aGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGc                                                                              <  1:1815982/150‑1 (MQ=255)
                                                                         ggCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCt                                                                             >  1:609439/1‑150 (MQ=255)
                                                                              cACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGc                                                                        >  2:1923986/1‑150 (MQ=255)
                                                                               aCCCTGCACCAAGCGCACCTGGACCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCt                                                                       <  2:723382/150‑1 (MQ=255)
                                                                                cccTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAg                                                                                                  <  2:473669/122‑1 (MQ=255)
                                                                                 ccTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTa                                                                            >  1:1129713/1‑143 (MQ=255)
                                                                                 ccTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTa                                                                            <  2:1129713/143‑1 (MQ=255)
                                                                                 ccTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGt                                                                     >  2:498166/1‑150 (MQ=255)
                                                                                   tGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGt                                                                                                                  <  1:1493857/103‑1 (MQ=255)
                                                                                   tGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGt                                                                                                                  >  2:1493857/1‑103 (MQ=255)
                                                                                   tGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTcc                                                                   >  2:598242/1‑150 (MQ=255)
                                                                                    gCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgtc                                                                                               <  1:496639/121‑1 (MQ=255)
                                                                                        cacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgtgt                                                              >  2:730603/1‑150 (MQ=255)
                                                                                           ctgNACCTGGGCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTAGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTc                                                           >  2:1497954/3‑150 (MQ=255)
                                                                                           ctgCACCTGGGCCAGGTGTAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgtcgtc                                                                                            <  1:1448528/115‑1 (MQ=255)
                                                                                                 cTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTggg                                                                                                                           <  1:2129650/80‑1 (MQ=255)
                                                                                                 cTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTggg                                                                                                                           >  2:2129650/1‑80 (MQ=255)
                                                                                                    gCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTg                                                  <  1:1649433/150‑1 (MQ=255)
                                                                                                      ccAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAg                                                                                                                                            <  2:887497/58‑1 (MQ=255)
                                                                                                      ccAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAg                                                                                                                                            >  1:887497/1‑58 (MQ=255)
                                                                                                      ccAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTggg                                                >  2:770647/1‑150 (MQ=255)
                                                                                                      ccAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTg                                                  >  1:276289/1‑148 (MQ=255)
                                                                                                      ccAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTg                                                  <  2:276289/148‑1 (MQ=255)
                                                                                                          gTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTgg                                                 >  1:2338602/1‑145 (MQ=255)
                                                                                                          gTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTgg                                                 <  2:2338602/145‑1 (MQ=255)
                                                                                                           tGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCcagc                                           >  2:1554863/1‑150 (MQ=255)
                                                                                                            gCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCcagca                                          <  1:598242/150‑1 (MQ=255)
                                                                                                                  cGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGa                                    <  2:287814/150‑1 (MQ=255)
                                                                                                                  cGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGa                                    <  2:408420/150‑1 (MQ=255)
                                                                                                                  cGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGa                                    <  1:770647/150‑1 (MQ=255)
                                                                                                                   gCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAg                                   <  2:1352357/150‑1 (MQ=255)
                                                                                                                   gCCGTGTACACGCATAAGCAGGCCACCCTGAGCATGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAg                                   <  1:2440514/150‑1 (MQ=255)
                                                                                                                      gtgtACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGc                                >  1:1535816/1‑150 (MQ=255)
                                                                                                                       tgtACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCg                               >  2:2199923/1‑150 (MQ=255)
                                                                                                                         tACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTg                             >  2:409689/1‑150 (MQ=255)
                                                                                                                          tctcGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGACATAAGGATTTATCAAGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTgg                            <  2:1880439/147‑1 (MQ=255)
                                                                                                                          acacGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTgg                            <  2:577640/150‑1 (MQ=255)
                                                                                                                              gCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgt                                                                <  1:1409476/110‑1 (MQ=255)
                                                                                                                              gCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtca                        >  1:1402389/1‑150 (MQ=255)
                                                                                                                              gCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtca                        >  2:1046787/1‑150 (MQ=255)
                                                                                                                              gCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtca                        >  2:2404333/1‑150 (MQ=255)
                                                                                                                              gCATAAGCAGGCCACCCTGAGCCTGCTTGCGCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgt                                                                >  2:1409476/1‑110 (MQ=255)
                                                                                                                                aTAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGc                                                    <  1:1237689/120‑1 (MQ=255)
                                                                                                                                aTAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGc                                                    >  2:1237689/1‑120 (MQ=255)
                                                                                                                                 tAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAt                                                                                                                                   >  2:81538/1‑40 (MQ=255)
                                                                                                                                 tAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTgg                                                                               >  2:2376142/1‑92 (MQ=255)
                                                                                                                                 tAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTgg                                                                               <  1:2376142/92‑1 (MQ=255)
                                                                                                                                    gCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTgg                            <  1:744747/140‑1 (MQ=255)
                                                                                                                                     cAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcagttcagt                 <  1:987559/150‑1 (MQ=255)
                                                                                                                                      aGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcagttccgac                <  1:1046787/150‑5 (MQ=255)
                                                                                                                                      aGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTg                <  1:498166/150‑1 (MQ=255)
                                                                                                                                        gacACCCTGAGCAAGCAGGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGt                              <  2:1728866/132‑1 (MQ=255)
                                                                                                                                        gCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCt              <  1:730603/150‑1 (MQ=255)
                                                                                                                                           aCCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGACAGTTCAGTGCTGAt           >  2:1830227/1‑150 (MQ=255)
                                                                                                                                             ccTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGTTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCt         >  2:1062219/1‑150 (MQ=255)
                                                                                                                                             ccTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGGGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCt         >  2:152934/1‑150 (MQ=255)
                                                                                                                                               cgAGCCTGCATGCACAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCtgt       <  2:508544/149‑1 (MQ=255)
                                                                                                                                                    cTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGcaca  >  1:2236149/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCACA  >  NC_002947/4580495‑4580792

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