Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,072,906 G→C 6.0% S450R (AGC→AGG) ‡ 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,9060GC 6.0% 269.0 / 6.2 84S450R (AGC→AGG) ‡PP_5327phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (48/31);  new base C (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 = 7.44e-01

CATGATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAGGGTCA  >  NC_002947/6072759‑6073054
                                                                                                                                                   |                                                                                                                                                    
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                                                                                                 >  6:27466/1‑149 (MQ=255)
                                                    gACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCaa                                                                                                 <  2:6009/149‑1 (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)
                                                                                                                                                   gCTGGCCATGGGCAGGACGATTTTCCCCCGGGGCGCGGCCTGGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAGGGTCa  <  6:169127/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CATGATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAGGGTCA  >  NC_002947/6072759‑6073054

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: