Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,406,053 C→G 5.3% A185A (GCG→GCC PP_4750 ← amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,406,0530CG 5.3% 245.0 / 3.7 75A185A (GCG→GCCPP_4750amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base G (2/2);  total (39/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.76e-01

ACAGCACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGG  >  NC_002947/5405906‑5406195
                                                                                                                                                    |                                                                                                                                              
aCAGCACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCAccgccg                                                                                                                                               >  1:42140/1‑149 (MQ=255)
       gaaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                        <  4:153953/148‑1 (MQ=255)
       aaaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGTCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                        >  3:160534/1‑149 (MQ=255)
       aaaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                        >  8:352607/1‑149 (MQ=255)
        aaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGTCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgtg                                                                                                                                       <  4:160534/149‑1 (MQ=255)
        aaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgtg                                                                                                                                       <  8:5155/149‑1 (MQ=255)
             cTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtg                                                                                                                                       <  8:272346/144‑1 (MQ=255)
             cTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtg                                                                                                                                       >  7:222719/1‑144 (MQ=255)
             cTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtg                                                                                                                                       <  8:222719/144‑1 (MQ=255)
             cTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtg                                                                                                                                       >  7:272346/1‑144 (MQ=255)
                 cacaACTCCGGCAACCATTGCATCCACATCGTTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCt                                                                                                                               >  1:231247/1‑149 (MQ=255)
                 cacaACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTc                                                                                                                              <  1:358781/149‑1 (MQ=255)
                       tCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGccc                                                                                                                        <  7:352607/149‑1 (MQ=255)
                         cGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCt                                                                                                                       >  5:179995/1‑148 (MQ=255)
                            cAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGaaa                                                                                                                   <  3:135004/149‑1 (MQ=255)
                               ccATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGt                                                                                                                <  2:42140/149‑1 (MQ=255)
                                  ttGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCGGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTg                                                                                                             <  3:181117/149‑1 (MQ=255)
                                     cATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGtttca                                                                                                            >  7:327909/1‑145 (MQ=255)
                                      aTCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCtgt                                                                                                         <  3:111353/149‑1 (MQ=255)
                                         cacaTCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGa                                                                                                       <  7:124326/148‑1 (MQ=255)
                                         cacaTCGTCAGCAGTTCGCTCATAGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACAGGCACGCTGAGCACCAGGTACACCACCGACGCCGCGGGGTAGACGATCAGGCCCTGAAAGGTTTGCTgaga                                                                                                       >  8:67153/1‑148 (MQ=255)
                                          acaTCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAAAGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGaaa                                                                                                     >  6:128549/1‑149 (MQ=255)
                                            aTCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAAcc                                                                                                   <  2:231247/149‑1 (MQ=255)
                                             tCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCt                                                                                                  >  3:255927/1‑149 (MQ=255)
                                             tCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCt                                                                                                  >  8:93467/1‑149 (MQ=255)
                                               gTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTgg                                                                                                <  6:179995/149‑1 (MQ=255)
                                                 tCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTgg                                                                                                >  1:241322/1‑148 (MQ=255)
                                                   agcagTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGt                                                                                               >  5:22267/1‑147 (MQ=255)
                                                    gcagTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAg                                                                                            >  1:197477/1‑149 (MQ=255)
                                                     cagTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGc                                                                                           >  7:190030/1‑149 (MQ=255)
                                                       gTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGACGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCTGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCt                                                                                         <  7:67153/149‑1 (MQ=255)
                                                        ttCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCt                                                                                         >  3:38366/1‑148 (MQ=255)
                                                          cGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGa                                                                                       <  4:255927/148‑1 (MQ=255)
                                                           gCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCt                                                                                                           >  1:154623/1‑127 (MQ=255)
                                                           gCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCt                                                                                                           <  2:154623/127‑1 (MQ=255)
                                                           gCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGaaa                                                                                     <  2:197477/149‑1 (MQ=255)
                                                             tCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGc                                                                                   >  3:41374/1‑149 (MQ=255)
                                                                tCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGa                                                                                <  4:41374/149‑1 (MQ=255)
                                                                  gggCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGAt                                                                               <  5:128549/148‑1 (MQ=255)
                                                                    gCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTg                                                                                                 >  2:94204/1‑128 (MQ=255)
                                                                        tCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCAcc                                                                         >  4:358012/1‑148 (MQ=255)
                                                                         ccTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCaa                                                                       >  7:64333/1‑149 (MQ=255)
                                                                           tctGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACg                                                                     >  4:227678/1‑149 (MQ=255)
                                                                               tttCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCaa                                                                 >  8:259622/1‑149 (MQ=255)
                                                                                ttCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCaa                                                                 <  3:227678/148‑1 (MQ=255)
                                                                                ttCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAAt                                                                >  3:48697/1‑149 (MQ=255)
                                                                                    aCGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAAtcgtc                                                            <  6:22267/149‑1 (MQ=255)
                                                                                       tGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAg                                                         <  7:93467/149‑1 (MQ=255)
                                                                                         gCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGc                                                       <  3:358012/149‑1 (MQ=255)
                                                                                               aGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGc                                                 >  7:63015/1‑149 (MQ=255)
                                                                                                  cAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATTGTCGAGGCCACGGcc                                                >  6:83510/1‑147 (MQ=255)
                                                                                                  cAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGcc                                                >  3:38545/1‑147 (MQ=255)
                                                                                                   aaCGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAg                                             <  2:241322/149‑1 (MQ=255)
                                                                                                      gagagaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAgg                                                        <  7:70954/135‑1 (MQ=255)
                                                                                                      gagagaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAgg                                                        >  8:70954/1‑135 (MQ=255)
                                                                                                        gagaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCtt                                        <  8:63015/149‑1 (MQ=255)
                                                                                                          gaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGa                                      <  7:259622/149‑1 (MQ=255)
                                                                                                                ggCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAg                                  <  4:94474/147‑1 (MQ=255)
                                                                                                                ggCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAg                                  >  3:94474/1‑147 (MQ=255)
                                                                                                                      tGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAg                                  >  2:230494/1‑141 (MQ=255)
                                                                                                                      tGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAg                                  <  1:230494/141‑1 (MQ=255)
                                                                                                                          aCGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAAGCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCg                        >  5:231153/1‑147 (MQ=255)
                                                                                                                          aCGCAGAGCACCGCGTACACCACCGCCGCCCCCCTGTAGACGCTCAGGCCCTGAAAGGTATGCACTCAAACCTGGTAACCCTGAAAGCAGATTTCACCAACGCCAATCCTCGAGGCCACCGCCGAGACCAAGACCAGGCCacagccg                        >  2:278292/1‑141 (MQ=255)
                                                                                                                              tGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTg                                                                                                 <  4:56579/70‑1 (MQ=255)
                                                                                                                              tGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTg                                                                                                 >  3:56579/1‑70 (MQ=255)
                                                                                                                                  cACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAgg                    <  5:69445/143‑1 (MQ=255)
                                                                                                                                  cACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAgg                    >  6:69445/1‑143 (MQ=255)
                                                                                                                                  cACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCt               >  4:204727/1‑148 (MQ=255)
                                                                                                                                     cAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAgg           >  8:81043/1‑149 (MQ=255)
                                                                                                                                     cAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAgg           >  8:86172/1‑149 (MQ=255)
                                                                                                                                           caccacCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGaa     <  8:64333/149‑1 (MQ=255)
                                                                                                                                            accaccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAAt    <  7:86172/149‑1 (MQ=255)
                                                                                                                                            accaccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAAt    <  7:81043/149‑1 (MQ=255)
                                                                                                                                              caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCACCGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg  >  7:355867/1‑149 (MQ=255)
                                                                                                                                              caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg  >  4:300616/1‑149 (MQ=255)
                                                                                                                                              caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg  <  8:190030/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ACAGCACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGG  >  NC_002947/5405906‑5406195

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: