Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,084,381 G→T 6.1% D175Y (GAC→TAC)  PP_3594 → amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,084,3810GT 6.1% 206.8 / 2.3 65D175Y (GAC→TAC) PP_3594amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base T (2/2);  total (30/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 = 7.83e-01

CTTCAAGACCCTGCCCAAGGGCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGTT  >  NC_002947/4084236‑4084526
                                                                                                                                                 |                                                                                                                                                 
cTTCAAGACCCTGCCCAAGGGCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCTACAACTGGCTGTCGCTGCTCAAGGaca                                                                                                                                                >  2:93358/1‑149 (MQ=255)
cTTCAAGACCCTGCCCAAGGGCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGaca                                                                                                                                                >  7:348537/1‑149 (MQ=255)
 ttCAAGACCCTGCCCAAGGGCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCGAGGacac                                                                                                                                               >  1:205812/1‑149 (MQ=255)
      gACCCTGCCCAAGGGCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGc                                                                                                                                          <  1:118495/149‑1 (MQ=255)
        cccTGCCCAAGGGCCAGCTGGAAGCCTGCTCGGCGCTTGTGCTTAGCAAGCGCACCTGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTg                                                                                                                                        <  5:418836/149‑1 (MQ=255)
         ccTGCCCAAGGGCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTgg                                                                                                                                       >  5:297251/1‑149 (MQ=255)
             cccAAGGGCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTc                                                                                                                                   >  1:270246/1‑149 (MQ=255)
              ccAAGGGCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTc                                                                                                                                   >  3:213988/1‑148 (MQ=255)
              ccAAGGGCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCg                                                                                                                                  >  6:439346/1‑149 (MQ=255)
              ccAAGGGCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCg                                                                                                                                  >  2:145195/1‑149 (MQ=255)
                aaGGGCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGAc                                                                                                                                <  6:130098/149‑1 (MQ=255)
                    gCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCAtc                                                                                                                            <  1:145195/149‑1 (MQ=255)
                    gCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCAtc                                                                                                                            <  8:355889/149‑1 (MQ=255)
                        gCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGc                                                                                                                        >  8:6582/1‑149 (MQ=255)
                         cTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCGAGGACACCTCGCTGGTGTCGACCATCTCGCt                                                                                                                       <  2:205812/149‑1 (MQ=255)
                          tGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGCCCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTg                                                                                                                      >  8:340457/1‑149 (MQ=255)
                                ccTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGAc                                                                                                                >  5:193067/1‑149 (MQ=255)
                                ccTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCAGCGAGATGGTCGAc                                                                                                                <  8:348537/149‑1 (MQ=255)
                                       ggCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGGTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGAt                                                                                                           <  7:406525/147‑1 (MQ=255)
                                            ttGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATgcgc                                                                                                       >  1:271466/1‑146 (MQ=255)
                                              gggCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAg                                                                                                    <  6:193067/147‑1 (MQ=255)
                                              gggCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGa                                                                                                   <  2:270246/148‑1 (MQ=255)
                                                gCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGAcc                                                                                                 >  3:367916/1‑148 (MQ=255)
                                                       aaGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGc                                                                                         >  7:17259/1‑149 (MQ=255)
                                                        aGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCg                                                                                        <  7:6582/149‑1 (MQ=255)
                                                           gcACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGAc                                                                                                  <  7:62895/136‑1 (MQ=255)
                                                           gcACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGAc                                                                                                  >  8:62895/1‑136 (MQ=255)
                                                            cACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCa                                                                                    <  4:171262/149‑1 (MQ=255)
                                                               cGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGTACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCg                                                                                 >  7:336500/1‑149 (MQ=255)
                                                               cGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGTACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCg                                                                                 >  7:330601/1‑149 (MQ=255)
                                                                 gTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGt                                                                                >  6:384436/1‑148 (MQ=255)
                                                                  tttCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTTGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTaa                                                                              <  7:230707/149‑1 (MQ=255)
                                                                      ttcAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAAcc                                                                            >  3:275924/1‑147 (MQ=255)
                                                                      ttcAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAAcc                                                                            >  1:310483/1‑147 (MQ=255)
                                                                         aaGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAgg                                                                        <  8:320928/148‑1 (MQ=255)
                                                                         aaGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAgg                                                                        >  7:299762/1‑148 (MQ=255)
                                                                           ggTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAAc                                                                     <  6:297251/149‑1 (MQ=255)
                                                                                 gctgcCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGt                                                               >  2:278620/1‑149 (MQ=255)
                                                                                 gctgcCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGt                                                               >  4:356358/1‑149 (MQ=255)
                                                                                  ctgcCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGt                                                               <  1:278620/148‑1 (MQ=255)
                                                                                     ccACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGTCGGTCAGCGTAAGCAAGGAACCGATGttctt                                                           <  6:177246/149‑1 (MQ=255)
                                                                                     ccACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGttct                                                            <  4:213988/148‑1 (MQ=255)
                                                                                      cacaACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACATCTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGGTGGTCGACCTGATGCGCCATACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGttcttc                                                          <  1:93358/149‑1 (MQ=255)
                                                                                      cacaACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGttcttc                                                          <  5:384436/149‑1 (MQ=255)
                                                                                       acaACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGttcttc                                                          >  5:14875/1‑148 (MQ=255)
                                                                                        caACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTa                                                        <  5:179008/149‑1 (MQ=255)
                                                                                                 cGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACg                                                      <  8:189569/142‑1 (MQ=255)
                                                                                                 cGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACg                                                      >  7:189569/1‑142 (MQ=255)
                                                                                                  gCATCGCCCTGCCCGTCTTGTCCAACAACTGGCTGTCGCTTCTCTAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCATCCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCgcc                                              <  1:316871/149‑1 (MQ=255)
                                                                                                    aTCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCgcctg                                            <  5:91902/149‑1 (MQ=255)
                                                                                                         ccTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTgg                                       <  2:271466/149‑1 (MQ=255)
                                                                                                         ccTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTgg                                       <  3:356358/149‑1 (MQ=255)
                                                                                                              ccGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTg                                                                                                                                        <  6:335277/47‑1 (MQ=255)
                                                                                                              ccGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTg                                                                                                                                        >  5:335277/1‑47 (MQ=255)
                                                                                                              ccGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTAc                                  <  8:17259/149‑1 (MQ=255)
                                                                                                              ccGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTAc                                  <  5:439346/149‑1 (MQ=255)
                                                                                                               cGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCCAcc                                 >  7:57321/1‑149 (MQ=255)
                                                                                                                     gTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGt                           >  1:181046/1‑149 (MQ=255)
                                                                                                                           caacTGGCTGTCGCTGCTCAAGTACACCTCGCTGGTGTCGACCATCTCGCTTGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCgg                     <  8:330601/149‑1 (MQ=255)
                                                                                                                           caacTGGCTGTCGCTGCTCAAGTACACCTCGCTGGTGTCGACAATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCgg                     <  8:336500/149‑1 (MQ=255)
                                                                                                                           caacTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCgg                     <  7:340457/149‑1 (MQ=255)
                                                                                                                                     tCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCAGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTcc             <  1:174829/147‑1 (MQ=255)
                                                                                                                                           ctCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGgcgc       <  2:310483/147‑1 (MQ=255)
                                                                                                                                            tCAAGGACACCTCGCTGGTGTCGACAATCTCGCTGGTCGACATGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGAGGTTCTTCTACGGCGTAGCCTGCCTGGGCTACCTGCAGTTCGCGGCACTGTCAGGGCGCgtg    >  4:62126/1‑149 (MQ=255)
                                                                                                                                              aaGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGtt  >  1:248989/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTTCAAGACCCTGCCCAAGGGCCAGCTGGAAGCCTGCTCGGCGCTTGGGCTTAGCAAGCGCACCGGTTTCTTCAAGGTGCTGCTGCCACAACTGACCCGCATCGCCCTGCCCGGCTTGTCCAACAACTGGCTGTCGCTGCTCAAGGACACCTCGCTGGTGTCGACCATCTCGCTGGTCGACCTGATGCGCCAGACCAACCTGGCGGTCAGCGTAACCAAGGAACCGATGTTCTTCTACGGCGTCGCCTGCCTGGGCTACCTGCTGTTCGCGGCACTGTCCGGGCGCGTGTT  >  NC_002947/4084236‑4084526

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: