Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,990,454 C→G 17.1% P71P (CCG→CCC ycaC‑III ← putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,990,4540CG17.1% 172.4 / 23.5 69P71P (CCG→CCCycaC‑IIIputative hydrolase
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base G (8/4);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 5.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGACAGCAGGC  >  NC_002947/5990307‑5990596
                                                                                                                                                   |                                                                                                                                              
cACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCg                                                                                                                                               <  1:72648/149‑1 (MQ=255)
cACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCg                                                                                                                                               >  3:245317/1‑149 (MQ=255)
          tCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCgg                                                                                                                                     >  5:77607/1‑149 (MQ=255)
           cGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCggg                                                                                                                                    >  3:386495/1‑149 (MQ=255)
                acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCAtt                                                                                                                               >  6:324430/1‑149 (MQ=17)
                 caacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTg                                                                                                                              >  1:25163/1‑149 (MQ=255)
                  aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCAttcg                                                                                                                             >  3:152974/1‑147 (MQ=17)
                  aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTgg                                                                                                                             >  4:383305/1‑149 (MQ=255)
                   acacCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg                                                                                                                            <  6:220060/149‑1 (MQ=255)
                     gccggCAATGATCCGCGGCTTCTTGCCGGTCGCCTTCTCCGCCTTGACGATATCCTCGTTGTCCGAGGCGGTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTGCAGCTCGGGCACCAGCGGGCCTTTcgg                                                                                                                            <  4:295587/143‑1 (MQ=14)
                       cTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCt                                                                                                                        <  7:76632/149‑1 (MQ=17)
                       cTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCAGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCt                                                                                                                        >  8:359423/1‑149 (MQ=255)
                       cTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGAGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCt                                                                                                                        >  4:325053/1‑149 (MQ=255)
                           aaTGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCt                                                                                                                    >  8:320704/1‑149 (MQ=14)
                               aTCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGaa                                                                                                                <  4:257665/149‑1 (MQ=255)
                                  aGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCt                                                                                                             <  3:171109/149‑1 (MQ=14)
                                  aGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCt                                                                                                             >  7:134456/1‑149 (MQ=14)
                                     tgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCAGCCGGGAACATATCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggt                                                                                                          >  7:77892/1‑149 (MQ=14)
                                        ttcttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTggtggt                                                                                                       >  2:240079/1‑149 (MQ=11)
                                        atcttGCCGGTCGCCTTCACCCCCTTCACGAATTCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGCCAATGTACGGCGCATCCGGGAACAGGCCTTTCAGTTCCGGCACCAGCCGGCCGTTGGTGCCTAGCTCGAAGCTggtggt                                                                                                       >  4:260478/2‑149 (MQ=17)
                                         tcttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTggaggt                                                                                                       >  5:30779/1‑148 (MQ=255)
                                             gCCGGTCGCCTTCCCCGCCTTGACGAAATCCTCGTTGGCCCAGGCGTTGATCTGGCCAGGGCTGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGa                                                                                                     <  2:153216/146‑1 (MQ=11)
                                             gCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGa                                                                                                     <  2:334273/146‑1 (MQ=255)
                                             gCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGGACGGCGCCTCCGGGAACATCTCTTTCAGCTCGGGCACCCGCGGGCCATTCTGGCCTTGCTCGAAGCTGGGGGTGAGGa                                                                                                  >  8:368906/1‑149 (MQ=11)
                                              ccGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGAt                                                                                                 <  4:386495/149‑1 (MQ=255)
                                                ggTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATgg                                                                                               <  1:377355/149‑1 (MQ=255)
                                                 gTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                                                                                              >  1:316964/1‑149 (MQ=11)
                                                 gTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGGTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                                                                                              <  6:30779/149‑1 (MQ=255)
                                                    gCCTTCACCGCCTTGACGAAATCCTCGTGGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCgg                                                                                           <  3:259737/149‑1 (MQ=255)
                                                      cTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACACGCCTAACAGTTCAGGCACCAGCGGGCCGATGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCgg                                                                                           >  1:322985/1‑147 (MQ=17)
                                                         cACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGt                                                                                      >  6:335820/1‑149 (MQ=255)
                                                             gCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTgaa                                                                                  <  5:335820/149‑1 (MQ=255)
                                                              ccTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTg                                                                                    >  7:364214/1‑146 (MQ=255)
                                                                  gACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAAc                                                                              <  2:25163/148‑1 (MQ=255)
                                                                    cGATATCCTCGTTGGCGCTGCCGTTCATCTGGCCAGGGCTGGCACTGTACGGCGCATCCGTGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGGCTTGCTGGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGATGAACTTg                                                                           <  4:249704/149‑1 (MQ=18)
                                                                     gAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTgg                                                                          >  1:333474/1‑149 (MQ=255)
                                                                     gAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTgg                                                                          >  1:381087/1‑149 (MQ=255)
                                                                     gAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCACGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTgg                                                                          >  3:175994/1‑149 (MQ=255)
                                                                           cTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTTGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAgg                                                                     <  4:245317/148‑1 (MQ=255)
                                                                             cGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGt                                                                    >  5:313100/1‑147 (MQ=255)
                                                                              gTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTca                                                                  <  6:313100/148‑1 (MQ=255)
                                                                                   cccAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGCCCAGGTcagccagc                                                            >  5:31719/1‑149 (MQ=255)
                                                                                    ccAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACtt                                                                            >  5:111913/1‑132 (MQ=255)
                                                                                    ccAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACtt                                                                            <  6:111913/132‑1 (MQ=255)
                                                                                    ccAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAgccagcg                                                           >  7:88732/1‑149 (MQ=255)
                                                                                          gTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCa                                                     >  8:182845/1‑149 (MQ=255)
                                                                                             gATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACgt                                                  >  4:141164/1‑149 (MQ=255)
                                                                                              aTCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCt                                                                                                                    >  5:349300/1‑82 (MQ=255)
                                                                                              aTCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCt                                                                                                                    <  6:349300/82‑1 (MQ=255)
                                                                                                cTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACgttg                                                <  7:359423/148‑1 (MQ=255)
                                                                                                   gCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACgttgtt                                              <  3:383305/147‑1 (MQ=255)
                                                                                                   gCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCt                                            >  5:342720/1‑149 (MQ=255)
                                                                                                     cAGGGCGGGCAATGTACGGCGCATCCGGGAAGCGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTGGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAAGTTGGCGAGGTCAGCCAGCGCCAGCACGTTGTTCtt                                           >  4:122731/1‑148 (MQ=17)
                                                                                                      aGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGa                                         <  7:335374/149‑1 (MQ=255)
                                                                                                          cgggcAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACtt                                     <  6:342720/149‑1 (MQ=255)
                                                                                                            ggcAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCAAGAAcaacg                                   <  4:175994/149‑5 (MQ=255)
                                                                                                             gcAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTc                                  >  1:204342/1‑149 (MQ=255)
                                                                                                              cAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAg                                                                                                    >  1:275495/1‑82 (MQ=255)
                                                                                                              cAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAg                                                                                                    <  2:275495/82‑1 (MQ=255)
                                                                                                              cAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCt                                 <  2:381087/149‑1 (MQ=255)
                                                                                                               aaTGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTg                                <  3:260478/149‑1 (MQ=255)
                                                                                                                aTGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTgg                               >  3:117240/1‑149 (MQ=255)
                                                                                                                         gcATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATc                      <  2:333474/149‑1 (MQ=255)
                                                                                                                         gcATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATc                      >  1:111085/1‑149 (MQ=255)
                                                                                                                               gggAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCAc                >  1:197760/1‑149 (MQ=255)
                                                                                                                                      gggCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGAc         <  3:325053/149‑1 (MQ=255)
                                                                                                                                            tCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACgttg                                                <  2:248571/104‑1 (MQ=255)
                                                                                                                                               gTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAgccagcgcca                                                        <  6:331208/93‑1 (MQ=255)
                                                                                                                                               gTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGACAGCAGGc  >  2:91802/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGACAGCAGGC  >  NC_002947/5990307‑5990596

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