Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 307,865:1 +C 91.3% intergenic (+880/‑989) hslO → / → PP_0254 molecular chaperone HslO/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947307,8651.C91.3% 218.8 / 14.7 69intergenic (+880/‑989)hslO/PP_0254molecular chaperone HslO/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base C (35/28);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

AAAGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTA  >  NC_002947/307717‑308006
                                                                                                                                                        |                                                                                                                                              
aaaGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAg                                                                                                                                                 <  6:233118/149‑1 (MQ=255)
  aGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGcc                                                                                                                                              <  8:8272/149‑2 (MQ=255)
    ccATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTTTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGccgt                                                                                                                                            <  7:3710/149‑4 (MQ=255)
    ccATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGTCAGCCGAAAGGGta                                                                                                                                           <  2:300111/149‑2 (MQ=255)
            tcGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGAt                                                                                                                                    >  5:205512/1‑149 (MQ=255)
                  cAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTg                                                                                                                               <  7:196605/148‑1 (MQ=255)
                          ccTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGAcgacga                                                                                                                        <  2:14440/147‑1 (MQ=255)
                          ccTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAg                                                                                                                       <  6:147738/148‑1 (MQ=255)
                            tgctgcCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGCAAGCCGTTACCTGATCGCTGACGACGACCAc                                                                                                                    >  6:246605/1‑149 (MQ=255)
                               tgccggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGt                                                                                                                 >  6:196969/1‑149 (MQ=255)
                                 ccggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGtt                                                                                                                <  5:196969/148‑1 (MQ=255)
                                      tgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGc                                                                                                           >  4:26943/1‑148 (MQ=255)
                                      cgccgatgttggccgcgcagtgcataggcagtacatcggcGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTga                                                                                                                              <  3:53694‑M2/89‑1 (MQ=255)
                                      cgccgatgttggccgcgcagtgcataggcagtacatcggcGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTga                                                                                                                              >  4:53694‑M2/41‑129 (MQ=255)
                                       gccgatgttggccgcgcagtgcataggcagtacatcggcGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCg                                                                                                          >  1:300662‑M2/40‑148 (MQ=255)
                                         CGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAg                                                                                                       >  3:176410‑M1/1‑149 (MQ=255)
                                              TACTGCCTATGGACTGCGCGGCCAACATCGGCTAAGAATGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                  >  8:99723‑M1/1‑149 (MQ=255)
                                                cTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTg                                                                                                 >  3:71553/1‑148 (MQ=255)
                                                   ccTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCAGTTACATGATCGGTGAAGACGAGCACGGTTGGGGCGAAGGCGTTGTCtt                                                                                             >  8:124748/1‑149 (MQ=255)
                                                   ccTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACAATATGGGGTGCTGCAGTGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTc                                                                                               >  1:111468/1‑149 (MQ=255)
                                                    cTATGCACTGCGCGGCCAACATCGGCGAAGTAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTTACGTAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTc                                                                                            <  5:246605/149‑1 (MQ=255)
                                                    cTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                  >  7:7241/1‑143 (MQ=255)
                                                    cTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                  <  8:7241/143‑1 (MQ=255)
                                                     tATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTc                                                                                            >  5:194002/1‑148 (MQ=255)
                                                        gCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACa                                                                                        >  1:6010/1‑149 (MQ=255)
                                                         cACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACAt                                                                                       <  4:71553/149‑1 (MQ=255)
                                                          aCTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATg                                                                                      <  6:205512/149‑1 (MQ=255)
                                                             ataggcagtacatcggcGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGaa                                                                                   <  2:300662‑M2/132‑1 (MQ=255)
                                                               gcgGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAgg                                                                                 <  4:312737/149‑1 (MQ=255)
                                                                  gCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAggcgg                                                                              >  5:97981/1‑149 (MQ=255)
                                                                     aaCATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTc                                                                                               <  5:122137/129‑1 (MQ=255)
                                                                     aaCATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTc                                                                                               >  6:122137/1‑129 (MQ=255)
                                                                     aaCATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGgctg                                                                           >  3:322025/1‑149 (MQ=255)
                                                                        atcggcGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTa                                                                        >  4:148440‑M2/7‑149 (MQ=255)
                                                                           ggcGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGc                                                                     >  4:263643‑M2/4‑149 (MQ=255)
                                                                                    ggCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCAt                                                            <  3:263643/149‑1 (MQ=255)
                                                                                       gATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGa                                                         >  5:265122/1‑149 (MQ=255)
                                                                                        aTGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGAc                                                        <  6:265122/149‑1 (MQ=255)
                                                                                        aTGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGAc                                                        <  6:194002/149‑1 (MQ=255)
                                                                                        aTGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGAc                                                        <  6:97981/149‑1 (MQ=255)
                                                                                                  gTTCTTCGGCCTGTCCTGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCt                                                      <  8:27447/141‑1 (MQ=255)
                                                                                                  gTTCTTCGGCCTGTCCTGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCt                                                      >  7:27447/1‑141 (MQ=255)
                                                                                                   ttcttcGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAggc                                                                                >  8:213142/1‑114 (MQ=255)
                                                                                                   ttcttcGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAggc                                                                                <  7:213142/114‑1 (MQ=255)
                                                                                                   ttcttcGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCgag                                               >  3:93717/1‑147 (MQ=255)
                                                                                                     cttcGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGagaag                                            >  5:258989/1‑148 (MQ=255)
                                                                                                        cGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAAc                                         <  6:167482/148‑1 (MQ=255)
                                                                                                          gCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGt                                                                >  8:139435/1‑123 (MQ=255)
                                                                                                          gCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGt                                                                <  7:139435/123‑1 (MQ=255)
                                                                                                             tGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCg                                   >  7:162782/1‑149 (MQ=255)
                                                                                                               tCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCgg                                  <  4:322025/148‑1 (MQ=255)
                                                                                                               tCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGt                                 >  7:301370/1‑149 (MQ=255)
                                                                                                                ccGGCACTGGCCACACAACCCTTTCGGCTGACGAAACCCGTTACCTCAGCGGTGACGGCGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCgggc                                <  4:255448/149‑3 (MQ=255)
                                                                                                                 cGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGGTACGCCAAGTGCATCGACCTGTGCGAGAAGAACGAGCCGGTcg                               >  8:189122/1‑148 (MQ=255)
                                                                                                                    cACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCt                            <  6:258989/149‑1 (MQ=255)
                                                                                                                         gCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAg                                                                                                       >  3:209419/1‑69 (MQ=255)
                                                                                                                         gCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAg                                                                                                       <  4:209419/69‑1 (MQ=255)
                                                                                                                         gCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGaaa                       <  8:162782/149‑1 (MQ=255)
                                                                                                                             gACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAgg                                                                                 >  4:197746/1‑87 (MQ=255)
                                                                                                                               cTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTTACGACGAGCACGGTTGGGGCTAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATc                 <  7:189122/149‑1 (MQ=255)
                                                                                                                               cTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCAt                  >  3:227371/1‑148 (MQ=255)
                                                                                                                               cTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCAt                  >  3:69162/1‑148 (MQ=255)
                                                                                                                                tACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATc                 >  7:282849/1‑148 (MQ=255)
                                                                                                                                  cccTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCaa                                                                  >  6:81848/1‑97 (MQ=255)
                                                                                                                                  cccTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAg              >  7:280748/1‑149 (MQ=255)
                                                                                                                                   ccTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGc             <  8:301370/149‑1 (MQ=255)
                                                                                                                                    cTTTCGGCTGACGAAATCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCa            <  7:124748/149‑1 (MQ=255)
                                                                                                                                       tCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCAc           <  8:280748/147‑1 (MQ=255)
                                                                                                                                            tGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAg    >  5:60562/1‑149 (MQ=255)
                                                                                                                                               aCGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAg              <  8:233546/137‑1 (MQ=255)
                                                                                                                                               aCGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAg              >  7:233546/1‑137 (MQ=255)
                                                                                                                                               aCGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTa  >  6:310532/1‑149 (MQ=255)
                                                                                                                                                        |                                                                                                                                              
AAAGCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTA  >  NC_002947/307717‑308006

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: