Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 966,599 A→T 6.9% L41M (TTG→ATG) ‡ PP_0828 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947966,5990AT 6.9% 510.8 / 21.2 159L41M (TTG→ATG) ‡PP_0828membrane protein
Reads supporting (aligned to +/- strand):  ref base A (78/70);  new base T (7/4);  total (85/74)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 165 total aligned reads displayed.

GCTGCGCGAACTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATG  >  NC_002947/966450‑966747
                                                                                                                                                     |                                                                                                                                                    
         aCTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCTTGAACAGCaa                                                                                                                                             <  1:4214043/150‑2 (MQ=255)
         aCTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCTTGAACAGCaa                                                                                                                                             <  2:3804458/150‑2 (MQ=255)
             ttCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGaaa                                                                                                                                         <  2:3953444/150‑1 (MQ=255)
               gaggctcgcgcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATc                                                                                                                                       >  2:268444/7‑150 (MQ=255)
               cACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATc                                                                                                                                       >  1:1480829/1‑150 (MQ=255)
                aCCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCg                                                                                                                                      >  1:4112168/1‑150 (MQ=255)
                aCCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCg                                                                                                                                      <  1:4587375/150‑1 (MQ=255)
                aCCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCg                                                                                                                                      >  1:2219997/1‑150 (MQ=255)
                  cAGCGCGCCCCACCGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGAt                                                                                                                                    <  2:1529821/150‑1 (MQ=255)
                  cAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGAt                                                                                                                                    <  2:4599606/150‑1 (MQ=255)
                  cAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGAt                                                                                                                                    >  2:3707561/1‑150 (MQ=255)
                  cAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGAt                                                                                                                                    >  2:3707639/1‑150 (MQ=255)
                   ctcgcgcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATg                                                                                                                                   <  1:268444/148‑1 (MQ=255)
                     cgcgcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGt                                                                                                                                 <  2:4895583/150‑1 (MQ=255)
                        gcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGAc                                                                                                                              >  1:527141/1‑150 (MQ=255)
                        gcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGAc                                                                                                                              >  1:607581/1‑150 (MQ=255)
                        gcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGAc                                                                                                                              >  1:1227780/1‑150 (MQ=255)
                        gcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGAc                                                                                                                              <  1:4294817/150‑1 (MQ=255)
                          cccACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCa                                                                                                                            >  1:3473499/1‑150 (MQ=255)
                          cccACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCa                                                                                                                            >  1:3635703/1‑150 (MQ=255)
                           ccACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAg                                                                                                                           >  2:255996/1‑150 (MQ=255)
                            cacaGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGc                                                                                                                          >  1:2832338/1‑150 (MQ=255)
                            cacaGGGTGCGAGCCTCCCACCAGCGGCGGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGc                                                                                                                          >  2:1896521/1‑150 (MQ=255)
                             acaGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGca                                                                                                                         >  2:2441120/1‑150 (MQ=255)
                                    gCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCa                                                                                                                  <  1:548757/150‑1 (MQ=255)
                                      gAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGa                                                                                                                <  1:4539158/150‑1 (MQ=255)
                                         ccTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATca                                                                                                             <  1:878718/150‑1 (MQ=255)
                                         ccTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATca                                                                                                             <  2:3137299/150‑1 (MQ=255)
                                          cTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATcag                                                                                                            >  2:2350625/1‑150 (MQ=255)
                                              caccaGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCagcagc                                                                                                        >  1:1533711/1‑150 (MQ=255)
                                              caccaGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCagcagc                                                                                                        >  1:1533513/1‑150 (MQ=255)
                                              caccaGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCagcagc                                                                                                        >  2:2960507/1‑150 (MQ=255)
                                                ggtatagCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCa                                                                                                      <  1:4024424/144‑1 (MQ=255)
                                                 caGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAg                                                                                                     >  2:599484/1‑150 (MQ=255)
                                                     ggCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACa                                                                                                 <  1:2313129/150‑1 (MQ=255)
                                                       cTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGc                                                                                               <  1:3049841/150‑1 (MQ=255)
                                                       cTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGc                                                                                               <  1:1338062/150‑1 (MQ=255)
                                                          tAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGAGGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGAc                                                                                            <  1:1550290/150‑1 (MQ=255)
                                                           aGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGAcc                                                                                           <  1:4427478/150‑1 (MQ=255)
                                                              cgcTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGc                                                                                        <  2:515731/150‑1 (MQ=255)
                                                                 tgctgTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGa                                                                                     <  2:2102896/150‑1 (MQ=255)
                                                                 cgctgTAGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGGGACCAGCACCTCCCAGCTCAGCAGCCAGAACAGCGACCAGCCGa                                                                                     <  2:1086893/149‑1 (MQ=255)
                                                                   ctgTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGata                                                                                   >  2:3279965/1‑150 (MQ=255)
                                                                    tgTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGatat                                                                                  >  2:4144843/1‑150 (MQ=255)
                                                                          gAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGa                                                                            <  2:396093/150‑1 (MQ=255)
                                                                          gAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGa                                                                            >  1:3585988/1‑150 (MQ=255)
                                                                            aaaCTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCGGATATAGCCGATg                                                                          <  2:1946094/150‑1 (MQ=255)
                                                                             aaCTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACa                                                                                                 <  2:4233621/126‑1 (MQ=255)
                                                                             aaCTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATgg                                                                         <  1:2477050/150‑1 (MQ=255)
                                                                              aCTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGt                                                                        >  1:1771031/1‑150 (MQ=255)
                                                                                tCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCt                                                                      >  1:3247647/1‑150 (MQ=255)
                                                                                  aCCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGCTGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGTGACCAGCCGATATAGCCGATGGTCTTc                                                                    <  1:1734941/150‑1 (MQ=255)
                                                                                  aCCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTc                                                                    >  1:2319312/1‑150 (MQ=255)
                                                                                   ccAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTca                                                                   <  1:2634259/150‑1 (MQ=255)
                                                                                    cAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTcac                                                                  <  2:3999595/150‑1 (MQ=255)
                                                                                      gCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTcacca                                                                >  1:1847907/1‑150 (MQ=255)
                                                                                          caccaGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAg                                                            >  2:559797/1‑150 (MQ=255)
                                                                                          caccaGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAg                                                            >  2:3895488/1‑150 (MQ=255)
                                                                                          caccaGCGCTGAGCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTTAACTTGCTGTTCAAGAACAGCATGAAATAGATGGTGACCAGCACATACCAGAACAGCAGCCAGAACAGCGACCAGCCGAAATAGCCGATGGTCATCACCAACAgg                                                           >  1:3754825/1‑150 (MQ=255)
                                                                                           accaGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAgg                                                           >  2:2585218/1‑150 (MQ=255)
                                                                                           accaGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCACGGTGCACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACCGCGACCAGCCGAGATAGCCGATGGTCTTCACCACCggg                                                           >  1:4335463/1‑150 (MQ=255)
                                                                                            ccaGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGCTGGTCTTCACCACCAGGc                                                          <  1:3279965/150‑1 (MQ=255)
                                                                                             caGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCg                                                         >  2:4528888/1‑150 (MQ=255)
                                                                                              aGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCgg                                                        >  2:233934/1‑150 (MQ=255)
                                                                                                 gcTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTAc                                                     >  1:4539691/1‑150 (MQ=255)
                                                                                                   tgtgCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACtt                                                   >  1:1932836/4‑150 (MQ=255)
                                                                                                   tGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATcagca                                                                                                          >  2:425465/1‑95 (MQ=255)
                                                                                                   tGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATcagca                                                                                                          <  1:425465/95‑1 (MQ=255)
                                                                                                     aGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGc                                                 <  2:3977020/150‑1 (MQ=255)
                                                                                                       cccAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCgg                                               <  1:1710323/150‑1 (MQ=255)
                                                                                                       cccAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCgg                                               >  2:3836882/1‑150 (MQ=255)
                                                                                                        ccAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCCGCAGCCAGAACAGCGACCCGCCGATATAGCCGATGGTCGTCACCACCAGGCGGGAcgggcggg                                              >  1:4711727/1‑150 (MQ=255)
                                                                                                        ccAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCggg                                              >  2:4064617/1‑150 (MQ=255)
                                                                                                          agcagGGTCAGCGGGATCAGTGGCAGGGTGAACTTGATGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCg                                            <  2:955699/150‑1 (MQ=255)
                                                                                                            cagGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGACATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGAt                                          <  2:3928153/150‑1 (MQ=255)
                                                                                                            cagGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGAt                                          >  1:876906/1‑150 (MQ=255)
                                                                                                            cagGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGAt                                          >  2:258074/1‑150 (MQ=255)
                                                                                                                   aGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCt                                   >  2:4474652/1‑150 (MQ=255)
                                                                                                                   aGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTAGCGGGCGATGACGTCt                                   >  2:331780/1‑150 (MQ=255)
                                                                                                                        gATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAg                                                                                         <  2:713803/91‑1 (MQ=255)
                                                                                                                        gATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAg                                                                                         >  1:713803/1‑91 (MQ=255)
                                                                                                                           cAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCAcgcg                           >  1:1560544/1‑150 (MQ=255)
                                                                                                                              tGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAgg                        <  1:3833055/150‑1 (MQ=255)
                                                                                                                                gCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGt                      >  2:84790/1‑150 (MQ=255)
                                                                                                                                  aGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGACCATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTAc                    >  2:4351436/1‑150 (MQ=255)
                                                                                                                                    ggTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCt                  >  2:1200727/1‑150 (MQ=255)
                                                                                                                                    ggTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCt                  >  2:1200725/1‑150 (MQ=255)
                                                                                                                                       tttcatGCTGTTCTTGAACAGCAAGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGa               >  1:3499884/6‑150 (MQ=255)
                                                                                                                                           ttGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAAcg           >  1:266161/1‑150 (MQ=255)
                                                                                                                                           atGCTGTTCTTGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAAcg           >  1:5008176/2‑150 (MQ=255)
                                                                                                                                            tGCTGTTCTTGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAAcgc          >  1:5024344/1‑150 (MQ=255)
                                                                                                                                               tGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCggg       <  2:4042084/150‑1 (MQ=255)
                                                                                                                                               tGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCggg       <  2:2832338/150‑1 (MQ=255)
                                                                                                                                               tGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCggg       <  2:3941353/150‑1 (MQ=255)
                                                                                                                                                 ttCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCggggg     >  1:2101790/1‑150 (MQ=255)
                                                                                                                                                   cAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCgg        >  1:4966377/1‑145 (MQ=255)
                                                                                                                                                   cAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGAt   <  2:266161/150‑1 (MQ=255)
                                                                                                                                                    aaGAACATCAAGAAATCGATGGTGACCAGAACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGTGGGATg  <  2:3635703/150‑1 (MQ=255)
                                                                                                                                                    aaGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATg  >  2:84951/1‑150 (MQ=255)
                                                                                                                                                    aaGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGACGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATg  >  2:4445925/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTGCGCGAACTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATG  >  NC_002947/966450‑966747

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