Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,290,379 T→A 7.9% intergenic (‑90/+94) PP_4660 ← / ← mdaB LysR family transcriptional regulator/drug activity modulator B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,290,3790TA 7.9% 407.3 / 22.1 150intergenic (‑90/+94)PP_4660/mdaBLysR family transcriptional regulator/drug activity modulator B
Reads supporting (aligned to +/- strand):  ref base T (58/69);  new base A (4/8);  total (69/81)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

Only 100 of 151 total aligned reads displayed.

AATCGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAACGCCG  >  NC_002947/5290230‑5290528
                                                                                                                                                     |                                                                                                                                                       
tgggcatgaccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCt                                                                                                                                                         >  2:3709809‑M1/135‑150 (MQ=255)
aaTCGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCT                                                                                                                                                         >  1:3126391‑M1/1‑150 (MQ=255)
 ctcGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTT                                                                                                                                                        <  2:451623‑M1/149‑1 (MQ=255)
 aTCGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATACCCGCTTGTCGAGGCCAAGGCCAGCGGCGCCTGTATGGGCCTCGTCGCGGGCTT                                                                                                                                                        <  2:450618‑M1/150‑1 (MQ=255)
  ggcatgaccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTg                                                                                                                                                       <  2:2768213‑M1/18‑1 (MQ=255)
         gCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTC                                                                                                                                                <  2:2858362‑M1/150‑1 (MQ=255)
         gCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTC                                                                                                                                                <  1:2120360‑M1/150‑1 (MQ=255)
         gCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGAGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGCga                                                                                                                                                <  2:4372706‑M1/150‑3 (MQ=255)
         ccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTc                                                                                                                                                >  2:3903681‑M1/126‑150 (MQ=255)
         ccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTc                                                                                                                                                >  2:1846592‑M1/126‑150 (MQ=255)
         ccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTc                                                                                                                                                <  2:784166‑M1/25‑1 (MQ=255)
          cggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTcc                                                                                                                                               >  2:1104595‑M1/125‑150 (MQ=255)
          cggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                               >  1:1403498‑M1/125‑147 (MQ=255)
          cggcccgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                               >  2:3303751‑M1/125‑147 (MQ=255)
          cAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCC                                                                                                                                               <  1:4131291‑M1/150‑1 (MQ=255)
          cAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGCgaa                                                                                                                                               <  1:1333514‑M1/150‑4 (MQ=255)
           ggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTccc                                                                                                                                              >  1:576870‑M1/124‑150 (MQ=255)
           aGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGCgaag                                                                                                                                              <  1:2131645‑M1/150‑5 (MQ=255)
            gcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCca                                                                                                                                             <  2:3126391‑M1/28‑1 (MQ=255)
            gcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCca                                                                                                                                             >  2:2679201‑M1/123‑150 (MQ=255)
            gTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGCgaaga                                                                                                                                             <  2:1663446‑M1/150‑6 (MQ=255)
               gATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACaa                                                                                                                                          <  2:1387397/150‑1 (MQ=255)
                  gACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTa                                                                                                                                       <  1:614997/150‑1 (MQ=255)
                   aCGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTAc                                                                                                                                      <  1:543537/150‑1 (MQ=255)
                    cGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACt                                                                                                                                     <  1:1626043/150‑1 (MQ=255)
                      gCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGc                                                                                                                                   <  2:3782396/150‑1 (MQ=255)
                       cgacgaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCa                                                                                                                                  <  1:434580/150‑1 (MQ=255)
                         acgaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGACTGTGGTTGGTCTGTGCAATGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgca                                                                                                                                <  2:4200920‑M1/150‑12 (MQ=255)
                          cgaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAg                                                                                                                               <  2:258410/150‑1 (MQ=255)
                           gaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTGCCCGCTCCCACAAGTACTGCAAAGt                                                                                                                              >  2:2105241/1‑150 (MQ=255)
                              tACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCt                                                                                                                           >  2:2466490/1‑150 (MQ=255)
                               aCCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCtt                                                                                                                          <  2:338842/150‑1 (MQ=255)
                               aCCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCtt                                                                                                                          <  2:2921223/150‑1 (MQ=255)
                               aCCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCtt                                                                                                                          <  2:1243420/150‑1 (MQ=255)
                                ccTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTg                                                                                                                         <  1:1882528/150‑1 (MQ=255)
                                      gTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCCCACCCACAAGTAATGCAAAGTCCTTGAAAccg                                                                                                                   >  1:1098658/1‑149 (MQ=255)
                                       ttcttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTg                                                                                                                  <  1:3242978/150‑1 (MQ=255)
                                            cGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTaa                                                                                                             <  1:2149304/150‑1 (MQ=255)
                                              gATCGGGTTTTCACCTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATa                                                                                                           <  1:977128/150‑1 (MQ=255)
                                               aTCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATAc                                                                                                          <  1:1707291/150‑1 (MQ=255)
                                                 cGGGTTTTCACATCAGTGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCt                                                                                                        >  1:2913963‑M2/1‑149 (MQ=255)
                                                  gggTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGGCCAGCGGCGCCTGTCTGGGCCTCTTCGCGGGCGTGCCCGCTCCCACAACTACATCCAAGTCCTGGTAACCGGCGTCATAccag                                                                                                       >  2:1235900/1‑148 (MQ=255)
                                                   ggTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgt                                                                                                      >  1:1082106‑M2/1‑147 (MQ=255)
                                                   ggTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgt                                                                                                      >  1:265090‑M2/1‑147 (MQ=255)
                                                    gTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                     >  2:599869‑M2/1‑146 (MQ=255)
                                                    gTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                     >  2:1808852‑M2/1‑146 (MQ=255)
                                                     ttttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgg                                                                                                    >  2:2420696‑M2/1‑145 (MQ=255)
                                                     ttttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgg                                                                                                    >  1:818401‑M2/1‑145 (MQ=255)
                                                     ttttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgg                                                                                                    >  1:851191‑M2/1‑145 (MQ=255)
                                                     ttttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCagtgg                                                                                                    >  1:3591882‑M2/1‑145 (MQ=255)
                                                     ttttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTTCCCCCGCCCCCAACTACAGCAAAGTACGTGAAACCGGCGGAATAACagtgg                                                                                                    <  1:3989597‑M2/150‑6 (MQ=255)
                                                       ttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtggga                                                                                                  <  2:1893697‑M1/150‑42 (MQ=255)
                                                          acaTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcg                                                                                               <  2:2178828‑M1/150‑45 (MQ=255)
                                                          acaTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcg                                                                                               <  2:2179076‑M1/150‑45 (MQ=255)
                                                          acaTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactctgcagtacttgtgggagcg                                                                                               <  1:1171017‑M1/150‑45 (MQ=255)
                                                          acaTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAAACTGCGTAATACCggtgggggcg                                                                                               <  2:2091487‑M2/150‑11 (MQ=255)
                                                           ccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCgg                                                                                            >  2:2676050‑M1/76‑150 (MQ=255)
                                                           ccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCgg                                                                                            >  2:1683001‑M1/76‑150 (MQ=255)
                                                           caTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcgg                                                                                              <  1:1779042‑M1/150‑46 (MQ=255)
                                                           caTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcgg                                                                                              <  1:3942985‑M1/150‑46 (MQ=255)
                                                           caTCAGGGCCTGGGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcgg                                                                                              <  2:2338870‑M1/150‑46 (MQ=255)
                                                            cggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCggg                                                                                           >  2:3596921‑M1/75‑150 (MQ=255)
                                                            aTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGc                                                                                                                                   <  2:80519/112‑1 (MQ=255)
                                                            aTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGc                                                                                                                                   >  1:80519/1‑112 (MQ=255)
                                                             tCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGGTTGCCCGCTACCACAAGTACTGCAAAGTCCTTTTAAACTGGGTAATACCtgtgggtgcgggc                                                                                            <  2:3480756‑M2/150‑14 (MQ=255)
                                                             tCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTTCCCGCTCCCACAACTACTGCAAAGTCCTTGAAACCTGCGGAATACCtgtgggagcgggc                                                                                            <  2:3879770‑M2/150‑14 (MQ=255)
                                                             tCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAAAAACtgtgggagcgggc                                                                                            <  2:2127008‑M2/150‑14 (MQ=255)
                                                             tCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACAGCACAGTCCTTGAAACCGGCGGAATACCagtgggagcgggc                                                                                            <  2:4194116‑M2/150‑14 (MQ=255)
                                                             tCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCGCCCCCAAGGACTGCAAAGTCCGTGAAACCTGCGTAATACCtgggggagcgggc                                                                                            <  2:3880444‑M2/150‑14 (MQ=255)
                                                             tCAGGGCCTGTGGTTGGTCTGTGCAAGGATGAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGTGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtactcgtgggagcgggc                                                                                            <  2:4085034‑M1/150‑48 (MQ=255)
                                                             ggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGc                                                                                          >  2:1810562‑M1/74‑150 (MQ=255)
                                                                  gCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgggagcgggcaagcc                                                                                       >  1:2338747‑M2/1‑132 (MQ=255)
                                                                    cTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATa                                                                                                           >  2:4378792/1‑128 (MQ=255)
                                                                    cTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCGGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATa                                                                                                           <  1:4378792/128‑1 (MQ=255)
                                                                          ttggtCTGTGCACGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcgggcaagcccgcgaaga                                                                               <  2:4444493‑M1/150‑61 (MQ=255)
                                                                          ttggtCTGTGCAAGGATTAGAGACCCGCTGGTCGCGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTCCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAACCCGGCGTAATACCtgtgggagcgggcaagcccgcgaaga                                                                               >  2:4248496‑M2/1‑124 (MQ=255)
                                                                             gtCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTAAGGGCCTCTTCGCGGGCGTGCCCCCTCCCACCAGAAAGGCCAAATGCCTGAAAACCCCGGAATAACAGTGGGGGCCGGCCAAGCCCCGAATAGGc                                                                          >  1:1404902/1‑150 (MQ=255)
                                                                                    cAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCAAAAATACTGCAAAGTCCTTGAAACCTGCGGAATACCTGTGGGAGCGGGCAAACCCGCGAAGAGGCCGGTACa                                                                   >  2:440587/1‑150 (MQ=255)
                                                                                    cAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCGCCCCCCAGGACGGCCAACTCCGTGGAACCGGCGTAAAAACTGTGTGGGCCGGCCAGCCCCCCAAGAGGCCggggca                                                                   >  1:2741862/1‑146 (MQ=255)
                                                                                      ctgATTAGATCACCGCTTGTCTCGGCCACGGCCAGCGGCGCCTGTTTGGGCCTCTCGGGGGGCTTCCCGGCTCCCAAACTTACTGAAAAGTCCTTGAAACCTGCTAATTCCGGTTGGTTGCGGGAAAGCCCGCGAAGAGGCCGGTACAGa                                                                 <  2:1662140/148‑1 (MQ=255)
                                                                                             gATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTTCCCGCTGCCACAAGTACTGCAAAGTCCTTGTAACCTGGGTAATACCTGTGGGTGCGGGGAAGCCCGCGTAGAGGCCGGTACAGACCTgaaa                                                          >  2:3593565/1‑150 (MQ=255)
                                                                                                         tCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGTCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACCTAAAATTAATGTGCCtt                                              >  2:612384/1‑150 (MQ=255)
                                                                                                         tCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGTCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCtt                                              >  2:612445/1‑150 (MQ=255)
                                                                                                          cGAGGCCAAGTCCAGCGGCGCCTGTAGGGGCCTCGACGCGGGCTGGCCCGCTCCCACACGTACTGCAACGTCATTGAACCCGGCTTAATACCTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTg                                             <  2:3256639/150‑1 (MQ=38)
                                                                                                                       aGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCCCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGGAAAACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAAATACTGTGCCGTGCCAAACACCTTGt                                >  2:1101353/1‑150 (MQ=37)
                                                                                                                        gcggcgCCTGTATGGGCCTCTTGGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTc                               <  2:3179066/150‑1 (MQ=255)
                                                                                                                             gcgtGTTTGGGCCTCTTCGGGGGGTTGCCCGCACCCACACTTACGGCAACGCCCTTGAAACCTGCTTAATACCGGTGGGAGGGGGAAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGAt                          <  2:390386/147‑1 (MQ=255)
                                                                                                                               gtgttggggCCTGTGGGCGGGGTTGCCCGCTCACACACTTACTCAAATGTCTTTTAAACCTGCGATCTTCCTGTGGTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATgg                        <  2:1703436/144‑1 (MQ=255)
                                                                                                                                tGTATGGGCCTCTTCGCGGGCTTGCCCCCTCCCACAAGTACTGCAAACTACTTGAAACCAGCGGAAAACCGGTGGGAGCGGGCAAGCCCCCGAAGAGGCCGGTACAGACATAAAAATAATGTGCCATTCCAAACACCTTGTCCAGATtgg                       >  2:261467/1‑150 (MQ=25)
                                                                                                                                          tctTCGCGGGCTTGCCCGCTCCCACAggtattacgccgggttcaaggactttgcagtacttgtgggagcgggcaagcccgcggagaggcccatacaggcgccgctggacttggccgcgaccagcgggtctctaatccttgcacagaccaa               >  2:4248496‑M1/1‑26 (MQ=255)
                                                                                                                                           CTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTcataatct            <  2:3057280‑M1/150‑4 (MQ=255)
                                                                                                                                            TTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGCGCCTTGCCAAACACCTTGTCCAGATGGGCCGCATAAGCTg           <  2:2509481‑M1/150‑1 (MQ=255)
                                                                                                                                             TCGCGGGCTTGCCCGCTCCAACAAGTACTCCAATGTCCTTGAACCCGGCTTAATACCTGTGGTAGCGGGCACGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTAGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGc          >  2:531599‑M1/1‑150 (MQ=255)
                                                                                                                                                cgGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                     >  1:1634764‑M2/1‑54 (MQ=255)
                                                                                                                                                cgGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                     <  2:1634764‑M2/58‑5 (MQ=255)
                                                                                                                                                 GGGCTTCCCCGCCCCCACATTTACTCCAAAGTCCTTCAAACCTGCTTAATACCGGGGGGCGCGGGCACGCCCGCTAAGAGGCCGGTACTGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAAc      >  1:290643‑M1/1‑150 (MQ=255)
                                                                                                                                                  ggCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcgggcaagcccgcgaagaggcccatacaggcgccgctggacttggcctcgacaagcggttatctaatccttgcacagaccaaccacaggc       >  1:2338747‑M1/1‑18 (MQ=255)
                                                                                                                                                   gCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagggggcaagcccgcgaagaggcccatacaggcgccgctggacttggcctcgacaagcggttatctaatccttgcacagaccaaccacaggcc      <  2:982049‑M1/150‑134 (MQ=255)
                                                                                                                                                     TTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAAcgccg  <  2:4115316‑M1/150‑1 (MQ=255)
                                                                                                                                                     TTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAAcgccg  <  1:592017‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AATCGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAACGCCG  >  NC_002947/5290230‑5290528

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