Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,519,495 C→T 8.7% E23E (GAG→GAA)  clpA ← ATP‑dependent serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,519,4950CT 8.7% 186.4 / 6.7 69E23E (GAG→GAA) clpAATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base C (33/29);  new base T (4/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

ACAGGGATCAACGGCGTGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGC  >  NC_002947/4519352‑4519644
                                                                                                                                               |                                                                                                                                                     
aCAGGGATCAACGGCGTGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGa                                                                                                                                                 >  1:1738797/1‑150 (MQ=255)
  aGGGATCAACGGCGTGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCATGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACg                                                                                                                                               >  2:823286/1‑150 (MQ=255)
  aGGGATCAACGGCGTGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCa                                                                                                                                                    >  2:1448850/1‑145 (MQ=255)
  aGGGATCAACGGCGTGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCa                                                                                                                                                    <  1:1448850/145‑1 (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                                                                                                                          >  1:296594/1‑150 (MQ=255)
                       gTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGaa                                                                                                                          <  2:1440648/150‑1 (MQ=255)
                       gTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGaa                                                                                                                          <  2:1283299/150‑1 (MQ=255)
                       gTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCACCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAAATCACTGCGCTTCGAACGGGCCTCCTTGaa                                                                                                                          <  1:324202/150‑1 (MQ=255)
                        tCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAg                                                                                                                         <  2:1086838/150‑1 (MQ=255)
                        tCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAg                                                                                                                         <  1:899313/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                                                                                                                 >  1:1792194/1‑138 (MQ=255)
                                            aTGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGAtt                                                                                                                 <  2:1792194/138‑1 (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                                                                                                           <  2:1985956/131‑1 (MQ=255)
                                                         tCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGt                                                                                                           >  1:1985956/1‑131 (MQ=255)
                                                           gagaTTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTaa                                                                                      >  2:1161164/1‑150 (MQ=255)
                                                                tGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAg                                                                                 <  2:1788371/150‑1 (MQ=255)
                                                                      cacaGGCGCGCAGAACGGGCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCAc                                                                           <  2:573339/150‑1 (MQ=255)
                                                                        caGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACggg                                                                                                                                    <  2:282025/91‑1 (MQ=255)
                                                                        caGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACggg                                                                                                                                    >  1:282025/1‑91 (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:823286/150‑1 (MQ=255)
                                                                                                                           agATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGc                      <  1:1065290/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)
                                                                                                                                  gaagcgTCATGAGTTCATGAC‑CGTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTg              >  2:2374186/6‑150 (MQ=255)
                                                                                                                                  gaagcgTCATGAGTTCATGAC‑CGTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTg              >  2:2374343/6‑150 (MQ=255)
                                                                                                                                  gaagcgTCATGAGTTCATGAC‑CGTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTg              <  1:2374343/145‑1 (MQ=255)
                                                                                                                                  gaagcgTCATGAGTTCATGAC‑CGTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTg              <  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  <  2:1531483/150‑1 (MQ=255)
                                                                                                                                               cTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGc  <  1:531532/150‑1 (MQ=255)
                                                                                                                                               cTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGc  <  1:281525/150‑1 (MQ=255)
                                                                                                                                               cTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGc  <  1:1161164/150‑1 (MQ=255)
                                                                                                                                               cTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACAGCCGGGCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGc  <  1:1001806/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
ACAGGGATCAACGGCGTGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGC  >  NC_002947/4519352‑4519644

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: ‑