Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 558,477 A→C 6.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 6.3% 177.5 / 4.3 64E21D (GAA→GAC) rplR50S ribosomal protein L18
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base C (2/2);  total (33/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.87e-01

GCCAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCG  >  NC_002947/558329‑558621
                                                                                                                                                    |                                                                                                                                                
gCCAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGaa                                                                                                                                                  <  2:5076/149‑1 (MQ=255)
       ccTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAg                                                                                                                                           >  1:208943/1‑149 (MQ=255)
          tACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGaa                                                                                                                                            <  5:140311/145‑1 (MQ=255)
          tACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGaa                                                                                                                                            >  6:140311/1‑145 (MQ=255)
            cAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGt                                                                                                                                       <  4:73618/148‑1 (MQ=255)
            cAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGt                                                                                                                                       <  8:293323/148‑1 (MQ=255)
             aaaggtaaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCAAGTCTCAAGATGCACGAACTCTAACTCGTgc                                                                                                                                     >  1:233361/1‑149 (MQ=255)
               aggtaaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTgcgc                                                                                                                                   <  2:208943/149‑1 (MQ=255)
                       gtgtGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCgt                                                                                                                           >  5:20133/1‑149 (MQ=255)
                          tgCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGGCAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGtt                                                                                                                        <  8:250194/149‑1 (MQ=255)
                              ttACGCTTACGAAGTATTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTTAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGc                                                                                                                    <  3:57160/149‑1 (MQ=255)
                                aCGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGctc                                                                                                                  >  6:15884/1‑149 (MQ=255)
                                  gcgGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGctcct                                                                                                                <  6:20133/149‑1 (MQ=255)
                                          agtagtCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCgcagca                                                                                                         >  8:26673/1‑148 (MQ=255)
                                           gtagtCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGccca                                                                                                       >  3:229531/1‑149 (MQ=255)
                                            tagtCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACAt                                                                                                      <  7:26673/149‑1 (MQ=255)
                                              gtCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATc                                                                                                     >  3:289059/1‑148 (MQ=255)
                                                 cgtcgtAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGTGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACg                                                                                                 >  1:301342/1‑149 (MQ=255)
                                                 cgtcgtAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTAc                                                                                                  <  5:119792/148‑1 (MQ=255)
                                                  gtcgtAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGc                                                                                                >  1:118073/1‑149 (MQ=255)
                                                     gtAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCCCGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGcc                                                                                               >  5:194582/1‑147 (MQ=255)
                                                      tAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAg                                                                                            >  7:77820/1‑149 (MQ=255)
                                                      tAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTATTCCGCTCCTCGCAGCACATCTACGCCCAg                                                                                            <  1:25995/149‑1 (MQ=255)
                                                        aagaagCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGt                                                                                          >  6:145890/1‑149 (MQ=255)
                                                        aagaagCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGt                                                                                          >  6:61522/1‑149 (MQ=255)
                                                          gaagCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCa                                                                                        >  7:210453/1‑149 (MQ=255)
                                                           aagCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAg                                                                                                                                           >  4:9747/1‑97 (MQ=255)
                                                           aagCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAg                                                                                                                                           <  3:9747/97‑1 (MQ=255)
                                                              ccAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTc                                                                                    <  6:194582/149‑1 (MQ=255)
                                                                      aagTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACg                                                                            >  6:249457/1‑149 (MQ=255)
                                                                      aagTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACg                                                                            >  3:165388/1‑149 (MQ=255)
                                                                        gTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACgg                                                                           >  6:62057/1‑148 (MQ=255)
                                                                          aGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACgg                                                                           <  6:46917/146‑1 (MQ=255)
                                                                          aGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACgg                                                                           >  5:46917/1‑146 (MQ=255)
                                                                               cTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAgg                                                                   >  5:74424/1‑149 (MQ=255)
                                                                                     aaTGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTg                                                              <  5:249457/148‑1 (MQ=255)
                                                                                     aaTGACAGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGATCTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCtt                                                              <  4:165388/148‑2 (MQ=255)
                                                                                             aCAAAAAAGTTAATCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTACTGACAAgcgc                                                      >  4:166246/1‑148 (MQ=255)
                                                                                                    aGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCt                                              <  5:145890/149‑1 (MQ=255)
                                                                                                      ttATTCGACTGCGTCTCGCTCGCAAAGCACGTCTCAATATGCACGAACTCTAAGTCGTGCGCCTGTGCGTGGTCCGCTCCTCGCCTCACATCTACGCACAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGAAAGCGCCTCGACCTTg                                            <  5:193931/149‑1 (MQ=255)
                                                                                                      ttATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCACGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTg                                            <  4:229531/149‑1 (MQ=255)
                                                                                                      ttATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTg                                            <  5:189474/149‑1 (MQ=255)
                                                                                                         ttCGACTGCGTCGCGCTCGCATAGCAGGTCTCTAGAAGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGGTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAATGGTCTGGCAAGCGGCTCGACCTTGGa                                          <  7:88737/148‑1 (MQ=255)
                                                                                                            gACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCAACCTTGGACa                                        >  8:95610/1‑147 (MQ=255)
                                                                                                                 cGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCt                                 <  5:62057/149‑1 (MQ=255)
                                                                                                                  gTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCACTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTg                                >  1:224165/1‑149 (MQ=255)
                                                                                                                   tCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGc                               <  5:61522/149‑1 (MQ=255)
                                                                                                                     gcgcTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGt                             >  6:168967/1‑149 (MQ=255)
                                                                                                                     gcgcTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGt                             <  8:77820/149‑1 (MQ=255)
                                                                                                                     gcgcTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCAACCTTGGACAAAGACCTGCGt                             <  7:95610/149‑1 (MQ=255)
                                                                                                                      cgcTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGACGTGCGCCTGTGCGTGTTCCGCTCCTCGCACCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAACCGCCTCGACCTTGGACAAAGACCTGCGt                             >  3:279606/1‑148 (MQ=255)
                                                                                                                           gCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCACTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATg                         <  2:224165/147‑1 (MQ=255)
                                                                                                                                  aCGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTg                <  5:15884/149‑1 (MQ=255)
                                                                                                                                      ctcAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCa             >  6:203974/1‑148 (MQ=255)
                                                                                                                                       tcAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCaa            >  5:141117/1‑148 (MQ=255)
                                                                                                                                         aaGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACAt         >  3:15019/1‑149 (MQ=255)
                                                                                                                                          aggcgCACGACTTCGAGTTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATc        >  6:109905/5‑149 (MQ=255)
                                                                                                                                           ggcgCACGACTTCGAGTTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTg                                                              >  4:298288/4‑94 (MQ=255)
                                                                                                                                           ggcgCACGACTTCGAGTTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTg                                                              <  3:298288/91‑1 (MQ=255)
                                                                                                                                           gATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCg       <  5:203974/149‑1 (MQ=255)
                                                                                                                                             cgCACGACTTCGAGTTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGAc     <  5:109905/148‑1 (MQ=255)
                                                                                                                                              gaaCGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGTCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGTTGCCAGTGGCGCCATcac      >  4:190334/3‑145 (MQ=255)
                                                                                                                                              gCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGTCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGTTGCCAGTGGCGCCATcac      <  3:190334/147‑3 (MQ=255)
                                                                                                                                               cACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGAcgc   <  6:74424/149‑1 (MQ=255)
                                                                                                                                                aCGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACGTCGAcgcg  >  3:241623/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCCAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCG  >  NC_002947/558329‑558621

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: ‑