Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,877,506 C→A 6.1% P199T (CCG→ACG)  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,5060CA 6.1% 205.4 / 3.0 66P199T (CCG→ACG) PP_3425RND family transporter MFP subunit
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base A (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01

AAGCCCGTGCCGGGGTCGATGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCGGTATACCTGGGCCTGACCAACG  >  NC_002947/3877360‑3877644
                                                                                                                                                  |                                                                                                                                          
aaGCCCGTGCCGGGGTCGATGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCCAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCg                                                                                                                                          >  8:189398/1‑149 (MQ=255)
         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)
                             gCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCg                                                                                                                                          >  6:135270/1‑120 (MQ=255)
                             gCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCg                                                                                                                                          <  5:135270/120‑1 (MQ=255)
                             gCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCg                                                                                                                                          >  5:117861/1‑120 (MQ=255)
                             gCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCg                                                                                                                                          <  6:117861/120‑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                         >  5:142827/1‑149 (MQ=255)
                                                                                                                 aCCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTcc                         >  3:165122/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)
                                                                                                                                                  |                                                                                                                                          
AAGCCCGTGCCGGGGTCGATGCCATCCAGGCCCAGCTCGACCTGGCGCGCCTGAACCTCAGCTTCACCCGCGTCACCGCGCCCATCAGCGGCCGAGTCAGCCGCGCCGAGTTCACCGCCGGCAACATCGTCACCGCCGATGTCACCCCGCTGACCAGCGTGGTGTCCACCGACAAGGTCTACGCCTATTTCGACGCCGACGAGCGCGTGTACCTCAAGTACACCCAGCTGGCCCGTGAAGGCCAGCGTAGCCAGAGCACTCCGGTATACCTGGGCCTGACCAACG  >  NC_002947/3877360‑3877644

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