Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,539,680 G→A 6.6% P297L (CCC→CTC)  PP_1350 ← two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,539,6800GA 6.6% 178.3 / 2.5 60P297L (CCC→CTC) PP_1350two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base A (2/2);  total (29/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 = 9.79e-01

GGAGGAACGGCTTCAACGGCTCGAACGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGCCCTGAGCCTGGTCGTAGGCA  >  NC_002947/1539540‑1539823
                                                                                                                                            |                                                                                                                                                
ggaggaACGGCTTCAACGGCTCGAACGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAg                                                                                                                                          >  4:340077/1‑149 (MQ=255)
     aaCGGCTTCAACGGCTCGAACGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCgat                                                                                                                                     <  8:102364/149‑1 (MQ=255)
      aCGGCTTCAACGGCTCGAACGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAgg                                                                                                                                                  <  1:63307/135‑1 (MQ=255)
      aCGGCTTCAACGGCTCGAACGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAgg                                                                                                                                                  >  2:63307/1‑135 (MQ=255)
       cGGCTTCAACGGCTCGAACGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCgatga                                                                                                                                   <  3:290915/149‑1 (MQ=255)
          cTTCAACGGCTCGAACGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGc                                                                                                                                >  3:295608/1‑149 (MQ=255)
              aaCGGCTCGAACGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGtt                                                                                                                            <  6:340460/149‑1 (MQ=255)
                 ggCTCGAACGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGAGCATCGAGTCGATGATGCTGtt                                                                                                                            >  3:64365/1‑146 (MQ=255)
                 ggCTCGAACGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGAGCATCGAGTCGATGATGCTGtt                                                                                                                            >  3:59518/1‑146 (MQ=255)
                      gAACGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATa                                                                                                                      <  5:211/147‑1 (MQ=255)
                       aaCGCAATGAACACCGTCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGtt                                                                                                                   <  5:198538/149‑1 (MQ=255)
                        aCGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTc                                                                                                                  >  1:14517/1‑149 (MQ=255)
                         cGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCt                                                                                                                 <  7:172329/149‑1 (MQ=255)
                          gCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTg                                                                                                                >  3:203393/1‑149 (MQ=255)
                                 acacCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCg                                                                                                         >  5:238170/1‑149 (MQ=255)
                                  cacCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCgg                                                                                                        <  2:14517/149‑1 (MQ=255)
                                   acCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGt                                                                                                       <  5:133519/149‑1 (MQ=255)
                                     cGGCTGGTTAAGTGGCTCGTCCAGTGGCGTGCCGGAAAGCTCCGTGGCTTGGTTGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGa                                                                                                     <  8:79218/149‑1 (MQ=255)
                                     cGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGt                                                                                                       >  7:3050/1‑147 (MQ=255)
                                     cGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGa                                                                                                     <  5:184710/149‑1 (MQ=255)
                                       gCTGGTTAAGTGCCTCGTCGAGTGGCGTGCCGGAAAGCACCGTGTCTTCGTGGTTGCATTGGGTCACGTAAAGCTGTTCGTCGTGGGCGATCAGTGCCTAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGc                                                                                                   <  2:74915/149‑1 (MQ=255)
                                                 tGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGa                                                                                        >  2:13201/1‑149 (MQ=255)
                                                   ccTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATc                                                                                      >  7:168262/1‑149 (MQ=255)
                                                   ccTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATc                                                                                      <  4:295608/149‑1 (MQ=255)
                                                    cTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCt                                                                                     >  2:301681/1‑149 (MQ=255)
                                                               gCGTGCCGGAAAGCACCGTGGCTTCGTTGTTCCATT‑GGTCACGTAATGCTGTTCGTCGAGGTCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCt                                                                         <  5:161188/149‑1 (MQ=255)
                                                                cGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCt                                                                         <  3:340077/149‑1 (MQ=255)
                                                                 gTGCCGGAAAGCACCGTGGCTTCGTGGTTCGATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTg                                                                        >  5:313127/1‑149 (MQ=255)
                                                                         aaGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACttcca                                                                <  4:203393/149‑1 (MQ=255)
                                                                                gTGGCTTAGTGGTTCCATTGGGACACCTAAAGCTGAACGTCGAGGGCGATCAGACCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCATTTT‑CTCGATCATGCTGCGTACCTGTACttccagttccag                                                         >  4:284734/1‑149 (MQ=255)
                                                                                 tGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGc                                                        >  3:226246/1‑149 (MQ=255)
                                                                                     ttCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGc                                                    >  7:13142/1‑149 (MQ=255)
                                                                                       cGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGc                                                    <  1:13201/147‑1 (MQ=255)
                                                                                       cGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCg                                                   <  8:3050/148‑1 (MQ=255)
                                                                                        gTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGt                                                 >  1:361189/1‑149 (MQ=255)
                                                                                           gTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTgg                                              >  4:150577/1‑149 (MQ=255)
                                                                                               cATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGgtgt                                           <  8:168262/148‑1 (MQ=255)
                                                                                                aTTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTgg                                         >  7:292128/1‑149 (MQ=255)
                                                                                                aTTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTgg                                         >  8:329732/1‑149 (MQ=255)
                                                                                                 ttGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGc                                        <  6:238170/149‑1 (MQ=255)
                                                                                                           tAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTttcttc                              <  2:199270/149‑1 (MQ=255)
                                                                                                             aaGGTGTTCGTCGATGGCGATCAGTGCCGAGGGCATCGAGTGGATGATGCTGTTGAGATAGTTCTGGATACCGGTGAGCTTTTT‑CTCGATCTTGCTTAGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCg                             <  3:284734/148‑1 (MQ=255)
                                                                                                               gCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGc                           >  8:48089/1‑148 (MQ=255)
                                                                                                               gCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGc                           >  5:371601/1‑148 (MQ=255)
                                                                                                                  gTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCaa                         >  4:263965/1‑147 (MQ=255)
                                                                                                                  gTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCaa                         <  3:263965/147‑1 (MQ=255)
                                                                                                                     cgtcgAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGcc                     <  4:226246/148‑1 (MQ=255)
                                                                                                                      gtcgAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGCCCt                   <  3:150577/149‑1 (MQ=255)
                                                                                                                       tcgAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTg                                               >  6:90933/1‑120 (MQ=255)
                                                                                                                       tcgAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTg                                               <  5:90933/120‑1 (MQ=255)
                                                                                                                                aTCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGCCCTGAGCCTGGTc         <  7:329732/149‑1 (MQ=255)
                                                                                                                                 tCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTTACTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGCCCTGAGCCTGGTc         >  1:289144/1‑149 (MQ=255)
                                                                                                                                   aGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCgg                                      <  7:311287/117‑1 (MQ=255)
                                                                                                                                   aGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCgg                                      >  8:311287/1‑117 (MQ=255)
                                                                                                                                   aGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGCCCTGAGCCTGGTCGTa      <  8:13142/149‑1 (MQ=255)
                                                                                                                                    gTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGCCCTGAGCCTGGTCGTAg     <  6:371601/149‑1 (MQ=255)
                                                                                                                                    gTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGCCCTGAGCCTGGTCGTAg     >  4:144313/1‑149 (MQ=255)
                                                                                                                                    gTGCCGAGAGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGCCCTGAGCCTGGTCGTAg     <  4:64365/149‑1 (MQ=255)
                                                                                                                                    gTGCCGAGAGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGCCCTGAGCCTGGTCGTAg     <  4:59518/149‑1 (MQ=255)
                                                                                                                                       ccGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGTGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGCCCTGAGCCTGGTCGTAGGCa  >  3:306424/1‑149 (MQ=255)
                                                                                                                                       ccGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGCCCTGAGCCTGGTCGTAGGc   >  3:77934/1‑148 (MQ=255)
                                                                                                                                       ccGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGCCCTGAGCCTGGTCGTAGGCa  >  7:31353/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                
GGAGGAACGGCTTCAACGGCTCGAACGCAATGAACACCGGCTGGTTAAGTGCCTCGTCCAGTGGCGTGCCGGAAAGCACCGTGGCTTCGTGGTTCCATTGGGTCACGTAAAGCTGTTCGTCGAGGGCGATCAGTGCCGAGGGCATCGAGTCGATGATGCTGTTGAGATAGTTCTGGAAACCGGTGAGCTTTTT‑CTCGATCTTGCTGCGTACCTGTACTTCCAGTTCCAGCTTGCGGTTGGTGTGGCGGGTTTCTTCGGCAAGGCCCTGAGCCTGGTCGTAGGCA  >  NC_002947/1539540‑1539823

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