Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,931,240 A→G 10.5% S59P (TCC→CCC)  cheY ← two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,931,2400AG10.5% 162.7 / 7.6 65S59P (TCC→CCC) cheYtwo‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base G (4/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

AACGGCTTGACCACATAACCATTGACGCCAGCCTGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCAACAGGTTCTTGATGATCCGCCGCATCGTCGAA  >  NC_002947/4931098‑4931388
                                                                                                                                              |                                                                                                                                                    
aaCGGCTTGACCACATAACCATTGACGCCAGCCTGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGccc                                                                                                                                               <  1:1405899/150‑1 (MQ=255)
aaCGACTTGACCACATAACCATTCACGCCAGCCTGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGccc                                                                                                                                               <  1:541353/150‑1 (MQ=255)
     cTTGACCACATAACCATTGACGCCAGCCTGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCAt                                                                                                                                          >  2:1022666/1‑150 (MQ=255)
            acaTAACCATTGACGCCAGCCTGGGCCGCTTCAATGAACTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAg                                                                                                                                   >  1:1153172/1‑150 (MQ=255)
              aTAACCATTGACGCCAGCCTGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTc                                                                                                                                 <  2:1187367/150‑1 (MQ=255)
               tAACCATTGACGCCAGCCTGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCg                                                                                                                                <  1:342540/150‑1 (MQ=255)
                aaCCATTGACGCCAGCCTGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCgg                                                                                                                               <  1:264387/150‑1 (MQ=255)
                   cATTGACGCCAGCCTGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCa                                                                                                                            >  1:507829/1‑150 (MQ=255)
                          gCCAGCCTGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGaa                                                                                                                     <  2:416818/150‑1 (MQ=255)
                            cagacTGTGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGGCATGTCCGGCAt                                                                                                                                          <  1:572345/123‑1 (MQ=255)
                            cAGCCTGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGt                                                                                                                   >  1:539963/1‑150 (MQ=255)
                            cAGACTGTGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGGCATGTCCGGCAt                                                                                                                                          >  2:572345/1‑127 (MQ=255)
                            cAGACTGTGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGGCATGTCCGGCAt                                                                                                                                          <  1:572114/127‑1 (MQ=255)
                            cAGACTGTGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGGCATGTCCGGCAt                                                                                                                                          >  2:572114/1‑127 (MQ=255)
                             aGCCTGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTc                                                                                                                  >  2:736990/1‑150 (MQ=255)
                               ccTGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGt                                                                                                                >  1:961221/1‑150 (MQ=255)
                                 tGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAg                                                                                                              <  2:795104/150‑1 (MQ=255)
                                            aaTGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTc                                                                                                   <  2:1154866/150‑1 (MQ=255)
                                            aaTGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACATAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTc                                                                                                   >  1:1095453/1‑150 (MQ=255)
                                                     gTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGCCATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCgg                                                                                          >  2:588544/1‑150 (MQ=255)
                                                     gTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCgg                                                                                          >  1:828576/1‑150 (MQ=255)
                                                      tCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGAATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGc                                                                                         >  1:86225/1‑150 (MQ=255)
                                                         cgcTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAAc                                                                                      >  2:366517/1‑150 (MQ=255)
                                                            ttGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGt                                                                                                                              >  1:194697/1‑107 (MQ=255)
                                                            ttGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGt                                                                                                                              <  2:194697/107‑1 (MQ=255)
                                                              ggCTTCAGCGGTCACCATCCAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCgt                                                                                 <  2:34515/150‑1 (MQ=255)
                                                                cTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGGCATGTCCGGCAtcgac                                                                                                                                      >  1:1252175/1‑91 (MQ=255)
                                                                cTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGGCATGTCCGGCAtcgac                                                                                                                                      <  2:1252175/95‑5 (MQ=255)
                                                                cTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCgtgg                                                                               >  1:265259/1‑150 (MQ=255)
                                                                 ttCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCgtggt                                                                              >  2:497916/1‑150 (MQ=255)
                                                                  tCAGCGGTCACCATCAAGACCGGCATGCTTTTTAGCTTGTGGCTGAGACGCACTTTACGCAGCATGTCGATGCCGGACTTACCCGGCCTGTTCCAGTCGGTCAGCAGGAATTCGTAGTGACCGTGGTCCAGCATACCCAACCCCgtggtg                                                                             >  2:537571/1‑150 (MQ=255)
                                                                        gTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCa                                                                       >  2:355480/1‑150 (MQ=255)
                                                                             cATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGc                                                                                         >  1:833636/1‑127 (MQ=255)
                                                                             cATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGc                                                                                         <  2:833636/127‑1 (MQ=255)
                                                                               tCAAGATCGGCATGCTTTTCAGCTTGTTGCTGGCACGCACTTTACGCAGCAGGTCGATGCCCGACATGCCCGGCATGTTCCAGTCGGTCACCTGTAAGTTGTATTGACCGTTTTCCAGCATCTGCAAATCCGTAGTTCCGTCATCgtttt                                                                <  2:479097/150‑5 (MQ=255)
                                                                                 aaGACCGGCATGCTTTTCAGCTTTGCGCTTGCGCGCACTTTACGCAGCAGGTCGATGCCGGACATGCGCTTCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTtcg                                                              <  1:1022666/150‑1 (MQ=255)
                                                                                  aGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTtcgt                                                             >  2:510536/1‑150 (MQ=255)
                                                                                    cccGGCACGCTTTTCAGCGTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTtcgtcg                                                           <  1:1035453/149‑1 (MQ=255)
                                                                                    aCCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTtcgtcg                                                           <  1:497916/150‑1 (MQ=255)
                                                                                         cATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGtt                                                      <  2:1153172/150‑1 (MQ=255)
                                                                                          aTGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTg                                                     <  1:355480/150‑1 (MQ=255)
                                                                                          aTGCTTTTCAGCTTGTCGCTGCAAAGTAGTTTTCGCAGCAGTTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTg                                                     <  1:91/150‑1 (MQ=255)
                                                                                           tGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTtcgtcg                                                           >  1:1130842/1‑143 (MQ=255)
                                                                                           tGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTtcgtcg                                                           <  2:1130842/143‑1 (MQ=255)
                                                                                            gCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGt                                                   <  1:755066/150‑1 (MQ=255)
                                                                                              ttttCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAg                                                                                                              >  2:988971/1‑89 (MQ=255)
                                                                                              ttttCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAg                                                                                                              <  1:988971/89‑1 (MQ=255)
                                                                                                   aGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGACGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCTGCTTCTTCGTTGTTTGGGAAAccc                                            <  2:797356/150‑1 (MQ=255)
                                                                                                      ttGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTGCCAGTCGGTCACCAGGAAGACGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAgg                                         >  2:1142485/1‑150 (MQ=255)
                                                                                                          cGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCAc                                     <  2:1017598/150‑1 (MQ=255)
                                                                                                             tGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCa                                  >  2:616248/1‑150 (MQ=255)
                                                                                                             tGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCa                                  >  2:182275/1‑150 (MQ=255)
                                                                                                                 acgcacTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAg                                                                                                              >  2:211652/1‑70 (MQ=255)
                                                                                                                 acgcacTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAg                                                                                                              <  1:211652/70‑1 (MQ=255)
                                                                                                                     acTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCAAc                                <  2:658711/144‑1 (MQ=255)
                                                                                                                     acTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCAAc                                >  1:658711/1‑144 (MQ=255)
                                                                                                                       tttACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCAACAGGTTCtt                        <  2:539963/150‑1 (MQ=255)
                                                                                                                         tACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCAACAGGTTCTtga                      <  2:828576/150‑1 (MQ=255)
                                                                                                                              agcagGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCAACAGGTTCTTGATGATc                 >  1:1287271/1‑150 (MQ=255)
                                                                                                                              agcagGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAAGCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCAACAGGTTCTTGATGATc                 >  1:36784/1‑150 (MQ=255)
                                                                                                                                  ggTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCAACAGGTTCTTGATGATccgcc             <  2:507829/150‑1 (MQ=255)
                                                                                                                                      gATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCAACAGGTTCTTGATGATCCGCCGCAt         >  1:1137600/1‑150 (MQ=255)
                                                                                                                                      caTGCCGGGCATGTCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCAACGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCAACAGGAGTTTGATGATCCGCCGCAt         >  1:1141676/2‑150 (MQ=255)
                                                                                                                                       aTGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGGAACCCAGGTCACGCAACAGGTTCTTGATGATCCGCCGCAtc        >  1:1084937/1‑150 (MQ=255)
                                                                                                                                          ccGGACCTGCCCGCCATCCGCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCAACAGGTTCTTGATGATCCGCCGCAtcgtc     <  2:36784/150‑1 (MQ=255)
                                                                                                                                            gaaCATGCCGGGCATGTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCAACAGGTTCTTGATGATCCGCCGCATCGTCGaa  >  2:938593/3‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
AACGGCTTGACCACATAACCATTGACGCCAGCCTGGGCCGCTTCAATGATCTGGTCGCGCTTGGCTTCAGCGGTCACCATCAAGACCGGCATGCTTTTCAGCTTGTCGCTGGCACGCACTTTACGCAGCAGGTCGATGCCGGACATGCCCGGCATGTTCCAGTCGGTCACCAGGAAGTCGTAGTGACCGTTTTCCAGCATCGGCAACGCCGTGGTGCCGTCATCGGCTTCGTCGGTGTTGGTGAAACCCAGGTCACGCAACAGGTTCTTGATGATCCGCCGCATCGTCGAA  >  NC_002947/4931098‑4931388

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