Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,427,963 A→T 6.0% L63Q (CTG→CAG)  rhlE‑I ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,427,9630AT 6.0% 280.1 / 5.2 91L63Q (CTG→CAG) rhlE‑IATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (46/39);  new base T (2/4);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCCGCCATACACCGCGTAAGTAGTCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGCTG  >  NC_002947/5427816‑5428102
                                                                                                                                                   |                                                                                                                                           
aCCCCGCCATACACCGCGTAATTAGTCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCCCAACCTATCCTCGGATCGATTTTCTGGCCACCTTCTCGCCCTCCAGCGCCAgg                                                                                                                                           <  1:824529/150‑1 (MQ=255)
aCCCCGCCATACACCGCGTAAGTAGTCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGTGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAgg                                                                                                                                           >  2:256598/1‑150 (MQ=255)
aCCCCGCCATACACCGCGTAAGTAGTCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAgg                                                                                                                                           >  2:1490507/1‑150 (MQ=255)
  cccGCCATACACCGCGTAAGTAGTCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGgcg                                                                                                                                         >  2:1940860/1‑150 (MQ=255)
          acacCGCGTAAGTAGTCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCTGGCGCTggaggg                                                                                                                                 <  1:1644026/150‑6 (MQ=255)
           cacCGCGTAAGTAGTCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGAGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGcac                                                                                                                                >  2:21380/1‑150 (MQ=255)
            acCGCGTAAGTAGTCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGcaca                                                                                                                               >  2:1905057/1‑150 (MQ=255)
             ccGCGTAAGTAGTCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCAcag                                                                                                                              >  2:516468/1‑150 (MQ=255)
             ccGCGTAAGTAGTCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCAcag                                                                                                                              <  1:934809/150‑1 (MQ=255)
                  tAAGTAGTCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCAcaggcagg                                                                                                                         >  2:1756816/1‑150 (MQ=255)
                    agtagtCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGc                                                                                                                       >  1:1891036/1‑150 (MQ=255)
                      tagtCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTgc                                                                                                                                     >  2:895151/1‑134 (MQ=255)
                      tagtCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTgc                                                                                                                                     <  1:895151/134‑1 (MQ=255)
                         tCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGTATTGCGCTGCAGCACAGGCAGGGCAAAGc                                                                                                                  <  1:1321864/150‑1 (MQ=255)
                         tCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGc                                                                                                                  <  1:1065317/150‑1 (MQ=255)
                                  ggTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGGTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCt                                                                                                         >  2:1857837/1‑150 (MQ=255)
                                   gTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCtt                                                                                                        >  1:1516229/1‑150 (MQ=255)
                                     tctcGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGc                                                                                                      >  1:696756/1‑150 (MQ=255)
                                      ctcGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGcc                                                                                                     <  2:1891036/150‑1 (MQ=255)
                                        cGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCgg                                                                                                   <  2:645206/150‑1 (MQ=255)
                                         ggCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCg                                                                                                    >  1:411978/1‑148 (MQ=255)
                                         ggCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCg                                                                                                    <  2:411978/148‑1 (MQ=255)
                                          gCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTa                                                                                                 <  1:1940860/150‑1 (MQ=255)
                                                  gcgcACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCAAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCAcaggcag                                                                                                                          >  2:632232/1‑117 (MQ=255)
                                                  gcgcACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCAAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCAcaggcag                                                                                                                          <  1:632232/117‑1 (MQ=255)
                                                    gcACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGc                                                                                                              <  2:1906556/127‑1 (MQ=255)
                                                    gcACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGc                                                                                                              >  1:1906556/1‑127 (MQ=255)
                                                        gttgttGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCTGGCGCTggag                                                                                                                                   <  2:1106865/102‑1 (MQ=255)
                                                        gttgttGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGc                                                                                   >  2:362521/1‑150 (MQ=255)
                                                        gttgttGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCCCCAGCGCCTGGCGCTggag                                                                                                                                   >  1:1106865/1‑99 (MQ=255)
                                                                gCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCg                                                                                                    >  2:1702229/1‑125 (MQ=255)
                                                                gCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCg                                                                                                    <  1:1702229/125‑1 (MQ=255)
                                                                gCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCa                                                                           <  1:1490507/150‑1 (MQ=255)
                                                                 cACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAg                                                                          >  1:608161/1‑150 (MQ=255)
                                                                   ccTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGt                                                                        >  1:537009/1‑150 (MQ=255)
                                                                    cTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGc                                                                                                      <  1:1263364/119‑1 (MQ=255)
                                                                    cTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGc                                                                                                      >  2:1263364/1‑119 (MQ=255)
                                                                    cTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTc                                                                       >  1:156668/1‑150 (MQ=255)
                                                                        tCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGcgg                                                                   >  1:962907/1‑150 (MQ=255)
                                                                         cGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGcggc                                                                  <  2:1930774/150‑1 (MQ=255)
                                                                            ccAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGACAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGcggccgg                                                               <  1:1871473/150‑1 (MQ=255)
                                                                             cAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGcggccggc                                                              <  2:156668/150‑1 (MQ=255)
                                                                             cAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGcggccggc                                                              >  2:841954/1‑150 (MQ=255)
                                                                             cAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGcggccggc                                                              >  1:880879/1‑150 (MQ=255)
                                                                                 tCGCGAGTGAGTACCAGCACCAAGGCGCGGATCGAGTGGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGc                                                          >  1:1568168/1‑150 (MQ=255)
                                                                                  cgcgAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGcggc                                                                  <  2:1148481/141‑1 (MQ=255)
                                                                                  cgcgAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGcggc                                                                  >  1:1148481/1‑141 (MQ=255)
                                                                                   gcgAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTg                                                                                                       >  1:1260870/1‑103 (MQ=255)
                                                                                   gcgAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCCTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTg                                                                                                       <  2:1260870/103‑1 (MQ=255)
                                                                                       gTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCAcggcg                                                    <  2:537009/150‑1 (MQ=255)
                                                                                        tGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCAcggcgg                                                   <  2:822503/150‑1 (MQ=255)
                                                                                           gTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGa                                                >  1:185364/1‑150 (MQ=255)
                                                                                           gTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGa                                                >  2:1817787/1‑150 (MQ=255)
                                                                                               cAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGc                                            >  2:1484883/1‑150 (MQ=255)
                                                                                                   aCCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCATCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTg                                        <  1:21380/150‑1 (MQ=255)
                                                                                                   aCCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTg                                        <  1:516468/150‑1 (MQ=255)
                                                                                                   aCCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTg                                        <  1:256598/150‑1 (MQ=255)
                                                                                                   aCCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTg                                        >  2:1693260/1‑150 (MQ=255)
                                                                                                    ccAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATc                                                                            <  1:1344960/113‑1 (MQ=255)
                                                                                                    ccAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATc                                                                            >  2:1344960/1‑113 (MQ=255)
                                                                                                      aaGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGcc                                     <  1:1905057/150‑1 (MQ=255)
                                                                                                      aaGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGcc                                     <  1:1023253/150‑1 (MQ=255)
                                                                                                       aGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCt                                    >  1:160576/1‑150 (MQ=255)
                                                                                                              gATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACcgg                             >  2:612532/1‑150 (MQ=255)
                                                                                                               aTCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCAt                                                                             >  2:1808581/1‑101 (MQ=255)
                                                                                                               aTCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCAt                                                                             <  1:1808581/101‑1 (MQ=255)
                                                                                                                tCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGa                                                <  2:1772187/129‑1 (MQ=255)
                                                                                                                tCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGa                                                >  1:1772187/1‑129 (MQ=255)
                                                                                                                   aGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACcggggtcg                        >  1:1291192/1‑150 (MQ=255)
                                                                                                                       gCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACcggggtcggggt                    <  2:880879/150‑1 (MQ=255)
                                                                                                                        cTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTg                                                                                         <  2:1748171/80‑1 (MQ=255)
                                                                                                                        cTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTg                                                                                         >  1:1748171/1‑80 (MQ=255)
                                                                                                                         tGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTgg                                                                                        >  1:71840/1‑80 (MQ=255)
                                                                                                                         tGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTgg                                                                                        <  2:71840/80‑1 (MQ=255)
                                                                                                                          ggCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGggtggt                 <  2:160576/150‑1 (MQ=255)
                                                                                                                           gCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGgtggtg                >  2:387184/1‑150 (MQ=255)
                                                                                                                           gCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGgtggtg                <  2:962907/150‑1 (MQ=255)
                                                                                                                           gCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGgtggtg                <  1:480264/150‑1 (MQ=255)
                                                                                                                             cACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTa              >  2:1563490/1‑150 (MQ=255)
                                                                                                                                cTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTc           <  1:1202803/150‑1 (MQ=255)
                                                                                                                                 ttCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTcc          >  2:1448542/1‑150 (MQ=255)
                                                                                                                                    tcGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGACATCAGGTCGAGGCCGCCCAGCACGGCGGTGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGt       <  1:1693260/150‑1 (MQ=255)
                                                                                                                                    tcGCCCTCCAGCGCCAGGCGCTGAAGAACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGt       <  1:1857837/150‑1 (MQ=255)
                                                                                                                                      gCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTGTTGCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTtg     >  2:1910544/1‑150 (MQ=255)
                                                                                                                                       cccTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTtgc    >  2:1671858/1‑150 (MQ=255)
                                                                                                                                        ccTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGggtgg                  <  1:1307372/135‑1 (MQ=255)
                                                                                                                                        ccTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGggtgg                  >  2:1307372/1‑135 (MQ=255)
                                                                                                                                         ctgcAGCGCCTGGCGCTGGAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGCGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGAAGTCCAGTtgctg  >  1:398609/4‑150 (MQ=255)
                                                                                                                                         cTCCAGCGACTGGGCCTTGAGCTCAGGCAGGGCAAAGCCTGCGGTCTTGCCAGTACGGCTGTTGCGAACCGGGTTCAGGTCGCGGCCGTTCACCACGGCGGGGATGGCTTTGGGTTGCGCCGGGGTCGGGGTGGTGTAGTCCAGTtgctg  <  1:1893600/150‑1 (MQ=255)
                                                                                                                                                gcCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACcg                              <  2:1558731/115‑1 (MQ=255)
                                                                                                                                                gcCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACcg                              >  1:1558731/1‑115 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
ACCCCGCCATACACCGCGTAAGTAGTCAGCGGCAGGTTCTCGGCGTATTCGCGCACGTTGTTGTGCACCTGCTCGGCCAGCTCGCGAGTGGGTACCAGCACCAAGGCGCGGATCGAGTTGCTGGCCACCTTCTCGCCCTCCAGCGCCAGGCGCTGCAGCACAGGCAGGGCAAAGCCTGCGGTCTTGCCGGTACCGGTCTGGGCTGCGGCCATCAGGTCGCGGCCGGCCAGCACGGCGGGGATGGCTTTGGCCTGCACCGGGGTCGGGGTGGTGTAGTCCAGTTGCTG  >  NC_002947/5427816‑5428102

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