Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,072,908 T→C 5.9% S450G (AGC→GGC) ‡ 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,9080TC 5.9% 253.5 / 4.2 84S450G (AGC→GGC) ‡PP_5327phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (47/32);  new base C (3/2);  total (50/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

TGATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAGGGTCAGC  >  NC_002947/6072761‑6073056
                                                                                                                                                   |                                                                                                                                                    
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                                                                                                                                               >  5:75987/1‑134 (MQ=255)
                     ggTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATg                                                                                                                                               <  6:75987/134‑1 (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                                                                                                                  >  7:329840/1‑149 (MQ=255)
                                   ccatccaCCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGc                                                                                                                  >  6:204799/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                                                                                                                                         >  6:145111/4‑118 (MQ=255)
                                           cgccAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAgg                                                                                                                                         <  5:145111/115‑1 (MQ=255)
                                             gCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGACCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTgcg                                                                                                        >  7:99766/1‑149 (MQ=255)
                                              cAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCCGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTgcgc                                                                                                       <  3:119522/149‑1 (MQ=255)
                                              cAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTgcgc                                                                                                       <  6:187942/149‑1 (MQ=255)
                                              cAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTgcgc                                                                                                       <  2:36228/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:32955/1‑149 (MQ=255)
                                                  gACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCaa                                                                                                   >  6:27466/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                                                                                  <  2:300182/149‑1 (MQ=255)
                                                                   cacaCCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTcac                                                                                  >  6:222098/1‑149 (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                                                                                   >  4:276426/1‑131 (MQ=255)
                                                                                    gcggcgCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTca                                                                                   <  3:276426/131‑1 (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              >  3:185228/1‑149 (MQ=255)
                                                                                                                                       tcatGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAAt              <  1:170469/149‑1 (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)
                                                                                                                                                   tGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAGGGTCAGc  >  6:376372/1‑149 (MQ=255)
                                                                                                                                                   tGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACCCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAGGGTCAGc  <  1:305148/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGATCTTCTGGTCCAGGTGCAGGTACGGGTAGTTTCCATCCACCGGCAGCGACGGCGCCAGCTTCACCACACCTACCAGCATCAGCGGCGCCACTTCGCCGGCGGCACGGGCCACAGCGAGGATCATGCCGGTCATCATGGCCGGGCTGGCCATGGGCAGGACGATTTTCCACAGGGTCTCGGCCTTGGTTGCGCCCAACGCCAGCGAGCCTTCACGCACGGTGCGCGGAATTCGCGCCAGGCCTTCTTCGGTGGCCACGATCACCACCGGTACCGCCAGCAATGCCAGGGTCAGC  >  NC_002947/6072761‑6073056

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