Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,260,670 C→G 7.1% R712R (CGG→CGC)  PP_3734 ← ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,260,6700CG 7.1% 196.2 / 5.2 65R712R (CGG→CGC) PP_3734ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (37/23);  new base G (3/2);  total (40/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TGCGACCAGGCACCAGGCCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGACCACCCGGCT  >  NC_002947/4260523‑4260819
                                                                                                                                                   |                                                                                                                                                     
tGCGACCAGGCACCAGGCCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCcg                                                                                                                                                     >  2:188489/1‑150 (MQ=255)
       aggcaccaggcCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTgg                                                                                                                                              <  2:775763/150‑1 (MQ=255)
        ggcaccaggcCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGc                                                                                                                                             >  1:491995/1‑150 (MQ=255)
                      aaCAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAgg                                                                                                                               >  2:303261/1‑150 (MQ=255)
                       aCAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGt                                                                                                                              >  1:492178/1‑150 (MQ=255)
                        cAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCTCCGCTGCGCCTGGCCCAGGGTCAGCAGGtt                                                                                                                             >  1:73/1‑150 (MQ=255)
                                    aagcGCCAGCAGTACGGTGAAGGCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGc                                                                                                                 <  1:346103/147‑1 (MQ=255)
                                     aTCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCAAAGCGGGGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGcc                                                                                                                >  1:12309/1‑150 (MQ=255)
                                        gCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGATGTCGCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCcaca                                                                                                             >  1:622674/1‑150 (MQ=255)
                                         ccAGCAGTAGGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCacac                                                                                                            >  2:494196/1‑150 (MQ=255)
                                         ccAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCacac                                                                                                            >  1:1140039/1‑150 (MQ=255)
                                          cagcagTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACAcc                                                                                                           >  2:1049246/1‑150 (MQ=255)
                                            gcagTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCgg                                                                                                         <  1:722340/150‑1 (MQ=255)
                                                tACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCTCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGCCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCAc                                                                                                     >  1:726367/1‑150 (MQ=255)
                                                     tGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCATTCGCATGAGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACCCCGGCCACGCCaa                                                                                                >  2:441154/1‑150 (MQ=255)
                                                       aaGCCGCGGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAgg                                                                                              <  2:219600/150‑1 (MQ=255)
                                                           cGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTaa                                                                                          >  2:615686/1‑150 (MQ=255)
                                                            gCTGAGCATGAGCGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAg                                                                                         >  2:318381/1‑150 (MQ=255)
                                                                  cATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGtt                                                                                   <  2:492178/150‑1 (MQ=255)
                                                                   aTGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTg                                                                                  >  1:692050/1‑150 (MQ=255)
                                                                    tGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTga                                                                                 <  2:692050/150‑1 (MQ=255)
                                                                          tCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTg                                                                           >  1:708870/1‑150 (MQ=255)
                                                                          ccTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTg                                                                           >  2:255043/2‑150 (MQ=255)
                                                                              ccccAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGt                                                                       >  2:986442/1‑150 (MQ=255)
                                                                                 cAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTggc                                                                    >  2:100520/1‑150 (MQ=255)
                                                                                   ggCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTCCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTggcgg                                                                  <  1:494196/150‑1 (MQ=255)
                                                                                        aaCCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAg                                                             <  1:1049246/150‑1 (MQ=255)
                                                                                         accagaccagTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAgg                                                            <  1:318381/150‑1 (MQ=255)
                                                                                             gaccagTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCg                                                        <  2:708870/150‑1 (MQ=255)
                                                                                              accagTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCgg                                                       >  1:153611/1‑150 (MQ=255)
                                                                                              accagTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCgg                                                       <  2:1140039/150‑1 (MQ=255)
                                                                                                cagTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGtt                                                     <  1:188489/150‑1 (MQ=255)
                                                                                                       aTGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCCCGCCAAGGGTAAGGCCGTTGAGAGCTGCGGTGGCCGCAAAGGTGCGGTTGAAGACg                                              >  2:199993/1‑150 (MQ=255)
                                                                                                       aTGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTATGGCTGtt                                                                                   <  1:1144888/113‑1 (MQ=255)
                                                                                                       aTGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTATGGCTGtt                                                                                   >  2:1144977/1‑113 (MQ=255)
                                                                                                       aTGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTATGGCTGtt                                                                                   >  2:1144888/1‑113 (MQ=255)
                                                                                                       aTGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTATGGCTGtt                                                                                   <  1:1144977/113‑1 (MQ=255)
                                                                                                       aTGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCg                                                                         >  1:129237/1‑123 (MQ=255)
                                                                                                       aTGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCg                                                                         <  2:129237/123‑1 (MQ=255)
                                                                                                       aTGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACg                                              >  2:1251239/1‑150 (MQ=255)
                                                                                                           gccgcACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCgg                                                       >  1:250827/1‑137 (MQ=255)
                                                                                                           gccgcACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCgg                                                       <  2:250827/137‑1 (MQ=255)
                                                                                                                aCCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGAc                                     >  2:310769/1‑150 (MQ=255)
                                                                                                                  ccccAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCa                                   >  1:625472/1‑150 (MQ=255)
                                                                                                                    ccAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTgag                                                                                >  2:826316/1‑103 (MQ=255)
                                                                                                                    ccAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTgag                                                                                <  1:826316/103‑1 (MQ=255)
                                                                                                                           ccACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAg                          >  2:858980/1‑150 (MQ=255)
                                                                                                                            cacaGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAg                                                             >  1:63058/1‑114 (MQ=255)
                                                                                                                            cacaGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAg                                                             <  2:63058/114‑1 (MQ=255)
                                                                                                                               aGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTCCAGgcgc                      <  2:961630/150‑1 (MQ=255)
                                                                                                                                 cGGTGCTAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGgcgcgc                    <  1:441154/150‑1 (MQ=255)
                                                                                                                                     gCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGc                <  1:100520/150‑1 (MQ=255)
                                                                                                                                       cAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTc              >  1:484828/1‑150 (MQ=255)
                                                                                                                                       cAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTc              <  1:303261/150‑1 (MQ=255)
                                                                                                                                       cAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGAAGGCGGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTc              <  2:622674/150‑1 (MQ=255)
                                                                                                                                        aaCTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCg             >  2:1062347/1‑150 (MQ=255)
                                                                                                                                         ccTGGGCCAGGCGCGGCTGGGCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTc                            <  2:1227709/133‑1 (MQ=255)
                                                                                                                                         ccTGGGCCAGGCGCGGCTGGGCCAGGGTCAGCAGGTTGATCAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTc                            >  1:1227709/2‑134 (MQ=255)
                                                                                                                                         aCTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGa            >  2:1214695/1‑150 (MQ=255)
                                                                                                                                           tgagccAGGCGCGGCTGGGCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAAACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGacc          >  1:593761/5‑150 (MQ=255)
                                                                                                                                           tGGGCCAGGCGCGGCTGGGCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGa                                      >  2:413438/1‑122 (MQ=255)
                                                                                                                                           tGGGCCAGGCGCGGCTGGGCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCGGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGa                                      <  1:413438/122‑1 (MQ=255)
                                                                                                                                               ccAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGACAGCACAGATGGGGTGAAAGGTGCGGTTGAAGAAGTAGGTCGCCCGGTTTGTCAGGCGTGCCTGCTCGACCAccc      >  1:348174/1‑150 (MQ=255)
                                                                                                                                                  gCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGACCACCCGGc   >  2:1097528/1‑150 (MQ=255)
                                                                                                                                                  gCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGACCACCCGGc   >  2:1111896/1‑150 (MQ=255)
                                                                                                                                                   ccGCGCGTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCAAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGACCACCCGGCt  <  2:491995/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGCGACCAGGCACCAGGCCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGACCACCCGGCT  >  NC_002947/4260523‑4260819

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑