Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,387,896 T→C 5.7% *198W (TGA→TGG)  PP_1209 ← cold‑shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,387,8960TC 5.7% 209.2 / 2.2 70*198W (TGA→TGG) PP_1209cold‑shock protein
Reads supporting (aligned to +/- strand):  ref base T (38/28);  new base C (2/2);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

GAGGCCCCGCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCG  >  NC_002947/1387751‑1388037
                                                                                                                                                  |                                                                                                                                                
gagGCCCCGCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCaa                                                                                                                                               <  5:182915/149‑1 (MQ=255)
    gccatccgcggtgagggccaccgcattctggtcgaaggccagcgcgtggagttttcggtgatgcaccgcgacaagggcctgcaggccgaagacgtggttgcagtaaaccGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAAcggc                                                                                                                                           >  6:55084‑M1/110‑149 (MQ=255)
     cccGCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCCCAGGTATGATGCTGCGTGCAGGTGCAGCCCGGCCGTTGTGGCCGGGGGTTTGCCGTTGATAGCTAG‑GGCCGGAACGCCGTTGGTAAAcggcg                                                                                                                                          >  8:308043/1‑149 (MQ=255)
      ccGCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAAcggcgg                                                                                                                                         <  6:77063/149‑1 (MQ=255)
        gCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGACCAAcggctgac                                                                                                                                       <  7:174496/149‑5 (MQ=255)
        gCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGACCAAcggctgac                                                                                                                                       <  2:96596/149‑5 (MQ=255)
               cATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGc                                                                                                                                 >  6:77001/1‑148 (MQ=255)
                   attGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCa                                                                                                                            >  5:143380/1‑149 (MQ=255)
                    ttGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCAc                                                                                                                           >  1:306720/1‑149 (MQ=255)
                    ttGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCAc                                                                                                                           >  2:250768/1‑149 (MQ=255)
                    ttGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCAc                                                                                                                           <  6:143380/149‑1 (MQ=255)
                       aGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACt                                                                                                                                   >  6:281090/1‑138 (MQ=255)
                       aGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACt                                                                                                                                   <  5:281090/138‑1 (MQ=255)
                       aGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTc                                                                                                                        >  3:264264/1‑149 (MQ=255)
                           gttCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTc                                                                                                                     <  5:77001/148‑1 (MQ=255)
                            ttCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCgg                                                                                                                   >  4:151998/1‑149 (MQ=255)
                               cAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCt                                                                                                                >  7:42261/1‑149 (MQ=255)
                                aGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTg                                                                                                               <  2:232082/149‑1 (MQ=255)
                                  gctgctATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACCGCAACCACGTCTTCGGCCTGCa                                                                                                             <  7:50729/149‑1 (MQ=255)
                                   ctgctATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAg                                                                                                            >  2:303615/1‑149 (MQ=255)
                                      ctATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAgg                                                                                                           <  3:151998/147‑1 (MQ=255)
                                       tATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGccc                                                                                                        >  5:146063/1‑149 (MQ=255)
                                       tATCGCCGGCAACCCATGCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGccc                                                                                                        <  7:308043/149‑1 (MQ=255)
                                              ggCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTcgc                                                                                                 <  8:26266/149‑1 (MQ=255)
                                              ggCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTcgc                                                                                                 >  8:269213/1‑149 (MQ=255)
                                                cAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCgg                                                                                               <  4:60968/149‑1 (MQ=255)
                                                 aaGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGt                                                                                              <  6:146063/149‑1 (MQ=255)
                                                   gccagccCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCACCCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑g                                                                                           >  5:99363/1‑148 (MQ=255)
                                                   gccagccCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑g                                                                                           >  1:155870/1‑148 (MQ=255)
                                                     cagccCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGGCGGGACAGCCTTGGGTCAACGGCGGTTTCATGAACCCCCGTGTGCGGCGTGCGGGCCATTGACGCGGG‑‑G‑c                                                                                         >  2:229279/1‑148 (MQ=255)
                                                        cccccACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCAc                                                                                    <  3:274087/149‑1 (MQ=255)
                                                               aGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAAc                                                                             >  2:44966/1‑149 (MQ=255)
                                                                  tATGATGCTGCCTTCAGGTGTAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTc                                                                           >  7:79517/1‑148 (MQ=255)
                                                                    tgatgCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTcgc                                                                                                 >  6:280857/1‑127 (MQ=255)
                                                                    tgatgCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTcgc                                                                                                 <  5:280857/127‑1 (MQ=255)
                                                                         gacaagggcctgcaggccgaagacgtggttgcagtaaaccGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAAcggcgg                                                                                                                                         <  2:227578‑M1/42‑1 (MQ=255)
                                                                         gacaagggcctgcaggccgaagacgtggttgcagtaaaccGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAAcggcgg                                                                                                                                         >  1:227578‑M1/41‑82 (MQ=255)
                                                                         ctgcCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCt                                                                   >  1:203157/1‑149 (MQ=255)
                                                                          tgcCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTg                                                                                                     >  3:76696/1‑117 (MQ=255)
                                                                          tgcCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTg                                                                                                     <  4:76696/117‑1 (MQ=255)
                                                                           gcCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTgg                                                                 >  4:286756/1‑149 (MQ=255)
                                                                             cTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGcc                                                               <  7:269213/149‑1 (MQ=255)
                                                                               tCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCtt                                                             >  1:153397/1‑149 (MQ=255)
                                                                                    tGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGAc                                                         >  8:306697/1‑148 (MQ=255)
                                                                                      cagcaCGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGAc                                                         >  8:83627/1‑146 (MQ=255)
                                                                                       agcaCGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGc                                                                                                                  <  4:174137/91‑1 (MQ=255)
                                                                                       agcaCGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGc                                                                                                                  <  4:169922/91‑1 (MQ=255)
                                                                                       agcaCGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGc                                                                                                                  >  3:174137/1‑91 (MQ=255)
                                                                                       agcaCGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGc                                                                                                                  >  3:169922/1‑91 (MQ=255)
                                                                                          aCGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATg                                                  >  8:253105/1‑149 (MQ=255)
                                                                                           cGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCAcc                                                                                   <  7:225824/115‑1 (MQ=255)
                                                                                           cGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCAcc                                                                                   >  8:225824/1‑115 (MQ=255)
                                                                                              cccTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑GCCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCgg                                               <  1:229279/149‑1 (MQ=255)
                                                                                               ccTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTg                                             <  8:79517/149‑1 (MQ=255)
                                                                                               ccTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTg                                             >  7:2096/1‑149 (MQ=255)
                                                                                                 ttGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAAc                                                                             >  1:10824/1‑115 (MQ=255)
                                                                                                 ttGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAAc                                                                             <  2:10824/115‑1 (MQ=255)
                                                                                                   gtgGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGcc                                          <  4:264264/148‑1 (MQ=255)
                                                                                                      ggCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGGGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCGCTc                                       <  2:153397/148‑1 (MQ=255)
                                                                                                             gTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCgg                                <  7:253105/148‑1 (MQ=255)
                                                                                                                  CCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGcc                          >  8:42261‑M1/1‑149 (MQ=255)
                                                                                                                   CGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGccc                         <  2:112928‑M1/149‑1 (MQ=255)
                                                                                                                            cGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTAGAGCAACCACGTCTTCGGCCTGCAGGCCCTTGTAGCGGT‑‑G‑CCTCACCGAAAACTCCCCGCCCTGGCCTTCGACCAGAATGCGGTTGCCCTCCCCGCGGATTGCCCGAAATTGCa                >  1:215784/1‑149 (MQ=255)
                                                                                                                            cGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCa                <  3:286756/149‑1 (MQ=255)
                                                                                                                                 ggCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTAGG‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCAc               >  1:266988/1‑148 (MQ=255)
                                                                                                                                  gCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACg              >  4:43904/1‑146 (MQ=255)
                                                                                                                                      gTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGACTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATg        >  2:141572/1‑148 (MQ=255)
                                                                                                                                        ccgCCGTTGACCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCa     >  2:240577/3‑149 (MQ=255)
                                                                                                                                        ccgCCGTTGACCAAAGGAGGTTTACTGCAACCACGTCTTCGGACTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCGCGCAAGTGCACGATGCTGTct     >  4:51246/3‑148 (MQ=255)
                                                                                                                                           ccGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCg  >  3:242871/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GAGGCCCCGCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTC‑GGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGT‑‑G‑CATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCG  >  NC_002947/1387751‑1388037

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: ‑