Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 935,949 C→T 31.7% E3225K (GAA→AAA)  PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947935,9490CT31.7% 70.4 / 45.4 60E3225K (GAA→AAA) PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (23/18);  new base T (9/10);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.47e-02

AGCCGCTTTGACCCTGGCCGGTAAGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCAT  >  NC_002947/935805‑936092
                                                                                                                                                |                                                                                                                                                
aGCCGCTTTGACCCTGGCCGGTAAGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTgg                                                                                                                                              >  2:292986/1‑149 (MQ=255)
 gCCGCTTTGACCCTGGCCGGTAAGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGc                                                                                                                                             >  4:349398/1‑149 (MQ=255)
 gCCGCTTTGACCCTGGCCGGTAAGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCt                                                                                                                                                <  5:177726/146‑1 (MQ=255)
 gCCGCTTTGACCCTGGCCGGTAAGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCt                                                                                                                                                >  6:177726/1‑146 (MQ=255)
  ccGCTTTGACCCTGGCCGGTAAGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGcc                                                                                                                                            >  7:128065/1‑149 (MQ=255)
    gCTTTGACCCTGGCCGGTAAGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGt                                                                                                                                          >  3:19818/1‑149 (MQ=255)
               ggCCGGTAAGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGAGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCt                                                                                                                               >  7:159921/1‑149 (MQ=255)
                    gTAAGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGtt                                                                                                                          >  7:128289/1‑149 (MQ=255)
                     tAAGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTgg                                                                                                                             <  7:245182/145‑1 (MQ=255)
                     tAAGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTgg                                                                                                                             >  8:245182/1‑145 (MQ=255)
                      aaGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCCCGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGtt                                                                                                                          >  2:54804/1‑147 (MQ=255)
                       aGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAg                                                                                                                       <  6:106315/149‑1 (MQ=255)
                       aGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAg                                                                                                                       <  8:128065/149‑1 (MQ=255)
                            aaTGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGa                                                                                                                     <  2:28788/146‑1 (MQ=255)
                            aaTGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGa                                                                                                                     >  1:28788/1‑146 (MQ=255)
                              tGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACg                                                                                                                <  5:234050/149‑1 (MQ=255)
                                ccAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCt                                                                                                              >  5:244231/1‑149 (MQ=255)
                                 cAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAg                                                                                                                       <  7:61892/139‑1 (MQ=255)
                                 cAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAg                                                                                                                       >  8:61892/1‑139 (MQ=255)
                                 cAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCt                                                                                                              >  1:42294/1‑148 (MQ=255)
                                       tATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATcaca                                                                                                       <  6:46391/149‑1 (MQ=255)
                                           gATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGt                                                                                                   >  5:64817/1‑149 (MQ=255)
                                           gATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACAACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGt                                                                                                   >  8:60524/1‑149 (MQ=255)
                                                 cAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCAc                                                                                             >  5:125355/1‑149 (MQ=255)
                                                 cAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCAc                                                                                             >  3:282216/1‑149 (MQ=255)
                                                     ccGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGAc                                                                                         >  7:304325/1‑149 (MQ=255)
                                                         aGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTgg                                                                                     <  1:343660/149‑1 (MQ=255)
                                                           ggTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGc                                                                                   >  8:207625/1‑149 (MQ=255)
                                                                cAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTg                                                                               <  8:128289/148‑1 (MQ=255)
                                                                  gCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGTGTTGccc                                                                            <  8:304325/149‑1 (MQ=255)
                                                                           cGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑c                                                                  <  4:282216/149‑1 (MQ=255)
                                                                              tATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCg                                                               <  8:159921/149‑1 (MQ=255)
                                                                                 cGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACTACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCg                                                            >  7:180920/1‑149 (MQ=255)
                                                                                   ggCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGgc                                                          <  7:60524/149‑1 (MQ=255)
                                                                                       cgacgaTGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCtg                                                      <  6:244231/149‑1 (MQ=255)
                                                                                        gacgaTGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCtga                                                     >  6:347499/1‑149 (MQ=255)
                                                                                        gacgaTGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCtga                                                     >  3:143782/1‑149 (MQ=255)
                                                                                            aTGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGAc                                                 <  3:349398/149‑1 (MQ=255)
                                                                                            aTGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGAc                                                 <  1:292986/149‑1 (MQ=255)
                                                                                                      tCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACggtg                                       <  5:347499/149‑1 (MQ=255)
                                                                                                         aGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGgtggtg                                    <  4:143782/149‑1 (MQ=255)
                                                                                                          ggTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGc                                   >  5:329965/1‑149 (MQ=255)
                                                                                                             gtTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCg                                >  1:286158/1‑149 (MQ=255)
                                                                                                                gCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGtt                             >  1:177597/1‑149 (MQ=255)
                                                                                                                 ccAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGtt                             <  2:286158/148‑1 (MQ=255)
                                                                                                                 ccAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTc                            >  5:355625/1‑149 (MQ=255)
                                                                                                                  cAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTATGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGCACTGGGCGTTGCCCAGCAAAGT‑CCCCTCGGCGCTGATGACCTTCACGGTGGTGCCCGGttag                           >  6:36799/1‑147 (MQ=255)
                                                                                                                   aGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTcgc                          >  7:231731/1‑149 (MQ=255)
                                                                                                                     cGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTcgccg                        <  8:231731/149‑1 (MQ=255)
                                                                                                                         gTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGAc                    >  3:85735/1‑149 (MQ=255)
                                                                                                                         gTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGAc                    >  2:271763/1‑149 (MQ=255)
                                                                                                                             cgcAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGt                <  4:85735/149‑1 (MQ=255)
                                                                                                                             cgcAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTGCAGGCCCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGt                <  1:54804/149‑1 (MQ=255)
                                                                                                                              gcAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGGTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTGCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGt                <  1:271763/149‑1 (MQ=255)
                                                                                                                                 ggACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGTTGGTGCCCGGTTCGCCGCGACCGGTGAGc            <  2:42294/149‑1 (MQ=255)
                                                                                                                                 ggACCACGTCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGc            <  6:64817/149‑1 (MQ=255)
                                                                                                                                   aCCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAg          <  7:207625/149‑1 (MQ=255)
                                                                                                                                        gTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGc     >  7:241592/1‑149 (MQ=255)
                                                                                                                                         tCCAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGcc    <  7:228791/149‑1 (MQ=255)
                                                                                                                                           cAATTTCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCAt  >  4:61589/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AGCCGCTTTGACCCTGGCCGGTAAGGTGAATGCCAGCGGTATCGATGGCCAGCCCGGAGGGTTGCAGCGGGCCGTCGGTATCGGGCGCGACGATGTTGGCGGTCGAGGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGT‑CCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCAT  >  NC_002947/935805‑936092

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