Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,072,904 G→A 5.9% P451L (CCG→CTG)  PP_5327 ← phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,072,9040GA 5.9% 258.4 / 3.9 84P451L (CCG→CTG) PP_5327phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (48/31);  new base A (3/2);  total (51/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 = 3.61e-01

CATGATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAGGG  >  NC_002947/6072759‑6073051
                                                                                                                                                 |                                                                                                                                                   
cATGATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGc                                                                                                                                                  >  3:164497/1‑149 (MQ=255)
  tGATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTg                                                                                                                                                >  6:291822/1‑149 (MQ=255)
  tGATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTg                                                                                                                                                >  2:170469/1‑149 (MQ=255)
   gATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATGAGCGGCGCCACATCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTgg                                                                                                                                               >  5:156150/1‑149 (MQ=255)
            gTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATggg                                                                                                                                        <  7:213625/147‑1 (MQ=255)
            gTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCAGCCCGGCCATga                                                                                                                                         <  7:335394/146‑2 (MQ=255)
                    gCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGAtt                                                                                                                              >  5:187942/1‑149 (MQ=255)
                      aGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGAtttt                                                                                                                            <  1:323803/149‑1 (MQ=255)
                       ggTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATg                                                                                                                                          <  6:75987/134‑1 (MQ=255)
                       ggTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATg                                                                                                                                          >  5:75987/1‑134 (MQ=255)
                       ggTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTc                                                                                                                           >  2:27794/1‑149 (MQ=255)
                        gTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTc                                                                                                                           >  1:36228/1‑148 (MQ=255)
                        gTACGCGTAGTTTCCATCCACCGCCAGCGACGGCGCCAGCTTCAACACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCAGGCCGGTCATCATGCCCGGGCTGGCCATGGGCAGGACGATTTTc                                                                                                                           >  8:268741/1‑148 (MQ=255)
                         tACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCACCATGGCCAGCCCGGCCATgat                                                                                                                                        >  7:285118/1‑132 (MQ=255)
                         tACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCACCATGGCCAGCCCGGCCATgat                                                                                                                                        <  8:285118/134‑3 (MQ=255)
                             ggTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGGGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAgg                                                                                                                     <  7:353493/149‑1 (MQ=255)
                               tAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGt                                                                                                                   >  4:283471/1‑149 (MQ=255)
                                  tttCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGtctc                                                                                                                >  2:59760/1‑149 (MQ=255)
                                   ttCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCg                                                                                                               >  6:382702/1‑149 (MQ=255)
                                     ccatccaCCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGc                                                                                                             >  6:204799/1‑149 (MQ=255)
                                     ccatccaCCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGc                                                                                                             >  7:329840/1‑149 (MQ=255)
                                       atccaCCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCt                                                                                                           <  1:27794/149‑1 (MQ=255)
                                          caCCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTgg                                                                                                        <  2:371060/149‑1 (MQ=255)
                                           aCCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCt                                                                                                           >  6:219593/1‑145 (MQ=255)
                                           aCCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCt                                                                                                           <  5:219593/145‑1 (MQ=255)
                                           aCCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGCTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGt                                                                                                       >  7:219556/1‑149 (MQ=255)
                                             cgccAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAgg                                                                                                                                    <  5:145111/115‑1 (MQ=255)
                                             cgccAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAgg                                                                                                                                    >  6:145111/4‑118 (MQ=255)
                                               gCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGACCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTgcg                                                                                                   >  7:99766/1‑149 (MQ=255)
                                                cAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCCGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTgcgc                                                                                                  <  3:119522/149‑1 (MQ=255)
                                                cAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTgcgc                                                                                                  <  2:36228/149‑1 (MQ=255)
                                                cAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTgcgc                                                                                                  <  6:187942/149‑1 (MQ=255)
                                                 aGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGcc                                                                                                 <  5:234570/149‑1 (MQ=255)
                                                  gCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGcc                                                                                                 >  1:300182/1‑148 (MQ=255)
                                                    gACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCaa                                                                                              <  2:6009/149‑1 (MQ=255)
                                                    gACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCaa                                                                                              >  6:27466/1‑149 (MQ=255)
                                                    gACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCaa                                                                                              >  6:32955/1‑149 (MQ=255)
                                                     aCGGCGCCCGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCCCCGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGCCTATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAAc                                                                                             >  5:169127/1‑149 (MQ=255)
                                                            cAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCg                                                                                      >  5:176600/1‑149 (MQ=255)
                                                             aGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGACCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGa                                                                                     <  8:99766/149‑1 (MQ=255)
                                                              gCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGc                                                                                       >  4:8402/1‑146 (MQ=255)
                                                              gCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGa                                                                                     <  3:8402/148‑1 (MQ=255)
                                                                ttCACCACACCTGCCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGcc                                                                                  >  2:352287/1‑149 (MQ=255)
                                                                ttCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGcc                                                                                  >  7:11400/1‑149 (MQ=255)
                                                                   accacACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCtt                                                                                >  6:81355/1‑148 (MQ=255)
                                                                    ccacACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTca                                                                              >  6:97333/1‑149 (MQ=255)
                                                                     cacaCCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTcac                                                                             >  6:222098/1‑149 (MQ=255)
                                                                     cacaCCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTcac                                                                             <  2:300182/149‑1 (MQ=255)
                                                                          cTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTcacgcacg                                                                        <  5:27466/149‑1 (MQ=255)
                                                                          cTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTcacgcacg                                                                        <  5:32955/149‑1 (MQ=255)
                                                                           tACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACgg                                                                       >  7:272645/1‑149 (MQ=255)
                                                                               aGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTgcg                                                                   >  5:110646/1‑149 (MQ=255)
                                                                                  aTCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGcgcg                                                                 >  7:258672/1‑148 (MQ=255)
                                                                                      gcggcgCCACTTCGCCGGCGGCACGGGCCACAGCGATGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAAtt                                                            <  4:196425/149‑1 (MQ=255)
                                                                                      gcggcgCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTca                                                                              <  3:276426/131‑1 (MQ=255)
                                                                                      gcggcgCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTca                                                                              >  4:276426/1‑131 (MQ=255)
                                                                                               cTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAgg                                                   <  1:352287/149‑1 (MQ=255)
                                                                                                  cGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGcc                                                 <  2:284613/148‑1 (MQ=255)
                                                                                                      ggcggcACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCcttctt                                            <  6:118002/149‑1 (MQ=255)
                                                                                                       gcggcACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCttcttc                                           <  5:291822/149‑1 (MQ=255)
                                                                                                         ggcACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCgg                                         >  7:73735/1‑149 (MQ=255)
                                                                                                            aCGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTgg                                      >  4:83284/1‑149 (MQ=255)
                                                                                                                gCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCAc                                  >  4:208167/1‑149 (MQ=255)
                                                                                                                    caGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATc                              <  8:73735/149‑1 (MQ=255)
                                                                                                                      gCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATcac                            <  6:110646/149‑1 (MQ=255)
                                                                                                                       cGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATcac                            >  3:97189/1‑148 (MQ=255)
                                                                                                                         aGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATcaccac                         <  5:382702/149‑1 (MQ=255)
                                                                                                                            aTCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCgg                      >  4:363575/1‑149 (MQ=255)
                                                                                                                              cATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTa                    >  5:134839/1‑149 (MQ=255)
                                                                                                                               aTGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCACCGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTAc                   >  4:74700/1‑149 (MQ=255)
                                                                                                                                tGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTAcc                  >  5:131242/1‑149 (MQ=255)
                                                                                                                                ggCCGGTCATCATCGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGCGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTAcc                  <  8:258672/148‑1 (MQ=255)
                                                                                                                                   tcctgcccatGGCCAGCCCGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGc                >  2:327440/8‑148 (MQ=255)
                                                                                                                                   cGGTCATCATGGCCGGGCTGGCCACGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGc                >  7:115076/1‑148 (MQ=255)
                                                                                                                                    ggTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACACCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGc                >  6:246861/1‑147 (MQ=255)
                                                                                                                                      tcctcatGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAg             >  2:305148/4‑148 (MQ=255)
                                                                                                                                         tcatGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAAt         <  1:170469/149‑1 (MQ=255)
                                                                                                                                         tcatGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAAt         >  3:185228/1‑149 (MQ=255)
                                                                                                                                         ccatGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCACCGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAAt         >  2:92551/2‑149 (MQ=255)
                                                                                                                                         ccatGGCCAGCCCGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAAt         >  6:314561/2‑149 (MQ=255)
                                                                                                                                          catGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATg        >  1:340656/1‑149 (MQ=255)
                                                                                                                                           atGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGTGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATg        <  6:374149/148‑1 (MQ=255)
                                                                                                                                               ccGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAgg   <  4:97189/149‑1 (MQ=255)
                                                                                                                                                cGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAggg  >  8:356130/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CATGATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAGGG  >  NC_002947/6072759‑6073051

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: