Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,709,377 G→C 6.8% R1132P (CGG→CCG) ‡ PP_5009 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,709,3770GC 6.8% 186.2 / 3.2 59R1132P (CGG→CCG) ‡PP_5009hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base C (2/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

CGGTGCAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGA  >  NC_002947/5709231‑5709520
                                                                                                                                                  |                                                                                                                                               
cGGTGCAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACgg                                                                                                                                                <  4:75530/148‑1 (MQ=255)
    gCAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACgg                                                                                                                                                <  7:165246/144‑1 (MQ=255)
    gCAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACgg                                                                                                                                                >  8:165246/1‑144 (MQ=255)
          aTCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGGACCCGTTCGCGAtaa                                                                                                                                     <  7:159518/149‑3 (MQ=255)
           tCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGATCGGGTCCGCGATTcc                                                                                                                                    <  8:184937/149‑1 (MQ=255)
              aaaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagc                                                                                                                                 >  7:281681/1‑149 (MQ=255)
              aaaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagc                                                                                                                                 >  5:207101/1‑149 (MQ=255)
              aaaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGACTACCTCCACGCGCTGAACGCTCAGTCCATCTGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACACCAATACCTCGATTATCGGGCAGGGGTCCGGGATTccagc                                                                                                                                 >  8:48646/1‑149 (MQ=255)
               aaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagc                                                                                                                                 >  1:113538/1‑148 (MQ=255)
               aaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagc                                                                                                                                 >  2:195621/1‑148 (MQ=255)
                 aCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagccag                                                                                                                              <  4:210346/149‑1 (MQ=255)
                 aCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagccag                                                                                                                              >  7:249528/1‑149 (MQ=255)
                   gACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCTgccaggc                                                                                                                            <  8:15553/149‑1 (MQ=255)
                      aTGCTGCCGTGCCGCCGCGTGCTGGAGCCGGCCGACCACCCCGCGCTCAGCGCACAGGCCATCGGCCTGCCCCGCGGTCCAGCGACTGCCGACCAGATACAGCACTACCTCGATTATCGCGAACGGGTCCGCGATTCCAgccaggccag                                                                                                                         <  8:124954/149‑1 (MQ=255)
                              gtgCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGAGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTTTAGCGAACGGGTCCGCGATTCCAGCCAGGCCAgccgccgc                                                                                                                  >  7:195059/1‑148 (MQ=255)
                               tgCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAgccgccgcc                                                                                                                 <  4:23880/148‑1 (MQ=255)
                                 ccgccgCTTGCTGGAGCAGGTCGACCAGCAGGCGCTGAGCGCTCAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTCCCTAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCTCGATTCCAGCCAGGCCAGCCGCCGCCTGc                                                                                                              <  7:48646/149‑1 (MQ=255)
                                 ccgccgCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCcgc                                                                                                                                          <  8:10147/121‑1 (MQ=255)
                                 ccgccgCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCcgc                                                                                                                                          >  7:10147/1‑121 (MQ=255)
                                 ccgccgCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTg                                                                                                               <  5:212295/148‑1 (MQ=255)
                                      gcTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTTTCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTgtc                                                                                                         >  7:230011/1‑149 (MQ=255)
                                                 cAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCg                                                                                              <  1:138413/149‑1 (MQ=255)
                                                    gTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCt                                                                                           >  8:243532/1‑149 (MQ=255)
                                                       gACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTgga                                                                                        >  8:84571/1‑149 (MQ=255)
                                                          caccaGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTggatcg                                                                                     >  2:217035/1‑149 (MQ=255)
                                                            ccaGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTgg                                                                                         <  4:191624/143‑1 (MQ=255)
                                                            ccaGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTgg                                                                                         >  3:191624/1‑143 (MQ=255)
                                                                  gcTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTgg                                                                             >  4:219974/1‑149 (MQ=255)
                                                                   cTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGc                                                                            <  3:173772/149‑1 (MQ=255)
                                                                   cTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGc                                                                            <  6:207101/149‑1 (MQ=255)
                                                                    tGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGCAAGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCGTGGCGTCAGTCAGCGGCTGGATCGGCTCCTGGCg                                                                           >  4:204107/1‑149 (MQ=255)
                                                                              gTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGc                                                                 <  5:127406/149‑1 (MQ=255)
                                                                                           ccccGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGa                                                     <  2:113538/148‑1 (MQ=255)
                                                                                           ccccGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGAt                                                    <  1:261993/149‑1 (MQ=255)
                                                                                           ccccGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGAt                                                    <  1:258649/149‑1 (MQ=255)
                                                                                            cccGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATc                                                   <  7:84571/149‑1 (MQ=255)
                                                                                             ccGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAg                                                               <  5:227850/136‑1 (MQ=255)
                                                                                             ccGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAg                                                               >  6:227850/1‑136 (MQ=255)
                                                                                                    ccAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTcc                                           >  6:83684/1‑149 (MQ=255)
                                                                                                     cAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCATGACCCCAAGATCGATTTTcc                                           >  2:89117/1‑148 (MQ=255)
                                                                                                     cAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTcc                                           <  5:83684/148‑1 (MQ=255)
                                                                                                            tGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAg                                   >  3:168175/1‑149 (MQ=255)
                                                                                                                    aGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGaaa                           >  1:236405/1‑149 (MQ=255)
                                                                                                                       tACAGCACTACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTa                        <  8:281681/149‑1 (MQ=255)
                                                                                                                         cAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAg                       <  1:217035/148‑1 (MQ=255)
                                                                                                                          aGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGcc                                 >  5:182694/1‑137 (MQ=255)
                                                                                                                          aGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGcc                                 <  6:182694/137‑1 (MQ=255)
                                                                                                                           gCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAg                       <  2:251058/146‑1 (MQ=255)
                                                                                                                           gCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAg                       >  1:251058/1‑146 (MQ=255)
                                                                                                                                  cTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCAc                               >  7:40409/1‑131 (MQ=255)
                                                                                                                                  cTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCAc                               <  8:40409/131‑1 (MQ=255)
                                                                                                                                  cTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAt             >  8:96323/1‑149 (MQ=255)
                                                                                                                                     tggaaTCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTa                        <  1:240968/131‑1 (MQ=255)
                                                                                                                                     tggaaTCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTa                        >  2:240968/5‑135 (MQ=255)
                                                                                                                                     tggaaTCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTg           >  4:16191/5‑148 (MQ=255)
                                                                                                                                     gATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAtt            <  7:38016/147‑1 (MQ=255)
                                                                                                                                        tATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACg       >  8:79438/1‑149 (MQ=255)
                                                                                                                                          tCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACg       >  2:111016/1‑147 (MQ=255)
                                                                                                                                            gcgAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAgg   <  7:3060/149‑1 (MQ=255)
                                                                                                                                             cgAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCATGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGa  <  1:89117/149‑1 (MQ=255)
                                                                                                                                             cgAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAAc        >  7:54407/1‑143 (MQ=255)
                                                                                                                                             cgAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAAc        <  8:54407/143‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CGGTGCAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGA  >  NC_002947/5709231‑5709520

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