Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,105,716 G→A 8.0% P815L (CCG→CTG)  ydiJ ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,105,7160GA 8.0% 213.3 / 7.6 76P815L (CCG→CTG) ydiJputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (30/39);  new base A (2/4);  total (33/43)
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.62e-01

GATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGCCGGCA  >  NC_002947/5105567‑5105863
                                                                                                                                                     |                                                                                                                                                   
gATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACg                                                                                                                                                     >  2:954424/1‑150 (MQ=255)
   acacCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGgcg                                                                                                                                                  >  1:1135081/1‑150 (MQ=255)
      ccAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCa                                                                                                                                               >  1:746319/1‑150 (MQ=255)
       cAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCGGGGCCTTGAACTGCGTGGCGTTGCGTATCGCTGCCTTGGCAAACGCGCCAAGGAAACACTGTACGTGCAACGACTTGCCGTTGGCGCTGTACCGCGcccg                                                                                                                                              >  2:376760/1‑148 (MQ=255)
       cAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGTGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGTTTAGTAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAg                                                                                                                                              <  1:822089/150‑1 (MQ=255)
       cAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAg                                                                                                                                              <  1:126970/150‑1 (MQ=255)
       cAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaa                                                                                                                                              <  2:636210/150‑2 (MQ=255)
        aGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTCGGCGCCGTACAGCGCCaac                                                                                                                                             <  1:1139496/150‑3 (MQ=255)
              aTCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATAcc                                                                                                                                       >  1:651294/1‑150 (MQ=255)
                  ccGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGt                                                                                                                                   <  1:1044208/150‑1 (MQ=255)
                  ccGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGt                                                                                                                                   <  1:984141/150‑1 (MQ=255)
                         cAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaa                                                                                                                                              <  1:127232/132‑2 (MQ=255)
                         cAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaa                                                                                                                                              >  2:127232/1‑131 (MQ=255)
                             cccACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCgg                                                                                                                        <  1:806334/150‑1 (MQ=255)
                              ccaccaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGt                                                                                                                       <  1:669396/150‑1 (MQ=255)
                               cctcaGCGGCCAACCGCCACCCGCCCGTGCTTCGAGCTGCCGGGTGATTCGGATTGCTGCCGTGGCAACCGCGCCAAGGGCACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGACCGCAATCCCCGGTGTGCCAGGCGGTg                                                                                                                      <  2:315688/147‑1 (MQ=255)
                                tggtggggCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTgg                                                                                                                     <  2:750981/144‑1 (MQ=255)
                                accaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTgg                                                                                                                     <  1:968990/150‑1 (MQ=255)
                                accaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTgg                                                                                                                     <  2:651294/150‑1 (MQ=255)
                                accaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTgg                                                                                                                     >  1:125788/1‑150 (MQ=255)
                                   aGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCa                                                                                                                  >  2:1068471/1‑150 (MQ=255)
                                    gCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCaa                                                                                                                 <  1:317194/150‑1 (MQ=255)
                                    gCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCaa                                                                                                                 <  1:1034755/150‑1 (MQ=255)
                                        cacaCCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTc                                                                                                             <  2:526269/150‑1 (MQ=255)
                                         acacCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCg                                                                                                            >  2:352614/1‑150 (MQ=255)
                                         acacCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCg                                                                                                            <  2:690780/150‑1 (MQ=255)
                                           acCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGAt                                                                                                          <  2:97838/150‑1 (MQ=255)
                                            ccGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATg                                                                                                         <  1:967455/150‑1 (MQ=255)
                                               cAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGACTACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAg                                                                                                      <  1:978606/150‑1 (MQ=255)
                                               cAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAg                                                                                                      <  2:606048/150‑1 (MQ=255)
                                                 aTCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGc                                                                                                    <  2:949489/150‑1 (MQ=255)
                                                     gCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGAc                                                                                                <  1:954424/150‑1 (MQ=255)
                                                     gCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGAc                                                                                                <  1:310065/150‑1 (MQ=255)
                                                      ccAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGAca                                                                                               >  2:294809/1‑150 (MQ=255)
                                                      ccAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGCTGGCCCAGGCGGTCGGCAAGGTCGCTGCAGGCCGAcc                                                                                               >  2:123037/1‑149 (MQ=255)
                                                        aGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGACAACGACAACATGTACGTGCAACGGGTTACCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGAcagc                                                                                             >  1:450403/1‑150 (MQ=255)
                                                          gcccTTGAGCGGCGTGGCATTGCGGATCGCTACCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGCCGCTGGGCAAGGTCGCTGAAGGCCGAcagcag                                                                                           <  1:897496/148‑1 (MQ=255)
                                                           gCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACagcagc                                                                                          <  2:633042/150‑1 (MQ=255)
                                                             cTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACagcagc                                                                                          <  2:412688/148‑1 (MQ=255)
                                                             cTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACagcagc                                                                                          >  1:412688/1‑148 (MQ=255)
                                                              ttGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAgcggc                                                                                       <  1:1123555/150‑1 (MQ=255)
                                                                  gCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGttt                                                                                   <  2:1135081/150‑1 (MQ=255)
                                                                       gTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAg                                                                              <  2:414576/150‑1 (MQ=255)
                                                                             ttGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCgg                                                                        >  1:1089493/1‑150 (MQ=255)
                                                                                 ggATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGaa                                                                                                       <  2:942859/115‑1 (MQ=255)
                                                                                 ggATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGaa                                                                                                       >  1:942859/1‑115 (MQ=255)
                                                                                  gATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGaa                                                                   <  2:746319/150‑1 (MQ=255)
                                                                                     cGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGc                                                                <  2:125788/150‑1 (MQ=255)
                                                                                        tGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCAGACAGCAGAGGCGGTTAGAAGTAGCGGGTGAAGGCGtg                                                             >  1:252438/1‑149 (MQ=255)
                                                                                            ttGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGc                                                         >  1:136935/1‑150 (MQ=255)
                                                                                                    tttgcCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCAc                                                 <  1:407229/147‑1 (MQ=255)
                                                                                                       gcCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGAt                                              >  2:98391/1‑150 (MQ=255)
                                                                                                         cAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGAtg                                             <  2:542496/149‑1 (MQ=255)
                                                                                                         cAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGAtg                                             >  1:542496/1‑149 (MQ=255)
                                                                                                         cAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCAGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTc                                                             >  2:179840/1‑133 (MQ=255)
                                                                                                         cAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCAGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTc                                                             <  1:179840/133‑1 (MQ=255)
                                                                                                             gAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCg                                                                         >  1:409322/1‑117 (MQ=255)
                                                                                                             gAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCg                                                                         <  2:409322/117‑1 (MQ=255)
                                                                                                               aaCCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCt                                      >  1:613252/1‑150 (MQ=255)
                                                                                                                 cccTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTcg                                    >  1:221448/1‑150 (MQ=255)
                                                                                                                     gTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCt                                >  1:940276/1‑150 (MQ=255)
                                                                                                                         gTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGc                            <  1:45565/150‑1 (MQ=255)
                                                                                                                           gCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTc                                     >  1:106055/1‑139 (MQ=255)
                                                                                                                           gCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTc                                     <  2:106055/139‑1 (MQ=255)
                                                                                                                           ccAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGGCAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCg                          >  1:1199982/2‑150 (MQ=255)
                                                                                                                            cAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGGACGATGCTGCGTTCGCGCTGGGCCgg                         >  1:1107867/1‑150 (MQ=255)
                                                                                                                             aaCGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCggg                        <  2:940276/150‑1 (MQ=255)
                                                                                                                                ggCTTGCCGTTGGCGCTGTACGGAGCAAGGAATACACGGTGGACCAGGCGGAGGGCAAGGGCGAGGAAGACCGACAGCAGCGGCGTTTAGAAGTAGAGGGTGAAAGCGTCCTGCACCAACACAATGCTGCGCACGAGCTTGGGCggggac                     >  2:92748/1‑148 (MQ=255)
                                                                                                                                 gCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGACCGGGGTCa                    >  1:799987/1‑150 (MQ=255)
                                                                                                                                     gCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTc                >  1:255554/1‑150 (MQ=255)
                                                                                                                                         ctGGCGCCGTACAGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCgcgc            >  2:26633/2‑150 (MQ=255)
                                                                                                                                            gcgcTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAgg         <  1:352614/150‑1 (MQ=255)
                                                                                                                                            gcgcTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGGTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAgg         >  2:290769/1‑150 (MQ=255)
                                                                                                                                              gcTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGc       <  1:98391/150‑1 (MQ=255)
                                                                                                                                                cgTACAGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGCCg     <  1:26633/149‑1 (MQ=255)
                                                                                                                                                  tACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGACGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGAGCCGGc   >  2:508151/1‑150 (MQ=255)
                                                                                                                                                   aCGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGCCGGCa  >  1:252484/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGCCGGCA  >  NC_002947/5105567‑5105863

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