Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,926,513 T→G 5.3% D118A (GAC→GCC)  PP_1726 ← ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,926,5130TG 5.3% 240.9 / 2.6 75D118A (GAC→GCC) PP_1726ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (40/31);  new base G (2/2);  total (42/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 = 9.87e-01

GACCTTGTACTTGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCGGTGTC  >  NC_002947/1926365‑1926661
                                                                                                                                                    |                                                                                                                                                      
gACCTTGTACTTGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGt                                                                                                                                                        <  6:321168/149‑1 (MQ=255)
 aCCTTGTACTTGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTc                                                                                                                                                       <  3:88630/149‑1 (MQ=255)
    ttGTACTTGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTcc                                                                                                                                                      >  2:337765/1‑147 (MQ=255)
     tGTACTTGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAg                                                                                                                                                  >  4:80795/1‑149 (MQ=255)
            gCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgcttg                                                                                                                                           >  6:124646/1‑149 (MQ=255)
            gCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATACGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGGACCAGGTgccg                                                                                                                                            <  6:299338/148‑3 (MQ=255)
              cTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgcttggc                                                                                                                                         <  5:133961/149‑1 (MQ=255)
               tttCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTCTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCACAGTGGCCGTCCTTGTCCTTGGCCCATTCCGTCACCTGGTCCCAGGTgcttggct                                                                                                                                        >  8:8967/1‑149 (MQ=255)
                             gCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctggct                                                                                                                           >  8:27974/1‑148 (MQ=255)
                              ccAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGa                                                                                                                         >  7:47748/1‑149 (MQ=255)
                                      tCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCtt                                                                                                                  <  5:124646/148‑1 (MQ=255)
                                         gCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGt                                                                                                               >  3:299207/1‑148 (MQ=255)
                                            gTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCAcc                                                                                                           >  6:327245/1‑149 (MQ=255)
                                             ttCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTg                                                                                                                 >  2:241283/1‑142 (MQ=255)
                                             ttCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTg                                                                                                                 <  1:241283/142‑1 (MQ=255)
                                                 tctTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGAt                                                                                                      >  2:168720/1‑149 (MQ=255)
                                                   ttCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCTCGAt                                                                                                      <  6:24455/147‑1 (MQ=255)
                                                   ttCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCg                                                                                                    >  5:61521/1‑149 (MQ=255)
                                                       ttGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgctt                                                                                                                                            <  1:330224/105‑1 (MQ=255)
                                                       ttGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgctt                                                                                                                                            >  2:330224/1‑105 (MQ=255)
                                                            cAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAgg                                                                                           <  6:61521/149‑1 (MQ=255)
                                                            cAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAgg                                                                                           <  4:299207/149‑1 (MQ=255)
                                                                  cTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCa                                                                                      <  3:80795/148‑1 (MQ=255)
                                                                        cgccatgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCAc                                                                                 <  6:214763/143‑1 (MQ=255)
                                                                         atgatgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCGTTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCAc                                                                                 <  6:113094/146‑1 (MQ=255)
                                                                         atgatgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCGTTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCAc                                                                                 >  5:113094/1‑146 (MQ=255)
                                                                         atgatgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGt                                                                               >  1:303126/1‑148 (MQ=255)
                                                                         atgatgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGt                                                                               >  2:133663/1‑148 (MQ=255)
                                                                            atgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTccc                                                                                                                                                     >  8:327935/1‑76 (MQ=255)
                                                                            atgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTccc                                                                                                                                                     <  7:327935/76‑1 (MQ=255)
                                                                            atgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCAcc                                                                           >  6:141767/1‑149 (MQ=255)
                                                                             tgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctgg                                                                                                                             >  8:147867/1‑98 (MQ=255)
                                                                             tgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctgg                                                                                                                             <  7:147867/98‑1 (MQ=255)
                                                                             tgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCAcc                                                                           <  8:47748/148‑1 (MQ=255)
                                                                               aaaGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGAt                                                                        <  1:337765/149‑1 (MQ=255)
                                                                                     aTGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGgc                                                                  <  1:168720/149‑1 (MQ=255)
                                                                                            cGGTATAGGCCAGGGGCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCg                                                           <  5:59693/149‑1 (MQ=255)
                                                                                               tataGGCCAGCGCCCAGTGTCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGTTCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTg                                                       <  8:291794/149‑1 (MQ=255)
                                                                                                   ggCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGCTGTCGGCGCTGGCGTGGTCCt                                                   >  8:273842/1‑149 (MQ=255)
                                                                                                   ggCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTcc                                                    <  1:133663/148‑1 (MQ=255)
                                                                                                    gCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCt                                                               <  2:58810/137‑1 (MQ=255)
                                                                                                    gCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCt                                                               >  1:58810/1‑137 (MQ=255)
                                                                                                          gcCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGc                                            >  2:178144/1‑149 (MQ=255)
                                                                                                             cAGTGGCCGTCCTTGTCCTAGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTCAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTc                                         >  6:260071/1‑149 (MQ=255)
                                                                                                                 ggCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAAc                                     >  8:358/1‑149 (MQ=255)
                                                                                                                  gCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACt                                    >  3:331214/1‑149 (MQ=255)
                                                                                                                   ccctccttgtccatgGCCCAATCCCGCACCTGGTCCCAGCTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGCTCACCCCGATCGGGACGAAGCCGGCACCCACGTCACCGATGTAGGCGCTGGACTTGTCCTTATCGGCATCGAACtt                                   >  4:82165/4‑149 (MQ=255)
                                                                                                                   ccGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAAGGTGCTTGGCTTGTAAGGATGGCTGACGCCCTTGGACACCGCGATCGGGCCGAAGCCGGCACCCACGAAACCGATGTCGGAGCTGGCGTTGTCCTTATCGGCTTCGACCt                                    >  1:245720/1‑149 (MQ=255)
                                                                                                                    cGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGcc                                                                                                >  1:204944/1‑88 (MQ=255)
                                                                                                                    cGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGcc                                                                                                <  2:204944/88‑1 (MQ=255)
                                                                                                                    cGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTc                                         <  4:154181/142‑1 (MQ=255)
                                                                                                                    cGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTc                                         >  3:154181/1‑142 (MQ=255)
                                                                                                                    cGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTg                                  <  2:303126/149‑1 (MQ=255)
                                                                                                                        cttgtccttgGCCCATTCCGGCAGCTGGTCCGAGGTGCGTGTGTTGTATGGCTGGCTGACGCGCTTGGTCAGGGCGATGGGGCCGAAGGCGGCACCCACGTGACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTgggc                               <  3:82165/148‑3 (MQ=255)
                                                                                                                           gtccttgGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGtt                           >  5:329596/1‑149 (MQ=255)
                                                                                                                               ttgGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTg                       >  7:300879/1‑149 (MQ=255)
                                                                                                                                 ggCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTggg                     >  4:56900/1‑149 (MQ=255)
                                                                                                                                 ggCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTgg                      <  7:27974/148‑1 (MQ=255)
                                                                                                                                   cccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTtct                                                >  6:215232/1‑120 (MQ=255)
                                                                                                                                   cccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTtct                                                <  5:215232/120‑1 (MQ=255)
                                                                                                                                   cccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGcc                   >  6:120369/1‑149 (MQ=255)
                                                                                                                                    ccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTg                       <  5:67877/144‑1 (MQ=255)
                                                                                                                                    ccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTg                       >  6:67877/1‑144 (MQ=255)
                                                                                                                                    ccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGcc                   >  6:262574/1‑148 (MQ=255)
                                                                                                                                    ccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCATCGAACTTGGCCAGTTCCTGGGcc                   >  4:68195/1‑148 (MQ=255)
                                                                                                                                      aTTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAg                >  2:307054/1‑149 (MQ=255)
                                                                                                                                         caagCACCTGGGACCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACtt                                   <  2:164403/124‑1 (MQ=255)
                                                                                                                                         caagCACCTGGGACCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCt              >  8:59945/4‑148 (MQ=255)
                                                                                                                                          cGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCa            >  6:184344/1‑149 (MQ=255)
                                                                                                                                          aagCACCTGGGACCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGATGGCGTACAGAATACACGGCTACGAACTAGGCAAGTTCCTGGGCCgccaac             >  1:302722/3‑144 (MQ=255)
                                                                                                                                           ggCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAg                <  2:190548/144‑1 (MQ=255)
                                                                                                                                           ggCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAg                >  1:190548/1‑144 (MQ=255)
                                                                                                                                           ggCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCAt           >  5:56743/1‑149 (MQ=255)
                                                                                                                                               ccTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCg       <  6:56743/149‑1 (MQ=255)
                                                                                                                                               ccTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCg       <  8:300879/149‑1 (MQ=255)
                                                                                                                                                   gTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCGgtgt   <  5:327245/149‑1 (MQ=255)
                                                                                                                                                    tCCCAGGGGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGTCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTTGCCAGTTCCTGGGCCGAGCTCATGTAGGTGTc  <  7:8967/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GACCTTGTACTTGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCGGTGTC  >  NC_002947/1926365‑1926661

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: