Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,519,503 Δ1 bp 5.4% coding (61/2271 nt) clpA ← ATP‑dependent serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,519,5030G. 5.4% 235.6 / 6.7 74coding (61/2271 nt)clpAATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base G (33/37);  new base . (2/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

AGGGATCAACGGCGTGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGCGTATTGGT  >  NC_002947/4519354‑4519652
                                                                                                                                                     |                                                                                                                                                     
aGGGATCAACGGCGTGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCATGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACg                                                                                                                                                       >  2:823286/1‑150 (MQ=255)
aGGGATCAACGGCGGGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCATGGAGTGCCAGCAGCCGATGCTCGACGGTCATGAACTCATGACg                                                                                                                                                       >  2:822401/1‑150 (MQ=255)
       aaCGGCGTGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGCTTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGaa                                                                                                                                                >  1:888289/1‑150 (MQ=255)
       aaCGGCGTGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGaa                                                                                                                                                >  2:493549/1‑150 (MQ=255)
                   gAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTg                                                                                                                                    <  1:1372010/150‑1 (MQ=255)
                   gAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGAGTGGCGCCACCGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGCCGCTTCGAACGGGCCTCCTTg                                                                                                                                    >  2:398956/1‑150 (MQ=255)
                     gTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGaa                                                                                                                                  <  2:1283299/150‑1 (MQ=255)
                     gTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGaa                                                                                                                                  <  2:1440648/150‑1 (MQ=255)
                     gTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGaa                                                                                                                                  >  1:296594/1‑150 (MQ=255)
                     gTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCACCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAAATCACTGCGCTTCGAACGGGCCTCCTTGaa                                                                                                                                  <  1:324202/150‑1 (MQ=255)
                      tCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAg                                                                                                                                 <  1:899313/150‑1 (MQ=255)
                      tCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAg                                                                                                                                 <  2:1086838/150‑1 (MQ=255)
                      gcGATGAACTCTTCCAGGTCATGCTTGCGCTTGTCGAGTCTGTGAAGTGATTTGCGTAGAACTGTCGCGGCAGCCTCATGGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAg                                                                                                                                 <  2:1843045/149‑1 (MQ=255)
                        gATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGc                                                                                                                               >  1:1788371/1‑150 (MQ=255)
                        gATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGc                                                                                                                               >  1:2030418/1‑150 (MQ=255)
                        gATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGGGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGc                                                                                                                               >  2:1038540/1‑150 (MQ=255)
                          tGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCaa                                                                                                                             <  2:1738797/150‑1 (MQ=255)
                               tctTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGTCaccctca                                                                                                                        >  2:2131110/1‑144 (MQ=255)
                               tctTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTg                                                                                                                        >  1:1445554/1‑150 (MQ=255)
                                 ttGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAg                                                                                                                      <  1:493549/150‑1 (MQ=255)
                                 ttGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTACCAGCAGCAGATGCTCGCCGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAg                                                                                                                      <  1:53927/150‑1 (MQ=255)
                                      ggTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCGTGAAGGCAAGATTGAGGGTGa                                                                                                                 <  2:619138/150‑1 (MQ=255)
                                       gTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGAc                                                                                                                >  2:2394474/1‑150 (MQ=255)
                                          aTGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGAtt                                                                                                                         <  2:1792194/138‑1 (MQ=255)
                                          aTGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGAtt                                                                                                                         >  1:1792194/1‑138 (MQ=255)
                                            gCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGa                                                                                                           >  1:1018599/1‑150 (MQ=255)
                                             cTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAg                                                                                                          >  1:372005/1‑150 (MQ=255)
                                                   cTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGCCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTcgcg                                                                                                    <  2:678556/150‑1 (MQ=255)
                                                    ttGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCgg                                                                                                   >  2:2076167/1‑150 (MQ=255)
                                                      gTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGtt                                                                                                 >  1:1660658/1‑150 (MQ=255)
                                                       tCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGt                                                                                                                   >  1:1985956/1‑131 (MQ=255)
                                                       tCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGt                                                                                                                   <  2:1985956/131‑1 (MQ=255)
                                                         gagaTTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTaa                                                                                              >  2:1161164/1‑150 (MQ=255)
                                                              tGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAg                                                                                         <  2:1788371/150‑1 (MQ=255)
                                                                    cacaGGCGCGCAGAACGGGCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCAc                                                                                   <  2:573339/150‑1 (MQ=255)
                                                                      caGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACggg                                                                                                                                            >  1:282025/1‑91 (MQ=255)
                                                                      caGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACggg                                                                                                                                            <  2:282025/91‑1 (MQ=255)
                                                                       aGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTc                                                                                >  2:384993/1‑150 (MQ=255)
                                                                         gcgcgcAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCAt                                                                              >  1:75793/1‑150 (MQ=255)
                                                                             gcAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCGTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATAccc                                                                          >  2:1001806/1‑150 (MQ=255)
                                                                                      tCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCg                                                                 >  1:1021833/1‑150 (MQ=255)
                                                                                          ggCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGa                                                             <  2:888289/150‑1 (MQ=255)
                                                                                              gCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTaa                                                         <  2:1018599/150‑1 (MQ=255)
                                                                                              gCCTCATTGTCAAGGAGTGCCAGCAGCAGATACTACACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGCCGGCGATTaa                                                         <  1:1440738/150‑1 (MQ=255)
                                                                                               ccTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCtt                                                                                      <  1:2152275/120‑1 (MQ=255)
                                                                                               ccTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCtt                                                                                      >  2:2152275/1‑120 (MQ=255)
                                                                                               ccTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGa                                                             >  1:2228669/1‑145 (MQ=255)
                                                                                               ccTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGa                                                             <  2:2228669/145‑1 (MQ=255)
                                                                                                 tCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCg                                                      >  1:868709/1‑150 (MQ=255)
                                                                                                     tGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCt                                                  <  2:2030418/150‑1 (MQ=255)
                                                                                                         aaGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTtctc                                              <  1:516895/150‑1 (MQ=255)
                                                                                                                         agATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGc                              <  1:1065290/150‑1 (MQ=255)
                                                                                                                         agATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGc                              <  1:823286/150‑1 (MQ=255)
                                                                                                                         agATGCGCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGc                              <  1:822401/150‑1 (MQ=255)
                                                                                                                          gATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCg                             >  2:100219/1‑150 (MQ=255)
                                                                                                                               tCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGc                        >  2:281525/1‑150 (MQ=255)
                                                                                                                                gaagcgTCATGAGTTCATGACCGTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTg                      >  2:2374186/6‑150 (MQ=255)
                                                                                                                                gaagcgTCATGAGTTCATGACCGTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTg                      >  2:2374343/6‑150 (MQ=255)
                                                                                                                                gaagcgTCATGAGTTCATGACCGTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTg                      <  1:2374343/145‑1 (MQ=255)
                                                                                                                                gaagcgTCATGAGTTCATGACCGTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTg                      <  1:2374186/145‑1 (MQ=255)
                                                                                                                                   gcgTCATGAGTTCATGACGCTTCGAACGGCCCTCCTTGAAGGCAAGATTGAGGGTTACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGGTGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGc                    >  2:909995/3‑150 (MQ=255)
                                                                                                                                    cgTCATGAGTTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCt                   >  2:2392419/2‑150 (MQ=255)
                                                                                                                                      tCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTc                                                                  >  2:1653922/1‑101 (MQ=255)
                                                                                                                                      tCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTc                                                                  <  1:1653922/101‑1 (MQ=255)
                                                                                                                                          gAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCt             >  2:2155684/1‑150 (MQ=255)
                                                                                                                                            aCTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTgg           <  1:1038540/150‑1 (MQ=255)
                                                                                                                                             cTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGc          <  1:281525/150‑1 (MQ=255)
                                                                                                                                             cTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGc          <  1:531532/150‑1 (MQ=255)
                                                                                                                                             cTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGc          <  2:1531483/150‑1 (MQ=255)
                                                                                                                                             cTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGc          <  1:1161164/150‑1 (MQ=255)
                                                                                                                                             cTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACAGCCGGGCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGc          <  1:1001806/150‑1 (MQ=255)
                                                                                                                                              tCATGACGCTTCGAACCGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGAATCACCGTCCTACTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGCg         >  2:1970201/1‑150 (MQ=255)
                                                                                                                                                aTGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGCGTa       <  1:384993/150‑1 (MQ=255)
                                                                                                                                                  gACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCCTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGCGTAtt     <  2:868709/150‑1 (MQ=255)
                                                                                                                                                   aCGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACCGAGTAGCGGATGCTGGCTTTCCCTGGCGTATTg    <  1:2076167/150‑1 (MQ=255)
                                                                                                                                                   aCGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGCGTATTg    >  2:299421/1‑150 (MQ=255)
                                                                                                                                                     gCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACaga                                   >  1:2393295/1‑117 (MQ=255)
                                                                                                                                                     gCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACaga                                   <  2:2393295/117‑1 (MQ=255)
                                                                                                                                                     gCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGCGTATTGGt  <  2:1021833/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGGATCAACGGCGTGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGCGTATTGGT  >  NC_002947/4519354‑4519652

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: