Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,933,122 G→T 6.4% A93A (GCC→GCA)  pepP ← proline aminopeptidase P II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,933,1220GT 6.4% 197.4 / 3.5 62A93A (GCC→GCA) pepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base T (2/2);  total (32/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 = 7.26e-01

CATCAGCCGGCGGTCGAACTCCGGGTTGCTGCCCATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCTGAGGTACTGGAAATCGCTGTCCTGGC  >  NC_002947/5932975‑5933267
                                                                                                                                                   |                                                                                                                                                 
cATCAGCCGGCGGTCGAACTCCGGGTTGCTGCCCATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTggccgg                                                                                                                                                  <  5:337105/149‑1 (MQ=255)
   cAGCCGGCGGTCGAACTCCGGGTTGCTGCCCATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTggccggcc                                                                                                                                                >  6:105209/1‑148 (MQ=255)
   cAGCCGGCGGTCGAACTCCGGGTTGCTGCCCATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTggccggcc                                                                                                                                                <  5:105209/148‑1 (MQ=255)
     gCCGGCGGTCGAACTCCGGGTTGCTGCCCATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCg                                                                                                                                              >  1:46809/1‑148 (MQ=255)
        ggcggTCGAACTCCGGGTTGCTGCCCATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAgg                                                                                                                                          <  6:203771/149‑1 (MQ=255)
           ggTCGAACTCCGGGTTGCTGCCCATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCg                                                                                                                                       >  5:257572/1‑149 (MQ=255)
             tCGAACTCCGGGTTGCTGCCCATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGCTATCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTc                                                                                                                                     >  5:265316/1‑149 (MQ=255)
               gAACTCCGGGTTGCTGCCCATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTcc                                                                                                                                    <  6:111433/148‑1 (MQ=255)
                   tCCGGGTTGCTGCCCATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAg                                                                                                                               >  6:49066/1‑149 (MQ=255)
                     cGGGTTGCTGCCCATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCt                                                                                                                             >  2:68072/1‑149 (MQ=255)
                             tccccATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCt                                                                                                                       >  3:246706/3‑147 (MQ=255)
                               cccATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCTGCGggccggcc                                                                                                                                        <  8:314012/128‑1 (MQ=255)
                               cccATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCTGCGggccggcc                                                                                                                                        >  7:314012/1‑128 (MQ=255)
                               cccATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAg                                                                                                                    <  5:49066/148‑1 (MQ=255)
                                 cATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCTGCGggccggccagga                                                                                                                                    <  8:198398/130‑5 (MQ=255)
                                 cATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCTGCGggccggccagga                                                                                                                                    >  7:198398/1‑126 (MQ=255)
                                 cATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTcc                                                                                                                                    <  6:266226/130‑1 (MQ=255)
                                 cATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTcc                                                                                                                                    >  5:266226/1‑130 (MQ=255)
                                 cATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGt                                                                                                                 <  7:130622/149‑1 (MQ=255)
                                     gcgcTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGTCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCg                                                                                                             <  4:246706/149‑1 (MQ=255)
                                      cgcTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGt                                                                                                            >  4:278772/1‑149 (MQ=255)
                                       gcTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGtt                                                                                                           <  1:68072/149‑1 (MQ=255)
                                                acGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAg                                                                                                   >  3:328595/1‑148 (MQ=255)
                                                 cGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGaa                                                                                                 <  4:259826/149‑1 (MQ=255)
                                                 cGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGaa                                                                                                 <  6:257572/149‑1 (MQ=255)
                                                      gcgGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTTATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGc                                                                                             <  3:86685/148‑1 (MQ=255)
                                                           cTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCACCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTAc                                                                                       >  2:120167/1‑149 (MQ=255)
                                                              cGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTc                                                                                     <  5:286083/148‑1 (MQ=255)
                                                               gATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCg                                                                                    <  6:38160/148‑1 (MQ=255)
                                                                     gccaggcAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCt                                                                             >  4:147363/1‑149 (MQ=255)
                                                                      ccaggcAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCGGAACAGCACGTACTCGCCATGCTc                                                                            <  2:71413/149‑1 (MQ=255)
                                                                       caggcaggATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCACTGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTcg                                                                           >  3:90006/1‑149 (MQ=255)
                                                                        aggcaggATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTcgc                                                                          <  8:100415/149‑1 (MQ=255)
                                                                          gcaggATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCgcgc                                                                        >  5:8291/1‑149 (MQ=255)
                                                                                    gtcgATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAAt                                                                <  4:328595/147‑1 (MQ=255)
                                                                                    gtcgATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCa                                                              <  3:147363/149‑1 (MQ=255)
                                                                                       gATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAAc                                                            <  8:17198/148‑1 (MQ=255)
                                                                                                    gggAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTc                                              <  2:46809/149‑1 (MQ=255)
                                                                                                    gggAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCACGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTc                                              >  7:239780/1‑149 (MQ=255)
                                                                                                      gAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTc                                                                            <  4:169092/117‑1 (MQ=255)
                                                                                                      gAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTc                                                                            >  3:169092/1‑117 (MQ=255)
                                                                                                       aaGGCATCGTCGGCGCCGAAATCACGCAACGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTAGCGCCCTGGAATCAACGCAATCACCGCCTCggg                                           >  4:26930/1‑149 (MQ=255)
                                                                                                           cATCCTCGGCTCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCACGCCGTCCCACAGCTCCCGCTCAGGGATGCGTACACGGCAGAACAGCACGTACTCGCCAAGAACGCGCCCTCGAAACAACGCAATCACCCCCTCGGCCTcc                                       >  2:28172/1‑149 (MQ=255)
                                                                                                            aTCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCg                                      >  5:115556/1‑149 (MQ=255)
                                                                                                             tcgtcgGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCgg                                     <  3:278772/149‑1 (MQ=255)
                                                                                                                   gcgcCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACg                                                           >  2:350495/1‑121 (MQ=255)
                                                                                                                   gcgcCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACg                                                           <  1:350495/121‑1 (MQ=255)
                                                                                                                     gcCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGCGCTCCGGGAAcccc                              >  8:307030/1‑147 (MQ=255)
                                                                                                                       cGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACgg                                                                                                      <  6:166655/74‑1 (MQ=255)
                                                                                                                       cGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACgg                                                                                                      >  5:93670/1‑74 (MQ=255)
                                                                                                                       cGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACgg                                                                                                      >  5:166655/1‑74 (MQ=255)
                                                                                                                       cGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACgg                                                                                                      <  6:93670/74‑1 (MQ=255)
                                                                                                                       cGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCt                            >  5:42576/1‑148 (MQ=255)
                                                                                                                       cGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCt                            >  3:204592/1‑148 (MQ=255)
                                                                                                                         aaaTCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCTGAg                         <  7:307030/149‑1 (MQ=255)
                                                                                                                         aaaTCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCTGAg                         >  1:42158/1‑149 (MQ=255)
                                                                                                                            tCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTTAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCTGAGGTa                      <  6:42576/149‑1 (MQ=255)
                                                                                                                             cACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCTGAGGTAc                     >  7:200373/1‑149 (MQ=255)
                                                                                                                                gCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGGCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCTGAGGTACTgg                  <  1:28172/149‑1 (MQ=255)
                                                                                                                                  aCCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCTGAGGTACTGGa                 >  5:153644/1‑148 (MQ=255)
                                                                                                                                      ccccTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCAACGCCTCGGGCTCCGGGAAGCCCCTGAGGTACTGGAAATCg            >  2:249832/1‑149 (MQ=255)
                                                                                                                                        ccTTGCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCTGAGGTACTGGAAATCGCt          <  6:265316/149‑1 (MQ=255)
                                                                                                                                         cTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGc                                                                                             <  7:232843/65‑1 (MQ=255)
                                                                                                                                         cTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGc                                                                                             >  8:232843/1‑65 (MQ=255)
                                                                                                                                              tGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCTGAGGTACTGGAAATCGCTGTCCTg    >  3:166589/1‑149 (MQ=255)
                                                                                                                                               ggccggccCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCTGAGGTACTGGAAATCGCTGTCCTgg   >  6:73364/1‑149 (MQ=255)
                                                                                                                                                gccggccCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTg                                                                    >  2:48415/1‑83 (MQ=255)
                                                                                                                                                gccggccCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTg                                                                    <  1:48415/83‑1 (MQ=255)
                                                                                                                                                gccggccCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCTGAGGTACTGGAAATCGCTGTCCTGGc  >  8:172453/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CATCAGCCGGCGGTCGAACTCCGGGTTGCTGCCCATGGCGCTGTAGACACGCTCGCGGCCTTCGATCAGGCCAGGCAGGATCTCGTCGATATCCGTGATGGGGAAGGCATCGTCGGCGCCGAAATCACGCACCGCCCCTTCCTGGCCGGCCCGCAGGCCGTCCCACAGCTCGCGCTCAGGGTTGCGTTCACGGCAGAACAGCACGTACTCGCCATGCTCGCGCCCTGGAATCAACGCAATCACCGCCTCGGGCTCCGGGAAGCCGCTGAGGTACTGGAAATCGCTGTCCTGGC  >  NC_002947/5932975‑5933267

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