Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,461,502 Δ1 bp 5.1% intergenic (+17/‑11) PP_3954 → / → PP_3955 substrate‑binding protein/permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,461,5020T. 5.1% 251.2 / 5.4 79intergenic (+17/‑11)PP_3954/PP_3955substrate‑binding protein/permease
Reads supporting (aligned to +/- strand):  ref base T (40/33);  new base . (2/2);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTTCCACGCCTGGTCACGGCCCACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCACCCAAGTACTGCGTGC  >  NC_002947/4461355‑4461649
                                                                                                                                                   |                                                                                                                                                   
aCCTTCCACGCCTGGTCACGGTCCACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATg                                                                                                                                                    <  8:41391/149‑1 (MQ=255)
  cTTCCACGCCTGGTCACGGCCCACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATg                                                                                                                                                    >  3:199152/1‑147 (MQ=255)
      cACGCCTGGTCACGGCCCACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCCTCCAGGTGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGagcag                                                                                                                                              >  3:236115/1‑149 (MQ=255)
       aCGCCTGGTCACGGCCCACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGGGCAGCGAGGTACTTCAAGGAGCAAAACCTGGTCCCCTGAGGACGTTCCTGAGCCATGGagcagc                                                                                                                                             >  8:314898/1‑149 (MQ=255)
       aCGCCTGGTCACGGCCCACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGagcagc                                                                                                                                             <  6:188282/149‑1 (MQ=255)
           cTGGTCACGGCCCACGCTTCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCCTAt                                                                                                                                         <  3:329739/149‑1 (MQ=255)
           cTGGTCACGGCCCACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGTCATGGAGCGGCGTAt                                                                                                                                         <  2:8639/149‑1 (MQ=255)
           cTGGTCACGGCCCACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGCGCAGCGTAt                                                                                                                                         >  7:71358/1‑149 (MQ=255)
                 aCGGCCCACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGa                                                                                                                                   >  4:165259/1‑149 (MQ=255)
                    gCCCACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGTCCTGGAGCAGCGTATGAGTGATAc                                                                                                                                <  7:55334/149‑1 (MQ=255)
                     cccACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATAcc                                                                                                                               >  3:285781/1‑149 (MQ=255)
                      ccACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACcg                                                                                                                              >  2:334301/1‑149 (MQ=255)
                       cACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACcgc                                                                                                                             <  4:285781/149‑1 (MQ=255)
                            tGCAGCGAAAGAGATGCAGCGAGAAATTGGCGCCAAGAACCTGGCGAATCCGCTGCATCCAGGAGCAGCGAGGTACTTCAACGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAg                                                                                                                        <  4:119628/149‑1 (MQ=255)
                                       atatCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGc                                                                                                              >  2:162242/1‑148 (MQ=255)
                                          tCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTTCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAg                                                                                                           <  7:236727/148‑1 (MQ=255)
                                              ggCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCAcc                                                                                                       >  4:326681/1‑148 (MQ=255)
                                                cAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGa                                                                                                                                       <  7:61795/114‑1 (MQ=255)
                                                cAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGa                                                                                                                                       >  8:61795/1‑114 (MQ=255)
                                                cAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGAt                                                                                                    >  8:287817/1‑149 (MQ=255)
                                                   aaaaTGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCAGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCt                                                                                                >  3:316778/1‑149 (MQ=255)
                                                   aaaaTGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCAGGAGCAGCGTATGAGTGATACAGAGCAAGGTCTTCACGCCAGAACCGATCGCt                                                                                                >  3:322176/1‑149 (MQ=255)
                                                     aaTGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTg                                                                                               <  5:167164/149‑1 (MQ=255)
                                                               aGAAGCTGCCCATTCCGCTTCATCCAGGAGCAGCGAGGTACTTCAATGAGAAAAACCTGCTCCCCTGATGACATTCATGAGCCAGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCgcgc                                                                                    <  4:316778/149‑1 (MQ=255)
                                                               aGAACCTGTCCATTCCGCTTCATACAGGAGCAGCGAGGTACTTCAATGAGAAAAACCTGCTCCCCTGATGACATTCATGAGCCAGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCgcgc                                                                                    <  4:322176/149‑1 (MQ=255)
                                                                gAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCAGGGAGCAGCGTATGAGTGATACCGCGCAAGGTCGTCACGCCAGCACCGATCGCGGGCCGGCCtcgc                                                                                    >  4:231805/1‑149 (MQ=255)
                                                                  aCCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGTTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCtt                                                                                  >  6:219640/1‑149 (MQ=255)
                                                                   ccTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCttt                                                                                 >  3:352281/1‑149 (MQ=255)
                                                                       cccATTCCGCTTCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAATGTCTTCACGCCAACACCCATCGCTGGCCTTCCGCGCtttttt                                                                              >  4:263410/1‑147 (MQ=255)
                                                                       cccATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTg                                                                              <  4:199152/148‑1 (MQ=255)
                                                                        ccATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGcc                                                                            >  3:81865/1‑149 (MQ=255)
                                                                          aTTCCGCTGCATCCAGGAGCAGCGAGGTACTTCATGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGt                                                                          <  3:264291/149‑1 (MQ=255)
                                                                               gCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCAt                                                                     >  2:28712/1‑149 (MQ=255)
                                                                                cTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCAtt                                                                    <  3:96920/149‑1 (MQ=255)
                                                                                   cATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTg                                                                   <  4:81865/147‑1 (MQ=255)
                                                                                    aTCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTc                                                                >  6:35497/1‑149 (MQ=255)
                                                                                     tCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCt                                                               >  7:226003/1‑149 (MQ=255)
                                                                                        aGGAGCAGCGAGGTTCTTCTAGGATAAAAAGCTGCTCCCCTGATGACATTCCTGTGCCATGGAGCTGCGTATGAGTGATACCGCGCAAGGTCTTCACGGCAGCACCGTTCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGTTCtttt                                                            <  3:231805/149‑1 (MQ=255)
                                                                                        aGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCttt                                                             >  2:366027/1‑148 (MQ=255)
                                                                                          gAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCg                                                          >  2:286992/1‑149 (MQ=255)
                                                                                           agcagcGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCgg                                                         <  7:287817/149‑1 (MQ=255)
                                                                                            gcagcGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGt                                                        >  8:105534/1‑149 (MQ=255)
                                                                                               gcGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGtttt                                                     >  1:221722/1‑149 (MQ=255)
                                                                                                cGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTc                                                    <  4:236115/149‑1 (MQ=255)
                                                                                                  aGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCa                                                  <  3:165259/149‑1 (MQ=255)
                                                                                                        ttCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCa                                                  >  2:189896/1‑143 (MQ=255)
                                                                                                        ttCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCa                                                  <  1:189896/143‑1 (MQ=255)
                                                                                                           aaGGAGAAAAACCTGCTCCCCTGATGACATACCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCAc                                         >  1:33796/1‑149 (MQ=255)
                                                                                                             ggAGAAAAACCTGCTCCCCTGATGACATACCTGAGCCAAGGAGCAGCGAATGAGTCATACCGCGCAAGGTCTTCACGACACCACCCATCCCAGCCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCAcgg                                       >  7:231934/1‑149 (MQ=255)
                                                                                                               agaAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCAcggcg                                     <  4:352281/149‑1 (MQ=255)
                                                                                                                gaAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCAcggcgg                                    >  1:230966/1‑149 (MQ=255)
                                                                                                                 ataaagatGCTCCCATTATGACATTCCTGAGCCATGGAGCAGGGTATGAGTGTTAGCGGGGAAGGTCTTCACGCCATCACGGATGGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCAcggcgg                                    <  3:263410/141‑1 (MQ=255)
                                                                                                                 aaaaaCCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACggtggc                                   >  1:114272/1‑149 (MQ=255)
                                                                                                                      ccTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTcc                              <  2:221722/149‑1 (MQ=255)
                                                                                                                       cTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTccc                             >  5:233137/1‑149 (MQ=255)
                                                                                                                          cTCCCCTGCTGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGGGGCATTGGGCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCt                          <  7:314898/149‑1 (MQ=255)
                                                                                                                          cTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTcccc                            <  6:272341/147‑1 (MQ=255)
                                                                                                                             cccTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGAt                       <  1:334301/149‑1 (MQ=255)
                                                                                                                             cccTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGAt                       >  2:111883/1‑149 (MQ=255)
                                                                                                                             cccTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGAt                       >  7:160547/1‑149 (MQ=255)
                                                                                                                               cTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGAttt                     >  3:291351/1‑149 (MQ=255)
                                                                                                                                tgatgaCATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCt                                 >  2:291042/1‑136 (MQ=255)
                                                                                                                                tgatgaCATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCt                                 <  1:291042/136‑1 (MQ=255)
                                                                                                                                tgatgaCATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTc                    >  6:151864/1‑149 (MQ=255)
                                                                                                                                 gatgaCATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGt                                                                          <  2:217079/94‑1 (MQ=255)
                                                                                                                                 gatgaCATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGt                                                                          >  1:217079/1‑94 (MQ=255)
                                                                                                                                   tgaCATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCAc                 <  2:33796/149‑1 (MQ=255)
                                                                                                                                     aCATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCAccc               <  2:230966/149‑1 (MQ=255)
                                                                                                                                       aTTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGTTTTTACATATCGCCACGGCGGCCTTCCCCCTGATTTCCAccc               >  8:147666/1‑147 (MQ=255)
                                                                                                                                            tGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCAAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCACCCAAGTACt        >  7:111803/1‑149 (MQ=255)
                                                                                                                                               gCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCACCCAAGTACTGCg     <  8:226003/149‑1 (MQ=255)
                                                                                                                                                ccATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCACCCAAGTACTGCGt    >  6:248169/1‑149 (MQ=255)
                                                                                                                                                ccATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCACCCAAGTACTGCGt    <  5:35497/149‑1 (MQ=255)
                                                                                                                                                 cATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTTGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCACCCAAGTACTGCGTg   <  1:111883/149‑1 (MQ=255)
                                                                                                                                                  aTGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCACCCAAGTACTGCGTg   >  3:341981/1‑148 (MQ=255)
                                                                                                                                                  aTGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCACCCAAGTACTGCGTg   >  3:347418/1‑148 (MQ=255)
                                                                                                                                                  aTGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCACCCAAGTACTGCGTg   >  3:262633/1‑148 (MQ=255)
                                                                                                                                                  aTGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCACCCAAGTACTGCGTGc  <  4:291351/149‑1 (MQ=255)
                                                                                                                                                  aTGGAGCAGCGTATGAGTGAAACCGCGCAAGGTCTTCACGACAGCAACGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCACCCAAGTACTGCGTg   >  6:55326/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCTTCCACGCCTGGTCACGGCCCACGCTGCAGCGAAAGATATCCAGGCAGAAAATGCCGCCAAGAACCTGCCCATTCCGCTGCATCCAGGAGCAGCGAGGTACTTCAAGGAGAAAAACCTGCTCCCCTGATGACATTCCTGAGCCATGGAGCAGCGTATGAGTGATACCGCGCAAGGTCTTCACGCCAGCACCGATCGCTGGCCTGCCGCGCTTTTTGCCGTGGCATTGGTCTTTTCGGTTTTCCAGATCGCCACGGCGGCCTTCCCCCTGATTTCCACCCAAGTACTGCGTGC  >  NC_002947/4461355‑4461649

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: