Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,260,680 C→G 7.5% G709A (GGC→GCC)  PP_3734 ← ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,260,6800CG 7.5% 205.1 / 7.9 67G709A (GGC→GCC) PP_3734ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (37/25);  new base G (3/2);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

GGCACCAGGCCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGACCACCCGGCT  >  NC_002947/4260531‑4260819
                                                                                                                                                     |                                                                                                                                           
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                                                       <  2:250827/137‑1 (MQ=255)
                                                                                                   gccgcACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCgg                                                       >  1:250827/1‑137 (MQ=255)
                                                                                                        aCCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGAc                                     >  2:310769/1‑150 (MQ=255)
                                                                                                          ccccAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCa                                   >  1:625472/1‑150 (MQ=255)
                                                                                                            ccAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTgag                                                                                <  1:826316/103‑1 (MQ=255)
                                                                                                            ccAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTgag                                                                                >  2:826316/1‑103 (MQ=255)
                                                                                                                   ccACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAg                          >  2:858980/1‑150 (MQ=255)
                                                                                                                    cacaGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAg                                                             <  2:63058/114‑1 (MQ=255)
                                                                                                                    cacaGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAg                                                             >  1:63058/1‑114 (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)
                                                                                                                                            cgcgCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGaa                                                  <  2:543719/101‑1 (MQ=255)
                                                                                                                                            cgcgCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGaa                                                  <  2:543698/101‑1 (MQ=255)
                                                                                                                                            cgcgCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGaa                                                  >  1:543698/1‑101 (MQ=255)
                                                                                                                                            cgcgCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGaa                                                  >  1:543719/1‑101 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GGCACCAGGCCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTGGCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGACCACCCGGCT  >  NC_002947/4260531‑4260819

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