Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,039,566 T→C 6.1% D141D (GAT→GAC PP_3560 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,039,5660TC 6.1% 191.0 / 2.2 65D141D (GAT→GACPP_3560LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base C (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CCTGGGGGCGATCCACACGGTGCTGGCGCGGCGTTTGCCCAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGC  >  NC_002947/4039422‑4039713
                                                                                                                                                |                                                                                                                                                   
ccTGGGGGCGATCCACACGGTGCTGGCGCGGCGTTTGCCCAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAg                                                                                                                                                 <  6:100267/149‑1 (MQ=255)
 cTGGGGGCGATCCACACGGTGCTGGCGCGGCGTTTGCCCAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGACGCAGc                                                                                                                                                >  7:265161/1‑149 (MQ=255)
     gggCGATCCACACGGTGCTGGCGCGGCGTTTGCCCAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGAt                                                                                                                                             >  1:220193/1‑148 (MQ=255)
     gggCGATCCACACGGTGCTGGCGCGGCGTTTGCCCAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATc                                                                                                                                            <  6:217639/149‑1 (MQ=255)
              acacGGTGCTGGCGCGGCGTTTGCCCAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGaa                                                                                                                                   <  2:105124/149‑1 (MQ=255)
                    tGCTGGCGCGGCGTTTGCCCAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTc                                                                                                                             <  2:220193/149‑1 (MQ=255)
                      cTGGCGCGGCGTTTGCCCAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTATACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCa                                                                                                                            <  3:257544/148‑1 (MQ=255)
                         gcgcGGCGGTTGCCCAAGGCGCTGGAGTCGAGCAAGGCCCAGCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGccc                                                                                                                        <  8:212129/149‑1 (MQ=255)
                          cgcgGCGTTTGCCCAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAAc                                                                                                                                  >  8:61158/1‑138 (MQ=255)
                          cgcgGCGTTTGCCCAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAAc                                                                                                                                  <  7:61158/138‑1 (MQ=255)
                              gcgTTTGCCCAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTAcc                                                                                                                   >  8:127384/1‑149 (MQ=255)
                                gTTTGCCCAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACGACCGAACCGGTCCGCCCTTACCCg                                                                                                                 >  3:95060/1‑149 (MQ=255)
                                     cccAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGAAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGaa                                                                                                            >  1:6150/1‑149 (MQ=255)
                                        aaGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCt                                                                                                         >  3:193157/1‑149 (MQ=255)
                                         aGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTg                                                                                                        <  4:193157/149‑1 (MQ=255)
                                             gcTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACt                                                                                                    >  7:232960/1‑149 (MQ=255)
                                               tGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTc                                                                                                  <  1:62739/149‑1 (MQ=255)
                                                  tgtgGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTcaca                                                                                               >  3:201322/1‑149 (MQ=255)
                                                    tgGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGAt                                                                                                                                                     >  3:230367/1‑93 (MQ=255)
                                                    tgGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGAt                                                                                                                                                     <  4:230367/93‑1 (MQ=255)
                                                     ggATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCg                                                                                            >  3:119086/1‑149 (MQ=255)
                                                         cAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCt                                                                                                         <  5:263731/132‑1 (MQ=255)
                                                         cAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCt                                                                                                         >  6:263731/1‑132 (MQ=255)
                                                           aGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTc                                                                                      >  2:82441/1‑149 (MQ=255)
                                                           aGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGACGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGAc                                                                                                     <  7:12093/134‑1 (MQ=255)
                                                           aGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGACGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGAc                                                                                                     >  8:12093/1‑134 (MQ=255)
                                                             gCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGa                                                                                    <  1:39498/149‑1 (MQ=255)
                                                                  tCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCGCCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGAcc                                                                               <  5:142991/149‑1 (MQ=255)
                                                                   cACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCg                                                                              >  2:308623/1‑149 (MQ=255)
                                                                      ccGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGttt                                                                           >  8:260100/1‑149 (MQ=255)
                                                                          aGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTgg                                                                       <  1:174363/149‑1 (MQ=255)
                                                                              acacATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCAGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTc                                                                                                  >  4:117436/1‑118 (MQ=255)
                                                                              acacATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCAGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTc                                                                                                  <  3:117436/118‑1 (MQ=255)
                                                                                ttctCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGGTAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCAtt                                                                 <  5:180601/146‑1 (MQ=255)
                                                                                acaTCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGa                                                                                 >  4:249155/1‑133 (MQ=255)
                                                                                acaTCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGa                                                                                 <  3:249155/133‑1 (MQ=255)
                                                                                acaTCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCAtt                                                                 >  4:183612/1‑149 (MQ=255)
                                                                                  aTCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGc                                                               >  6:58607/1‑149 (MQ=255)
                                                                                  aTCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGc                                                               >  8:180355/1‑149 (MQ=255)
                                                                                   tCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGcc                                                              <  3:183612/149‑1 (MQ=255)
                                                                                    cAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGcc                                                              <  8:232960/148‑1 (MQ=255)
                                                                                     aGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAg                                                            >  5:37068/1‑149 (MQ=255)
                                                                                     aGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTAGAGGACCGTTTCTGGGCCAtt                                                                 >  6:180601/1‑144 (MQ=255)
                                                                                            ccTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTTAGCTCGATGCAGCGATCACCACCGGACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAg                                                     >  1:98877/1‑149 (MQ=255)
                                                                                              tCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGc                                                               <  7:89726/137‑1 (MQ=255)
                                                                                              tCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGc                                                               >  8:89726/1‑137 (MQ=255)
                                                                                                 ggCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCCGAGCTGGc                                                >  2:87900/1‑149 (MQ=255)
                                                                                                  gCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGcc                                               >  5:133811/1‑149 (MQ=255)
                                                                                                    aTGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCgg                                             <  2:6150/149‑1 (MQ=255)
                                                                                                      gTCGGCGCAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCACTTCGAGGACCTTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCggg                                            >  4:53172/1‑148 (MQ=255)
                                                                                                       tCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCTCAGGGCTGGcaatggc                                          <  8:20749/149‑7 (MQ=255)
                                                                                                       tCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGACGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCTCAGGGCTGGcaatggc                                          <  8:265161/149‑7 (MQ=255)
                                                                                                         ggcggAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGtata                                        <  7:127384/149‑1 (MQ=255)
                                                                                                                cTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGcc                                 <  7:260100/149‑1 (MQ=255)
                                                                                                                  ggcgcggcgcGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGa                                                                                    <  5:98244/96‑1 (MQ=255)
                                                                                                                  ggcgcggcgcGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGa                                                                                    >  6:98244/1‑96 (MQ=255)
                                                                                                                    cgcggcgcgTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAg                             >  4:213441/1‑149 (MQ=255)
                                                                                                                     gcggcgcgTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGt                            >  4:39662/1‑149 (MQ=255)
                                                                                                                        gcgcgTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCCGAGCTGGCCGGGTATAGCGTGCCGCAGTTGc                         <  1:87900/149‑1 (MQ=255)
                                                                                                                          gcgTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTg                       >  1:293847/1‑149 (MQ=255)
                                                                                                                               gAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGgccag                  >  6:317272/1‑149 (MQ=255)
                                                                                                                                aaGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTGCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGgccagc                 <  5:85397/149‑1 (MQ=255)
                                                                                                                                 aGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGgccagcc                <  6:134457/149‑1 (MQ=255)
                                                                                                                                         aGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTc        <  1:308623/149‑1 (MQ=255)
                                                                                                                                              gATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTg     <  2:293847/147‑1 (MQ=255)
                                                                                                                                               aTGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTgcgc  <  4:201322/149‑1 (MQ=255)
                                                                                                                                               aTGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGCACCGTTTCTGGGCCATTGCCAGCCCCGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTgcgc  <  1:82441/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CCTGGGGGCGATCCACACGGTGCTGGCGCGGCGTTTGCCCAAGGCGCTGGTGTGGATCAAGGCCCATCACCCGCAGCTACACATCAGCGTGGCCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGC  >  NC_002947/4039422‑4039713

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: