Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,519,505 T→G 5.8% K20T (AAG→ACG)  clpA ← ATP‑dependent serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,519,5050TG 5.8% 219.6 / 4.0 78K20T (AAG→ACG) clpAATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base T (33/40);  new base G (3/2);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

AACGGCGTGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGCGTATTGGTT  >  NC_002947/4519361‑4519653
                                                                                                                                                |                                                                                                                                                    
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)
                                                                                                                         gaagcgTCATGAGTTCATGACCGTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTg                       <  1:2374186/145‑1 (MQ=255)
                                                                                                                         gaagcgTCATGAGTTCATGACCGTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTg                       <  1:2374343/145‑1 (MQ=255)
                                                                                                                         gaagcgTCATGAGTTCATGACCGTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTg                       >  2:2374343/6‑150 (MQ=255)
                                                                                                                         gaagcgTCATGAGTTCATGACCGTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTg                       >  2:2374186/6‑150 (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:531532/150‑1 (MQ=255)
                                                                                                                                      cTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGc           <  1:281525/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                                    <  2:2393295/117‑1 (MQ=255)
                                                                                                                                              gCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACaga                                    >  1:2393295/1‑117 (MQ=255)
                                                                                                                                              gCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGCGTATTGGt   <  2:1021833/150‑1 (MQ=255)
                                                                                                                                               cTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTCACGCACCATCTCAATTTCACATATTACCTGATGCTGGCTTTCCCTGGCGTATTGGtt  <  1:299421/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AACGGCGTGGTGGAGTCGATGAACTCTTGCAGGTCATGCTTGAGCTTGTCGAGATTGGCGCCACAGGCGCGCAGAACGGTCGCGGCAGCCTCATTGTCAAGGAGTGCCAGCAGCAGATGCTCGACGGTCATGAACTCATGACGCTTCGAACGGGCCTCCTTGAAGGCAAGATTGAGGGTGACTTCGAGCTCGCGGTTTAACATAGCTTCACCTCATACCCAAGTGGTCGGCGATTAACCGTCCTTCTCGATTTCACAGAGTAGCGGATGCTGGCTTTCCCTGGCGTATTGGTT  >  NC_002947/4519361‑4519653

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: