Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,568,696 T→G 8.6% E284D (GAA→GAC)  PP_2252 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,568,6960TG 8.6% 208.0 / 7.5 69E284D (GAA→GAC) PP_2252hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/23);  new base G (3/3);  total (43/26)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CTCCACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTC  >  NC_002947/2568549‑2568842
                                                                                                                                                    |                                                                                                                                                  
ctCCACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑Agttgtt                                                                                                                                                   <  2:223492/149‑1 (MQ=255)
   cACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGa                                                                                                                                                >  7:2433/1‑149 (MQ=255)
   cACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGa                                                                                                                                                >  2:190468/1‑149 (MQ=255)
             ggcgCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGATGGTCGAACAACTGCgg                                                                                                                                      <  2:338305/149‑1 (MQ=255)
              gcgcACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCggg                                                                                                                                     <  8:270670/149‑1 (MQ=255)
                 cACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGggtgg                                                                                                                                  <  6:286324/149‑1 (MQ=255)
                  aCTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGggtggt                                                                                                                                 >  7:45571/1‑149 (MQ=255)
                     ggCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCAc                                                                                                                              >  8:253762/1‑149 (MQ=255)
                       ccAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACgg                                                                                                                            <  8:2433/149‑1 (MQ=255)
                        cAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGt                                                                                                                           >  2:263250/1‑149 (MQ=255)
                                  tCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtca                                                                                                                  >  3:156598/1‑148 (MQ=255)
                                  tCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtc                                                                                                                   <  1:190468/147‑1 (MQ=255)
                                    cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtcatc                                                                                                                <  8:45571/148‑1 (MQ=255)
                                     tCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAc                                                                                                              >  7:105816/1‑149 (MQ=255)
                                     tCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAc                                                                                                              >  3:112050/1‑149 (MQ=255)
                                       aaCCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAcc                                                                                                             <  1:263250/148‑1 (MQ=255)
                                         ccTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCg                                                                                                          <  3:150365/149‑1 (MQ=255)
                                          cTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCg                                                                                                          >  5:304940/1‑148 (MQ=255)
                                              cccGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAAcc                                                                                                     >  6:291128/1‑149 (MQ=255)
                                                      gTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCgg                                                                                             <  6:304940/149‑1 (MQ=255)
                                                        aaTAGACCCCCGCCATGGCCCCGCTCTGGCGGGTCTCCATCAACTCCTGCCATTCCTGGTTCACCGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGGCATCACCCCGTAACCGGCTGCGGc                                                                                            >  4:31673/1‑148 (MQ=255)
                                                         aTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGc                                                                                            <  5:149481/147‑1 (MQ=255)
                                                         aTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGc                                                                                            >  6:149481/1‑147 (MQ=255)
                                                           agccccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCCCCCCGTAACCGGCTGCGGCAATc                                                                                        >  5:301961/4‑149 (MQ=255)
                                                           aGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATc                                                                                        >  2:166556/1‑149 (MQ=255)
                                                             aCCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCg                                                                                       >  3:34306/1‑148 (MQ=255)
                                                              cccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGa                                                                                      >  7:76536/1‑148 (MQ=255)
                                                               ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACg                                                                                                                             <  8:199928/108‑1 (MQ=255)
                                                               ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACg                                                                                                                             <  8:194697/108‑1 (MQ=255)
                                                               ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACg                                                                                                                             >  7:199928/1‑108 (MQ=255)
                                                               ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACg                                                                                                                             >  7:194697/1‑108 (MQ=255)
                                                               ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGaaa                                                                                    <  7:253762/149‑1 (MQ=255)
                                                               ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGaaa                                                                                    >  8:214345/1‑149 (MQ=255)
                                                                      aTGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGc                                                                             >  2:245184/1‑149 (MQ=255)
                                                                           cccGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTg                                                                        >  4:136358/1‑149 (MQ=255)
                                                                                    gCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCtt                                                                 <  7:214345/147‑1 (MQ=255)
                                                                                     cGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGcc                                                              <  8:76536/149‑1 (MQ=255)
                                                                                          ctcCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTg                                                         >  2:187864/1‑149 (MQ=255)
                                                                                             cATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGcc                                                       >  1:266078/1‑149 (MQ=255)
                                                                                             cATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGcc                                                       >  5:301811/1‑148 (MQ=255)
                                                                                             cATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGcc                                                       >  7:178100/1‑148 (MQ=255)
                                                                                                    tCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCAc                                               >  3:149152/1‑149 (MQ=255)
                                                                                                      cTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCAc                                               <  4:156598/147‑1 (MQ=255)
                                                                                                      cTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACgg                                             >  2:82704/1‑149 (MQ=255)
                                                                                                       tGCCATTCCTGGTTCAACGACACGTTCAGGATAACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACTGTCATGGGTCATCACCCCGTAAACGGCAGCGGCAATCGAACTCATCGCCCCTGCCACCTTGCACACTGCCGATGCCAAGGc                                            >  4:104497/1‑149 (MQ=255)
                                                                                                         ccATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCAc                                          >  3:275051/1‑149 (MQ=255)
                                                                                                               cTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATg                                                   <  6:122196/134‑1 (MQ=255)
                                                                                                               cTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATg                                                   >  5:122196/1‑134 (MQ=255)
                                                                                                               cTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGa                                    >  7:274057/1‑149 (MQ=255)
                                                                                                                 ggTTCAACGCCACGTTCAGGATCACCCGC‑AGTAGTTCGACCATCTGCTGGTGGTCACGGTCATGGGTCATAACCCCGTAACCGCCTACGGCAATCGAAGTCATAGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTTGCCATc                                  >  7:190341/1‑149 (MQ=255)
                                                                                                                     cAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCcgc                              >  4:77637/1‑149 (MQ=255)
                                                                                                                      aaCGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCcgc                              <  4:149152/148‑1 (MQ=255)
                                                                                                                       aCGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGt                            <  2:26787/149‑1 (MQ=255)
                                                                                                                        cGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTg                           <  6:332830/149‑1 (MQ=255)
                                                                                                                              gTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGc                       >  1:184170/1‑147 (MQ=255)
                                                                                                                                tCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCTTGGCGCCacc                   >  5:111533/1‑147 (MQ=255)
                                                                                                                                   tgaccACCCGC‑AGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTg                <  8:249516/145‑1 (MQ=255)
                                                                                                                                    gATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCAcc                                                                   >  4:247523/1‑97 (MQ=255)
                                                                                                                                    gATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCAcc                                                                   <  3:247523/97‑1 (MQ=255)
                                                                                                                                    gATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTgg               >  8:186319/1‑149 (MQ=255)
                                                                                                                                      tCACGCTC‑AGTTGTTCTACCATCTGCGTTTGGTCACGGTCATGGGTCATCAGCCCTTATCCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACATTGCCCACTGCAGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgc              <  8:190341/148‑1 (MQ=255)
                                                                                                                                      tCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgcg             <  4:112050/149‑1 (MQ=255)
                                                                                                                                       cACCCGC‑AGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgc              >  2:90715/1‑147 (MQ=255)
                                                                                                                                       cACCCGC‑AGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgc              <  1:90715/147‑1 (MQ=255)
                                                                                                                                        aCCCGC‑AGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgcgc            >  1:41304/1‑148 (MQ=255)
                                                                                                                                         cccGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGTCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGcc           >  6:147197/1‑148 (MQ=255)
                                                                                                                                         cccGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGcc           >  2:149325/1‑148 (MQ=255)
                                                                                                                                          ccGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGcc           >  1:106070/1‑147 (MQ=255)
                                                                                                                                          ccGC‑AGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCCCTGGCGATCCCGCGGGGGGCCAGGGTGGCGcc           >  4:255232/1‑147 (MQ=255)
                                                                                                                                            gC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGc       <  8:105816/149‑1 (MQ=255)
                                                                                                                                                   gttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTc  <  4:275051/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTCCACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGC‑AGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTC  >  NC_002947/2568549‑2568842

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑