Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,877,513 C→T 6.4% T201I (ACC→ATC) ‡ PP_3425 → RND family transporter MFP subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,877,5130CT 6.4% 183.1 / 2.1 62T201I (ACC→ATC) ‡PP_3425RND family transporter MFP subunit
Reads supporting (aligned to +/- strand):  ref base C (27/31);  new base T (2/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

CCGGGGTCGATGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCGGTATACCTGGGCCTGACCAACGAAAGCGGCAACCCG  >  NC_002947/3877369‑3877658
                                                                                                                                                |                                                                                                                                                 
ccGGGGTCGATGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGc                                                                                                                                               <  4:77043/149‑1 (MQ=255)
 cGGGGTCGATGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCg                                                                                                                                              >  8:296170/1‑149 (MQ=255)
   gggTCGATGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTGAGCTTCACCCGCCTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCgtg                                                                                                                                            <  3:282958/149‑1 (MQ=255)
      tCGATGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCgtggtg                                                                                                                                         >  3:79941/1‑149 (MQ=255)
        gATGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCCCCGCCGAAGTCACCCCGCAGACCACCGTGGTGTc                                                                                                                                       >  4:202733/1‑149 (MQ=255)
        gATGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCACGCTGGTCAGCggggt                                                                                                                                          >  4:53191/1‑146 (MQ=255)
        gATGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCACGCTGGTCAGCggggt                                                                                                                                          <  3:53191/146‑1 (MQ=255)
         aTGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGTCCAGCGTGGTGTc                                                                                                                                       >  6:147821/1‑148 (MQ=255)
         aTGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGACGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTc                                                                                                                                       >  4:133299/1‑148 (MQ=255)
            ccATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCAcc                                                                                                                                   <  5:192917/149‑1 (MQ=255)
            ccATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCAcc                                                                                                                                   <  5:283313/149‑1 (MQ=255)
            ccATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCAcc                                                                                                                                   <  2:303950/149‑1 (MQ=255)
                        aGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTAc                                                                                                                       >  5:66699/1‑149 (MQ=255)
                                cTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGTCCAGCGTGGTGTCCACCGACAAGGTCTACGCCTAttt                                                                                                               <  5:147821/149‑1 (MQ=255)
                                cTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTAttt                                                                                                               <  3:133299/149‑1 (MQ=255)
                                     gcgcCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACg                                                                                                          >  6:159325/1‑149 (MQ=255)
                                      cgcCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGc                                                                                                         <  2:168002/149‑1 (MQ=255)
                                           gAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCcgac                                                                                                     >  5:221803/1‑148 (MQ=255)
                                               cTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCACGCTGGTCAGCggggt                                                                                                                                          <  7:146362/107‑1 (MQ=255)
                                               cTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCACGCTGGTCAGCggggt                                                                                                                                          >  8:146362/1‑107 (MQ=255)
                                                tCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAgcgc                                                                                               >  3:224907/1‑149 (MQ=255)
                                                tCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAgc                                                                                                 >  4:51299/1‑147 (MQ=255)
                                                  aGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCgt                                                                                             <  4:224907/149‑1 (MQ=255)
                                                    cTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCgtgt                                                                                           <  7:296170/149‑1 (MQ=255)
                                                    cTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCgtg                                                                                            >  8:132987/1‑148 (MQ=255)
                                                        aCCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCt                                                                                       >  3:6564/1‑149 (MQ=255)
                                                         cccGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTc                                                                                      <  7:189398/149‑1 (MQ=255)
                                                          ccGCGTCACCGCGCCCATCAGCGGCCCAGTCAGCCCCGCAAATTTCACCGCCGCCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGAGTCCACCGACAAGGTCTACGCCTATTTAGACGCAGAAGACCGCCTCTAcaac                                                                                      >  8:6656/1‑145 (MQ=255)
                                                            gcgTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAg                                                                                   <  6:66699/149‑1 (MQ=255)
                                                             cgTCAACGCGCCCAGCAGCGGCCGAGTAAGCCGCGCCGAGTTCCCCGACGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGt                                                                                  <  7:103399/149‑1 (MQ=255)
                                                                      gcCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGc                                                                         <  5:159325/149‑1 (MQ=255)
                                                                        ccATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTg                                                                       <  7:132987/149‑1 (MQ=255)
                                                                           tCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGcc                                                                    >  2:294298/1‑149 (MQ=255)
                                                                            cAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGccc                                                                   >  1:289415/1‑149 (MQ=255)
                                                                                  ccGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAg                                                             <  3:51299/149‑1 (MQ=255)
                                                                                  ccGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAACTACACCCAGCTGGCCCGTGAAg                                                             <  4:79941/149‑1 (MQ=255)
                                                                                   cGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAgg                                                            >  3:268233/1‑149 (MQ=255)
                                                                                           ccGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCa                                                         <  8:281845/144‑1 (MQ=255)
                                                                                           ccGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCa                                                         >  7:281845/1‑144 (MQ=255)
                                                                                           ccGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCg                                                      >  7:2054/1‑147 (MQ=255)
                                                                                                 cGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTTGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCaga                                              <  3:202733/149‑1 (MQ=255)
                                                                                                 cGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCaga                                              <  4:6564/149‑1 (MQ=255)
                                                                                                  gAGTTCATCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCagag                                             <  6:103101/149‑1 (MQ=255)
                                                                                                  gAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCagag                                             <  4:261981/149‑1 (MQ=255)
                                                                                                  gAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCagag                                             <  2:257548/149‑1 (MQ=255)
                                                                                                     ttCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGACTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCAc                                          >  8:233844/1‑149 (MQ=255)
                                                                                                       cACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTc                                        >  4:10747/1‑149 (MQ=255)
                                                                                                        aCCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTcc                                       >  3:165122/1‑149 (MQ=255)
                                                                                                        aCCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTcc                                       >  5:142827/1‑149 (MQ=255)
                                                                                                        aCCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTc                                        <  8:2054/148‑1 (MQ=255)
                                                                                                         ccgccgGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCg                                      <  7:6656/149‑1 (MQ=255)
                                                                                                         ccgccgGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCg                                      <  2:289415/149‑1 (MQ=255)
                                                                                                          cgccgGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCgg                                     <  1:294298/149‑1 (MQ=255)
                                                                                                          cgccgGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCg                                      >  6:209998/1‑148 (MQ=255)
                                                                                                            ccgGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCGGta                                   <  6:221803/149‑1 (MQ=255)
                                                                                                                      cGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCGGTATACCTGGGcc                         >  6:203166/1‑149 (MQ=255)
                                                                                                                        tCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCGGTATACCTGGGCCt                        >  8:188361/1‑148 (MQ=255)
                                                                                                                        tCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCGGTATACCTGGGCCTg                       <  5:203166/149‑1 (MQ=255)
                                                                                                                         cACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCGGTATACCTGGGCCTGa                      <  3:10747/149‑1 (MQ=255)
                                                                                                                           ccgccgATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCACGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCGGTATACCTGGGCCTGAc                     >  8:259500/1‑148 (MQ=255)
                                                                                                                               cgATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCGGTATACCTGGGCCTGACCAACg                <  7:188361/149‑1 (MQ=255)
                                                                                                                                             tGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCGGTATACCTGGGCCTGACCAACGAAAGCGGCAACCCg  >  1:138904/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CCGGGGTCGATGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCGGTATACCTGGGCCTGACCAACGAAAGCGGCAACCCG  >  NC_002947/3877369‑3877658

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