Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,406,056 C→G 8.1% A184A (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,0560CG 8.1% 189.5 / 6.8 62A184A (GCG→GCCPP_4750amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base G (2/3);  total (31/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 = 3.37e-01

AGCACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCA  >  NC_002947/5405908‑5406204
                                                                                                                                                    |                                                                                                                                                    
aGCACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCAccgccgcc                                                                                                                                                      >  5:51372/1‑149 (MQ=255)
  cacaAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACcgccgccgc                                                                                                                                                    >  2:129896/1‑149 (MQ=255)
     aaaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                                 >  5:133088/1‑149 (MQ=255)
      aaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCTCTCATCTGTAAATACTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAATGCCACTGTCACGCTGATCACCAGGTACACCACCGCGGCGGCGgtg                                                                                                                                                <  2:189774/149‑1 (MQ=255)
      aaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtg                                                                                                                                                <  1:14236/149‑1 (MQ=255)
         cccTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtggtg                                                                                                                                             <  4:190850/149‑4 (MQ=255)
                    cTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGAc                                                                                                                                           <  6:129283/140‑1 (MQ=255)
                    cTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGAc                                                                                                                                           >  5:129283/1‑140 (MQ=255)
                    cTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGcc                                                                                                                                  >  1:298782/1‑149 (MQ=255)
                       cGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGccc                                                                                                                                 <  6:29379/147‑1 (MQ=255)
                       cGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCt                                                                                                                                >  1:244403/1‑148 (MQ=255)
                       cGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCt                                                                                                                                <  1:262119/148‑1 (MQ=255)
                                ttGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtggt                                                                                                                                              >  8:113954/1‑123 (MQ=255)
                                      ccACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCtgtg                                                                                                                 >  3:93660/1‑148 (MQ=255)
                                       cacaTCGTCAGCAGTTCGCTCATTGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGa                                                                                                                >  3:253099/1‑148 (MQ=255)
                                        acaTCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTg                                                                                                                               <  2:68796/132‑1 (MQ=255)
                                        acaTCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTg                                                                                                                               >  1:68796/1‑132 (MQ=255)
                                        acaTCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAgg                                                                                                                          >  3:11142/1‑137 (MQ=255)
                                        acaTCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAgg                                                                                                                          <  4:11142/137‑1 (MQ=255)
                                        acaTCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGaaa                                                                                                              >  4:294967/1‑149 (MQ=255)
                                         caTCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAAc                                                                                                             >  7:89586/1‑149 (MQ=255)
                                                  cagTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCTCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGc                                                                                                    <  7:270747/149‑1 (MQ=255)
                                                    gTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCt                                                                                                  <  2:244403/149‑1 (MQ=255)
                                                    gTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCt                                                                                                  >  8:107735/1‑149 (MQ=255)
                                                        gCTCTTCGGGCAATCCTCTGTTTCAACTTGGCCGACAGCCAACGAGAGTACTAGGTCAATGGCACGCTGAGCACCATGTAGACGACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTg                                                                                                                      <  1:227825/125‑1 (MQ=255)
                                                        gCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGaaa                                                                                              <  5:282308/149‑1 (MQ=255)
                                                        gCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGaa                                                                                               <  1:129896/148‑1 (MQ=255)
                                                        gCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGa                                                                                                >  8:255039/1‑147 (MQ=255)
                                                           cATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCt                                                                                           <  3:294967/149‑1 (MQ=255)
                                                               gggCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGAt                                                                                        <  8:89586/148‑1 (MQ=255)
                                                                 gCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTc                                                                                     <  6:51372/149‑1 (MQ=255)
                                                                        tctGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTCCACCAACg                                                                              >  8:271084/1‑149 (MQ=255)
                                                                               cAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTAAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAAtcg                                                                       >  4:84097/1‑149 (MQ=255)
                                                                                    tGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCt                                                                                           <  8:218481/124‑1 (MQ=255)
                                                                                    tGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCt                                                                                           >  7:218481/1‑124 (MQ=255)
                                                                                             gCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGc                                                                >  1:303491/1‑142 (MQ=255)
                                                                                             gCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGc                                                                <  2:303491/142‑1 (MQ=255)
                                                                                              ccAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGcc                                                         >  3:149029/1‑148 (MQ=255)
                                                                                                aaCGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAg                                                      >  1:43904/1‑149 (MQ=255)
                                                                                                aaCGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAg                                                      >  6:70512/1‑149 (MQ=255)
                                                                                                 aCGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGt                                                     >  2:176398/1‑149 (MQ=255)
                                                                                                  cGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTc                                                    >  2:157045/1‑149 (MQ=255)
                                                                                                  cGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTc                                                    <  7:107735/149‑1 (MQ=255)
                                                                                                  cGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTCCACCAACGCCAATCGTCGAGGCCACGGCCGAGTc                                                    <  7:271084/149‑1 (MQ=255)
                                                                                                     gagaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCtt                                                 <  1:176398/149‑1 (MQ=255)
                                                                                                        aaCAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAg                                              >  7:8136/1‑149 (MQ=255)
                                                                                                          cAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAg                                              >  7:79439/1‑147 (MQ=255)
                                                                                                               ccACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGcc                                        >  6:65211/1‑148 (MQ=255)
                                                                                                                  cTGGCACGCTAAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATg                                    <  3:84097/149‑1 (MQ=255)
                                                                                                                     gCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCg                                 <  6:133088/149‑1 (MQ=255)
                                                                                                                       aCGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTa                               <  7:255039/149‑1 (MQ=255)
                                                                                                                        cGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAg                              >  2:12320/1‑149 (MQ=255)
                                                                                                                          cTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGt                            <  5:65211/149‑1 (MQ=255)
                                                                                                                                 ccAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGc                                  >  6:240812/1‑136 (MQ=255)
                                                                                                                                 ccAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGc                                  <  5:240812/136‑1 (MQ=255)
                                                                                                                                 ccAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAg                     <  4:149029/149‑1 (MQ=255)
                                                                                                                                  cAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATg                                    >  1:79000/1‑133 (MQ=255)
                                                                                                                                  cAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATg                                    <  2:79000/133‑1 (MQ=255)
                                                                                                                                    ggTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTg                  <  4:93660/149‑1 (MQ=255)
                                                                                                                                     gTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTgg                 >  6:135755/1‑149 (MQ=255)
                                                                                                                                      gtctacaccGCGGCGGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTgg                 >  7:183002/6‑148 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg           <  1:12320/149‑1 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg           >  8:127464/1‑149 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg           <  5:70512/149‑1 (MQ=255)
                                                                                                                                             ccgccgccgCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCg         <  2:43904/149‑1 (MQ=255)
                                                                                                                                              cgccgccgcGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGa        <  7:127464/149‑1 (MQ=255)
                                                                                                                                                    cgcgGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCa  >  2:263775/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCA  >  NC_002947/5405908‑5406204

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: