Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,105,711 T→C 6.5% S817G (AGC→GGC)  ydiJ ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,105,7110TC 6.5% 226.6 / 4.5 77S817G (AGC→GGC) ydiJputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (32/40);  new base C (2/3);  total (34/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 = 5.94e-01

TGGCGATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGC  >  NC_002947/5105563‑5105858
                                                                                                                                                    |                                                                                                                                                   
tGGCGATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTg                                                                                                                                                    >  2:984141/1‑150 (MQ=255)
    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)
            aGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTa                                                                                                                                                  <  2:942220/140‑1 (MQ=255)
            aGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTa                                                                                                                                                  >  1:942220/1‑140 (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                                                                                                                                         >  2:127232/1‑131 (MQ=255)
                             cAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaa                                                                                                                                         <  1:127232/132‑2 (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                                                                                                                >  1:125788/1‑150 (MQ=255)
                                    accaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTgg                                                                                                                <  2:651294/150‑1 (MQ=255)
                                       aGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCa                                                                                                             >  2:1068471/1‑150 (MQ=255)
                                        gCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCaa                                                                                                            <  1:1034755/150‑1 (MQ=255)
                                        gCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCaa                                                                                                            <  1:317194/150‑1 (MQ=255)
                                            cacaCCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTc                                                                                                        <  2:526269/150‑1 (MQ=255)
                                             acacCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCg                                                                                                       <  2:690780/150‑1 (MQ=255)
                                             acacCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCg                                                                                                       >  2:352614/1‑150 (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                                                                                     >  1:412688/1‑148 (MQ=255)
                                                                 cTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACagcagc                                                                                     <  2:412688/148‑1 (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                                                                                                  >  1:942859/1‑115 (MQ=255)
                                                                                     ggATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGaa                                                                                                  <  2:942859/115‑1 (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)
                                                                                                                                                    |                                                                                                                                                   
TGGCGATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGC  >  NC_002947/5105563‑5105858

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