Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,029,111 T→C 7.2% L284P (CTG→CCG)  wbpM → polysaccharide biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,029,1110TC 7.2% 202.2 / 4.7 69L284P (CTG→CCG) wbpMpolysaccharide biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base C (2/3);  total (33/36)
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.98e-01

CAGCAGATGATCAAGGAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCATCCACAAGCAAACGGTACTGGTCACCGGGGCCGGTGGCTCC  >  NC_002947/2028972‑2029256
                                                                                                                                            |                                                                                                                                                 
cagcagATGATCAAGGAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTg                                                                                                                                           <  7:204285/148‑1 (MQ=255)
cagcagATGATCAAGGAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTg                                                                                                                                           >  8:204285/1‑148 (MQ=255)
   cagATGATCAAGGAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCg                                                                                                                                       >  2:260422/1‑149 (MQ=255)
    agATGATCAAGGAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCgg                                                                                                                                      <  1:172309/149‑1 (MQ=255)
        gatCAAGGAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTc                                                                                                                                   >  7:255610/1‑148 (MQ=255)
         atCAAGGAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCCGGCCACTTGCCAGGTcgg                                                                                                                                 <  4:337461/149‑3 (MQ=255)
          tCAAGGAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCa                                                                                                                                  >  1:17701/1‑147 (MQ=255)
            aaGGAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGt                                                                                                                              >  1:279305/1‑149 (MQ=255)
            aaGGAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGt                                                                                                                              <  2:220844/149‑1 (MQ=255)
             aGGAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCaa                                                                                                                                              >  8:74697/1‑132 (MQ=255)
             aGGAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCaa                                                                                                                                              <  7:74697/132‑1 (MQ=255)
               gAAACCGGTGCATAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGa                                                                                                                           >  6:225852/1‑149 (MQ=255)
               gAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTTCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGa                                                                                                                           <  4:10708/149‑1 (MQ=255)
               gAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGa                                                                                                                           <  3:168242/149‑1 (MQ=255)
               gAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTCGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGa                                                                                                                           >  6:293961/1‑149 (MQ=255)
                aaaCCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGAt                                                                                                                          <  2:103382/149‑1 (MQ=255)
                      gTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATc                                                                                                                    <  4:51839/149‑1 (MQ=255)
                       tGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATcc                                                                                                                   <  2:17701/149‑1 (MQ=255)
                         cacaGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATcc                                                                                                                   >  7:84614/1‑147 (MQ=255)
                         cacaGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCaa                                                                                                                 <  2:259925/149‑1 (MQ=255)
                         cacaGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCCGGCCACTTGCCAGGTcg                                                                                                                                  <  8:160732/132‑2 (MQ=255)
                         cacaGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCCGGCCACTTGCCAGGTcg                                                                                                                                  >  7:160732/1‑131 (MQ=255)
                          acaGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATcc                                                                                                                   >  5:257041/1‑146 (MQ=255)
                          acaGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAg                                                                                                                <  1:260422/149‑1 (MQ=255)
                           caGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGa                                                                                                               <  1:219242/149‑1 (MQ=255)
                            aGGAAATCCTCCTCGCCCTGCCTTCGGTATGGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGAc                                                                                                              >  1:34748/1‑149 (MQ=255)
                                 aTCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTCGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGa                                                                                                         <  5:293961/149‑1 (MQ=255)
                                   cctcctCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCaaa                                                                                                                                >  6:104078/1‑124 (MQ=255)
                                   cctcctCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCaaa                                                                                                                                <  5:104078/124‑1 (MQ=255)
                                     tcctcGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACAtt                                                                                                     >  5:335969/1‑149 (MQ=255)
                                       ctcGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTg                                                                                                    <  6:324755/148‑1 (MQ=255)
                                          gCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGcc                                                                                                  <  8:255610/147‑1 (MQ=255)
                                                 cTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCtt                                                                                         <  6:335969/149‑1 (MQ=255)
                                                  ttCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTg                                                                                        >  6:118287/1‑149 (MQ=255)
                                                     ggTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGc                                                                                     <  4:176932/149‑1 (MQ=255)
                                                       tATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGt                                                                                   >  7:348142/1‑149 (MQ=255)
                                                       tATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGt                                                                                   <  8:84614/149‑1 (MQ=255)
                                                           gcgcgcCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGAc                                                                                <  5:225852/148‑1 (MQ=255)
                                                                  cGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTAcc                                                                        <  6:257041/149‑1 (MQ=255)
                                                                        aaGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCaa                                                                  <  2:279305/149‑1 (MQ=255)
                                                                           gagaTTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGAc                                                                                >  4:252834/1‑132 (MQ=255)
                                                                           gagaTTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGAc                                                                                <  3:252834/132‑1 (MQ=255)
                                                                           gagaTTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCaagg                                                                >  5:8759/1‑148 (MQ=255)
                                                                            agaTTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCaaggaa                                                              <  8:348142/149‑1 (MQ=255)
                                                                               ttCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAggg                                                           <  5:118287/149‑1 (MQ=255)
                                                                                tctcGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGgc                                                          >  2:58791/1‑149 (MQ=255)
                                                                                    gCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCAGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGgctgc                                                       >  1:252731/1‑149 (MQ=255)
                                                                                    gCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGgctgct                                                      >  4:202535/1‑149 (MQ=255)
                                                                                               aCCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCa                                           >  8:300305/1‑149 (MQ=255)
                                                                                                  ttACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCATc                                         >  7:56318/1‑148 (MQ=255)
                                                                                                   tACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACAt                                                                                                                     >  8:255714/1‑71 (MQ=255)
                                                                                                   tACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACAt                                                                                                                     <  7:255714/71‑1 (MQ=255)
                                                                                                    aCCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCATCca                                       >  2:44651/1‑148 (MQ=255)
                                                                                                        gTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCATCCACAAGc                                  <  2:34748/149‑1 (MQ=255)
                                                                                                             aCGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCaagg                                                                <  5:264183/114‑1 (MQ=255)
                                                                                                             aCGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCaagg                                                                >  6:264183/1‑114 (MQ=255)
                                                                                                             aCGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCATCCACAAGCAAACg                             <  5:78392/149‑1 (MQ=255)
                                                                                                                 gcgcAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCaaggaag                                                             <  2:300575/113‑1 (MQ=255)
                                                                                                                             cGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCATCCACAAGCAAACGGTACTGGTCACCggg              >  1:24211/1‑149 (MQ=255)
                                                                                                                                 cTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTTCCCGCCAAGGAAGGGCTGCTCGACCGCTGCGTCCACAAGCAAACGGTACTGGTCACCGGGGGCg          >  4:245891/1‑149 (MQ=255)
                                                                                                                                  tttGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCATCCACAAGCAAACGGTACTGGTCACCGGGGCCgg         >  2:351738/1‑149 (MQ=255)
                                                                                                                                  tttGCCGACCTGGAAAGTGGCCGCGTCAAAGTGTATGAAATCCAATACGAGGACATTGCCCATATCCTTGTGCGTGACACTGTAAAAGACAAGTAAGGGCTGCTCGAACGCTGCATCCACAAGCAAACGGTACTAGTCACCTGTGccta         >  4:50123/1‑147 (MQ=255)
                                                                                                                                   actttgACCCGGCCACTTGCCAGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTc                                                     >  3:131760/6‑104 (MQ=38)
                                                                                                                                   actttgACCCGGCCACTTGCCAGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTc                                                     <  4:131760/99‑1 (MQ=38)
                                                                                                                                     gCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCATCCACAAGCAAACGGTACTGGTCACCGGGGCCGGTgg      >  1:251890/1‑149 (MQ=255)
                                                                                                                                       cGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCATCCACAAGCAAACGGTACTGGTCa                   <  7:257917/134‑1 (MQ=255)
                                                                                                                                       cGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCATCCACAAGCAAACGGTACTGGTCa                   >  8:257917/1‑134 (MQ=255)
                                                                                                                                       cGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCATCCACAAGCAAACGGTACTGGTCACCGGGGCCGGTGGCt    <  6:8759/149‑1 (MQ=255)
                                                                                                                                         aCCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCATCCACAAGCAAACGGTACTGGTCACCGGGGCCGGTGGCTcc  <  1:44651/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
CAGCAGATGATCAAGGAAACCGGTGCACAGGAAATCCTCCTCGCCCTGCCTTCGGTATCGCGCGCCCGGCGCAAGGAGATTCTCGCCTTGCTCGAACCTTACCCGTTGCACGTGCGCAGCGTGCCGGCCTTTGCCGACCTGGCAAGTGGCCGGGTCAAAGTGGATGACATCCAAGACGTGGACATTGCCGATTTGCTTGGGCGTGACCCGGTACCCGCCAAGGAAGGGCTGCTCGAACGCTGCATCCACAAGCAAACGGTACTGGTCACCGGGGCCGGTGGCTCC  >  NC_002947/2028972‑2029256

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: