Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,687,705 C→G 6.7% V227V (GTC→GTG)  PP_4148 → ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,687,7050CG 6.7% 190.5 / 2.8 60V227V (GTC→GTG) PP_4148ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base G (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 = 8.86e-01

TGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACGTGTTCCGCAACGCCATGCTGCTGG  >  NC_002947/4687569‑4687853
                                                                                                                                         |                                                                                                                                                     
tGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTc                                                                                                                                           >  3:256800/1‑149 (MQ=255)
tGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTc                                                                                                                                           <  5:168317/149‑1 (MQ=255)
 gCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCg                                                                                                                                          >  8:7321/1‑149 (MQ=255)
  cTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGc                                                                                                                                         <  1:299883/149‑1 (MQ=255)
    gTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGc                                                                                                                                                 >  1:276147/1‑139 (MQ=39)
    gTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGc                                                                                                                                                 <  2:276147/139‑1 (MQ=39)
       cGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGccacc                                                                                                                                     <  5:18758/148‑1 (MQ=255)
         ccGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAGCAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGCa                                                                                                                                                <  2:292963/135‑1 (MQ=255)
         ccGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAGCAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGCa                                                                                                                                                >  1:292963/1‑135 (MQ=255)
                      cTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGAcc                                                                                                                      <  1:6223/148‑1 (MQ=255)
                            cTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACa                                                                                                               <  5:54927/149‑1 (MQ=255)
                             tGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTc                                                                                                                                                       <  8:308861/108‑1 (MQ=255)
                             tGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTc                                                                                                                                                       >  7:308861/1‑108 (MQ=255)
                                  cccGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTc                                                                                                          <  5:169877/148‑1 (MQ=255)
                                  cccGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTc                                                                                                          <  6:371925/148‑1 (MQ=255)
                                    cGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTc                                                                                                       >  2:156141/1‑149 (MQ=255)
                                     gCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGc                                                                                                                                         >  1:48539/1‑114 (MQ=255)
                                     gCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGc                                                                                                                                         <  2:48539/114‑1 (MQ=255)
                                     gCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTc                                                                                                       >  1:100744/1‑148 (MQ=255)
                                               atCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGACGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGAACCTGCTGACAAAGAACAGCTTCCTCGACGAAATcaa                                                                                            >  1:242205/2‑149 (MQ=255)
                                                tCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATcaa                                                                                            >  3:177477/1‑148 (MQ=255)
                                                   ccTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCaacaaa                                                                                        <  2:100744/149‑1 (MQ=255)
                                                    cTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCaacaaa                                                                                        >  7:108345/1‑148 (MQ=255)
                                                      ccAACAACTTCGACGAACTCAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAg                                                                                     <  4:307969/149‑1 (MQ=255)
                                                       caacaaCTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTAATCGGCAACTTATCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGt                                                                                    <  5:228996/149‑1 (MQ=255)
                                                        aacaacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTa                                                                                   >  6:37469/1‑149 (MQ=255)
                                                        aacaacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTa                                                                                   <  8:108345/149‑1 (MQ=255)
                                                          caacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACg                                                                                 >  3:345372/1‑149 (MQ=255)
                                                           aacTTCGACGCACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTAGACACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAAGCAACAACCAGTACGt                                                                                >  6:279898/1‑149 (MQ=39)
                                                           aacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCAGCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGt                                                                                <  5:228487/149‑1 (MQ=255)
                                                           aacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGt                                                                                >  8:301187/1‑149 (MQ=255)
                                                            acTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTg                                                                               <  3:155424/149‑1 (MQ=255)
                                                              ttCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGt                                                                             >  1:146550/1‑149 (MQ=255)
                                                              ttCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGt                                                                             >  6:193/1‑149 (MQ=255)
                                                                 gacgaACTTAGCACCACCTGCATGGTGCTGGACTACTTCTGGCACGTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAATAACAGCTTCGTCGACGAAATCAACAAACAGTACGTGGTCAc                                                                          <  5:279898/149‑1 (MQ=255)
                                                                     aaCTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTTACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTTGTCACCgcc                                                                      <  4:177477/149‑1 (MQ=255)
                                                                     aaCTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgc                                                                       >  4:371648/1‑148 (MQ=255)
                                                                          aGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaaggc                                                                 >  3:367183/1‑149 (MQ=255)
                                                                             accaccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaaggcc                                                                >  8:205464/1‑147 (MQ=255)
                                                                               caccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGccaaggccaagg                                                            >  5:203673/1‑149 (MQ=255)
                                                                               caccGGCAAGGTGCTGGACTACTTCTGGCACATGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGATATCAACAAACAGTACGTGGTCACCCccaaggccaagg                                                            >  6:283744/1‑149 (MQ=255)
                                                                                   ggCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCt                                                        <  1:66033/149‑1 (MQ=255)
                                                                                    gCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTTGTCACCGCCAAGGCCAAGGGCCTg                                                       <  6:203673/149‑1 (MQ=255)
                                                                                       aGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGc                                                    <  4:367183/149‑1 (MQ=255)
                                                                                        ggTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGcc                                                   <  7:163670/149‑1 (MQ=255)
                                                                                         gTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGcc                                                   >  8:86352/1‑148 (MQ=255)
                                                                                             tGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCa                                              >  3:291719/1‑149 (MQ=255)
                                                                                               gACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACg                                            >  7:286037/1‑149 (MQ=255)
                                                                                                actactTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACgg                                           <  2:23528/149‑1 (MQ=255)
                                                                                                actactTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACgg                                           <  5:193/149‑1 (MQ=255)
                                                                                                     ttCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACAAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCt                                      >  7:317586/1‑149 (MQ=255)
                                                                                                         ggCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTAc                                 >  5:181231/1‑149 (MQ=255)
                                                                                                          gCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTAc                                 <  3:371648/148‑1 (MQ=255)
                                                                                                                 gTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCAc                           <  6:181231/147‑1 (MQ=255)
                                                                                                                 gTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCAACATGACCCTGCTGACCTAGAACAGCTTCGTCGACGAAATCAACATTCAGTACGGTGTCACCGCCAAGGCCAAGGTCCTTAGCCTATCACGGGTTCT‑GTACTGCCACgt                         <  2:242205/149‑1 (MQ=255)
                                                                                                                        cGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACgt                         >  3:126207/1‑142 (MQ=255)
                                                                                                                        cGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACgt                         <  4:126207/142‑1 (MQ=255)
                                                                                                                              cc‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGATGCT‑GTACGGCCACGTGTTCCGCAACGcc            <  1:156141/149‑1 (MQ=255)
                                                                                                                               c‑gcgCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACGTGTTCCGCAACGCCa           >  3:232101/1‑149 (MQ=255)
                                                                                                                                   gcTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACGTGTTCCGCAACGCCAtgc        <  4:256800/149‑1 (MQ=255)
                                                                                                                                     tGGTCATCGGAAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACGTGTTCCGCAACGCCa           >  2:3213/1‑144 (MQ=255)
                                                                                                                                     tGGTCATCGGAAACATCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAACGCCAAGGGCCTGAGCCGATCACGGGTGCTCGTACGGCCACGTGTTCCGCAACGCCa           <  1:3213/145‑1 (MQ=255)
                                                                                                                                         cATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCTATCACGGGTGCT‑GTACGGCCACGTGTTCCGCAACGCCATgctgcttg  >  7:141049/1‑147 (MQ=255)
                                                                                                                                         cATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACGTGTTCCGCAACGCCATGCTGCTgg  >  7:146029/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                     
TGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCT‑GTACGGCCACGTGTTCCGCAACGCCATGCTGCTGG  >  NC_002947/4687569‑4687853

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑