Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 558,477 A→C 5.3% E21D (GAA→GAC rplR → 50S ribosomal protein L18

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947558,4770AC 5.3% 226.0 / 3.5 76E21D (GAA→GACrplR50S ribosomal protein L18
Reads supporting (aligned to +/- strand):  ref base A (33/39);  new base C (2/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

CCAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGCCA  >  NC_002947/558330‑558625
                                                                                                                                                   |                                                                                                                                                    
ccAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAAc                                                                                                                                                     <  3:20649/149‑1 (MQ=255)
ccAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAAc                                                                                                                                                     <  4:214042/149‑1 (MQ=255)
     gCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGaa                                                                                                                                                <  6:199005/149‑1 (MQ=255)
      ccTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAg                                                                                                                                               <  7:250263/149‑1 (MQ=255)
        ttACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTc                                                                                                                                             >  7:337997/1‑149 (MQ=255)
           cAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTg                                                                                                                                          >  2:319375/1‑149 (MQ=255)
            aaaggtaaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTgc                                                                                                                                         >  4:63525/1‑149 (MQ=255)
                gtaaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGAGTTCGTgcatc                                                                                                                                      >  6:100673/1‑145 (MQ=255)
                gtaaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGAGTTCGTgcatc                                                                                                                                      <  5:100673/148‑4 (MQ=255)
                gtaaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCg                                                                                                                                                  >  6:245304/1‑136 (MQ=255)
                gtaaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCg                                                                                                                                                  <  5:245304/136‑1 (MQ=255)
                      gtgtGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCgt                                                                                                                               <  8:51844/149‑1 (MQ=255)
                          gCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGtt                                                                                                                            <  1:319375/148‑1 (MQ=255)
                           cGTTACGCGGACGAAGTAGTCCGTGGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTcc                                                                                                                          >  1:14829/1‑149 (MQ=255)
                                 gcgGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGAGTTCGTgcat                                                                                                                                       <  2:281119/130‑3 (MQ=255)
                                 gcgGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGAGTTCGTgcat                                                                                                                                       >  1:281119/1‑128 (MQ=255)
                                 gcgGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGctcct                                                                                                                    <  7:283106/149‑1 (MQ=255)
                                 gcgGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGctcct                                                                                                                    >  7:77796/1‑149 (MQ=255)
                                  cgGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGctcctc                                                                                                                   <  4:35019/149‑1 (MQ=255)
                                     aCGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCgca                                                                                                                <  7:170836/149‑1 (MQ=255)
                                     aCGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCgcag                                                                                                               >  8:243079/1‑149 (MQ=255)
                                      cGAAGTAGTCCGTCGTAAAGACGCCAAGAAGAAGTTGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTGTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCgcag                                                                                                               <  2:14829/149‑1 (MQ=255)
                                         agtagtCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGcac                                                                                                            <  5:232709/149‑1 (MQ=255)
                                          gtagtCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGcaca                                                                                                           <  1:158305/149‑1 (MQ=255)
                                             gtCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATc                                                                                                         >  3:190482/1‑148 (MQ=255)
                                                cgtcgtAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTAc                                                                                                      <  8:339189/148‑1 (MQ=255)
                                                     tAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAg                                                                                                >  4:141384/1‑149 (MQ=255)
                                                              cAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATc                                                                                                         >  2:11782/1‑131 (MQ=255)
                                                              cAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATc                                                                                                         <  1:11782/131‑1 (MQ=255)
                                                               aagaagaagTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCgg                                                                                      <  3:63525/149‑1 (MQ=255)
                                                                    gaagTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTc                                                                                                                                             >  4:114087/1‑89 (MQ=255)
                                                                    gaagTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTc                                                                                                                                             <  3:114087/89‑1 (MQ=255)
                                                                          gggCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTgc                                                                                                                                         <  5:108546/87‑1 (MQ=255)
                                                                          gggCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTgc                                                                                                                                         >  6:108546/1‑87 (MQ=255)
                                                                          gggCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGgcag                                                                            >  8:318572/1‑148 (MQ=255)
                                                                          gggCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGgcag                                                                            <  7:100440/148‑1 (MQ=255)
                                                                          gggCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGgcag                                                                            <  3:275289/148‑1 (MQ=255)
                                                                          gggCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGgcag                                                                            >  8:100440/1‑148 (MQ=255)
                                                                          gggCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCACCTCGCAGCACATCTACCCCCAGGTCATTACGGCCGACGgcag                                                                            >  4:275289/1‑148 (MQ=255)
                                                                             ccTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAg                                                                        >  3:204451/1‑149 (MQ=255)
                                                                              cTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAgg                                                                       >  5:69703/1‑149 (MQ=255)
                                                                                    aaTGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTgg                                                                 <  7:318572/149‑1 (MQ=255)
                                                                                    aaTGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTgg                                                                 >  4:175404/1‑149 (MQ=255)
                                                                                    aaTGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTgg                                                                 <  7:246072/149‑1 (MQ=255)
                                                                                    aaTGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTg                                                                  <  8:51587/148‑1 (MQ=255)
                                                                                     aTGACCTACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGc                                                                <  6:341201/149‑1 (MQ=255)
                                                                                          cGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAgcg                                                           >  7:286398/1‑149 (MQ=255)
                                                                                          cGACAAAAAAGTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAgcgc                                                          <  7:243079/149‑1 (MQ=255)
                                                                                           gACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAgcgc                                                          >  4:8497/1‑149 (MQ=255)
                                                                                              aaaaaaGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTc                                                       <  4:204451/149‑1 (MQ=255)
                                                                                               aaaaaGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCg                                                      <  5:94265/149‑1 (MQ=255)
                                                                                                aaaaGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGa                                                     >  8:45663/1‑149 (MQ=255)
                                                                                                 aaaGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCAGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTACTGGCAACCGCCTAGa                                                     >  7:122570/1‑148 (MQ=255)
                                                                                                    gTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCtt                                                 >  5:175697/1‑149 (MQ=255)
                                                                                                       aTTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGa                                              >  4:92960/1‑149 (MQ=255)
                                                                                                            aCTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCtt                                                 >  8:340298/1‑141 (MQ=255)
                                                                                                            aCTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCtt                                                 <  7:340298/141‑1 (MQ=255)
                                                                                                                 gTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTg                                    >  6:282034/1‑149 (MQ=255)
                                                                                                                  tCGCGCTCTCAAAGGACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGc                                   <  8:73545/149‑1 (MQ=255)
                                                                                                                  tCGCGCTCGGATAGCACGTCTCAAGATTCAAGAACTCTAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGGACATCTACTCCCAGGTCATTTCGGCCGACGTCAGCAAGGTTCTGGCAAGCGCCTAGACCTTGGACAAAGACCTGc                                   <  8:122570/149‑1 (MQ=255)
                                                                                                                     cgcTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTg                                <  3:92960/149‑1 (MQ=255)
                                                                                                                      gcTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGa                               <  4:190482/149‑1 (MQ=255)
                                                                                                                            aaaGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCACGCGCCTCGACCTTGGACAAAGACCTGCGTGATGgcgc                         >  1:303571/1‑149 (MQ=255)
                                                                                                                            aaaGAACGTCTCAAGATGCACGAACTAGAAGTCGTGCGCCTGTGCGTGTTACGCTCCTCCAAGCACATCTAAGCCCAGGTCATTTCGGCAGACGGCAGCAAGTTTCAGGCAAGACACTCGACCTTGGACAAATACCTTCGTGATGgcg                          >  3:349767/1‑148 (MQ=255)
                                                                                                                             aaGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGgcgc                         >  4:295295/1‑148 (MQ=255)
                                                                                                                                cACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACt                     <  3:295295/149‑1 (MQ=255)
                                                                                                                                 aCGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTg                    <  6:69703/149‑1 (MQ=255)
                                                                                                                                    tctcAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCgg                                                                                      <  8:249164/80‑1 (MQ=255)
                                                                                                                                    tctcAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCgg                                                                                      >  7:249164/1‑80 (MQ=255)
                                                                                                                                       cAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGc                                   <  2:240715/128‑1 (MQ=255)
                                                                                                                                       cAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGc                                   >  1:240715/1‑128 (MQ=255)
                                                                                                                                             gCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGAcg        >  5:357714/1‑149 (MQ=255)
                                                                                                                                              cACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTg                                                >  2:323687/1‑108 (MQ=255)
                                                                                                                                              cACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTg                                                <  1:323687/108‑1 (MQ=255)
                                                                                                                                                 gAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGc    <  8:77796/149‑1 (MQ=255)
                                                                                                                                                 gAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTAATTTCGGCCGACGTCAGCAAGGTTCTGGCAAGCGCCTCGACCTTTGACAAAGAACTGCGTGATGGCGCCACTGGCAACATCGACGCGGc    <  4:349767/149‑1 (MQ=255)
                                                                                                                                                   aCTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGCCa  <  7:45663/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGCCA  >  NC_002947/558330‑558625

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