Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,614,321 A→G 5.0% intergenic (+61/+359) yjjV → / ← yjjW putative DNase/putative pyruvate formate lyase activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,614,3210AG 5.0% 1531.9 / 32.1 474intergenic (+61/+359)yjjV/yjjWputative DNase/putative pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base A (164/275);  new base G (17/7);  total (188/286)
Fisher's exact test for biased strand distribution p-value = 1.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 480 total aligned reads displayed.

CCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATA  >  NC_000913/4614172‑4614470
                                                                                                                                                     |                                                                                                                                                     
          cGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTa                                                                                                                                             <  1:621232/150‑1 (MQ=255)
          cGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTa                                                                                                                                             <  2:5595893/150‑1 (MQ=255)
           gAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTAt                                                                                                                                            <  2:2942832/150‑1 (MQ=255)
           gAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTAt                                                                                                                                            <  1:5950454/150‑1 (MQ=255)
              ttGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCAt                                                                                                                                         >  1:906603/1‑150 (MQ=255)
                gcgcCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGc                                                                                                                                       <  1:422204/150‑1 (MQ=255)
                gcgcCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGc                                                                                                                                       <  2:3201564/150‑1 (MQ=255)
                 cgcCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGcc                                                                                                                                      >  2:2347688/1‑150 (MQ=255)
                  gcCGGGAACCGGCGGATGATATTGCTCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACTTGCCGTAGTCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCTAACAATGCCTGATAAGACGCGTTAAGAGTCTCATCAGGcat                                                                                                                                     <  1:5560564/150‑3 (MQ=25)
                  gcCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCt                                                                                                                                     >  1:3084386/1‑150 (MQ=255)
                   ccGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTa                                                                                                                                    <  1:5039910/150‑1 (MQ=255)
                       gAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAgg                                                                                                                                <  2:1579822/150‑1 (MQ=255)
                        aaCCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGt                                                                                                                               <  1:4945444/150‑1 (MQ=255)
                          ccGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTAGGCGCAAAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGttt                                                                                                                             >  2:3489533/1‑150 (MQ=255)
                           cggcggATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCAACGCTTGTCGCGTCTTATCATGCCTACAGGTTtg                                                                                                                            >  1:2151421/1‑150 (MQ=255)
                             gcggATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTtgtg                                                                                                                          <  2:541621/150‑1 (MQ=255)
                              cggATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGc                                                                                                                         <  2:5716402/150‑1 (MQ=255)
                              cggATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGc                                                                                                                         <  2:5168353/150‑1 (MQ=255)
                              cggATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGATTAACGAGCCGTAGGCCGGATAAGGCGTACACGCCGCATCCGGCAGTTGGCGCACAATGCCAGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGc                                                                                                                         <  2:460508/150‑1 (MQ=255)
                               ggATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGcc                                                                                                                        <  1:2285994/150‑1 (MQ=255)
                               ggATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGcc                                                                                                                        <  1:1104949/150‑1 (MQ=255)
                               ggATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGcc                                                                                                                        <  2:2989658/150‑1 (MQ=255)
                                 aTGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTCTGTGCCGa                                                                                                                      >  1:749418/1‑150 (MQ=255)
                                  tGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGaa                                                                                                                     >  1:5821802/1‑150 (MQ=255)
                                      aTTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGt                                                                                                                 >  2:4681743/1‑150 (MQ=255)
                                          cgcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGc                                                                                                             <  1:5439858/150‑1 (MQ=255)
                                            cAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCg                                                                                                           <  1:1497183/150‑1 (MQ=255)
                                              aGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGa                                                                                                         <  2:1438795/150‑1 (MQ=255)
                                               gCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGAt                                                                                                        <  1:5693645/150‑1 (MQ=255)
                                                 gTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATaa                                                                                                      <  2:5861166/150‑1 (MQ=255)
                                                  ttGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGAGGCTTAACGCGTCTTATCATGCCTACTGGTTTGTGCCGAACCGTAGGCCGGATAAg                                                                                                     >  2:669448/1‑150 (MQ=255)
                                                  ttGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAg                                                                                                     <  1:4148199/150‑1 (MQ=255)
                                                   tGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGaa                                                                                                                     <  1:1581780/133‑1 (MQ=255)
                                                    gCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGc                                                                                                   <  2:2786758/150‑1 (MQ=255)
                                                     cTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCg                                                                                                  >  1:3519055/1‑150 (MQ=255)
                                                      ttAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGt                                                                                                 >  1:276775/1‑150 (MQ=255)
                                                       taataaCACGTATACGGTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGtt                                                                                                <  1:2884916/150‑1 (MQ=255)
                                                           aaCACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcg                                                                                            <  2:5361102/150‑1 (MQ=255)
                                                             cacGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgcc                                                                                          <  2:1097597/150‑1 (MQ=255)
                                                             cacGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgcc                                                                                          <  2:5956650/150‑1 (MQ=255)
                                                             cacGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgcc                                                                                          <  1:1040258/150‑1 (MQ=255)
                                                              acGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgccg                                                                                         >  1:4118091/1‑150 (MQ=255)
                                                               cGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                        <  1:5214897/150‑1 (MQ=255)
                                                                 tataCGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTCACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCTCGCCGCAt                                                                                      <  2:552723/150‑1 (MQ=255)
                                                                    aCGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                   <  2:790128/150‑1 (MQ=255)
                                                                     cGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                  <  2:1811290/150‑1 (MQ=255)
                                                                      gttgttTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGACCCGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                 <  2:2006934/150‑1 (MQ=17)
                                                                      gttgttTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGc                                                                                 >  1:5103763/1‑150 (MQ=255)
                                                                      gttgttTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                 >  1:2067852/1‑150 (MQ=255)
                                                                      gttgttTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                 >  1:3707405/1‑150 (MQ=255)
                                                                       ttgttTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                <  2:1015783/150‑1 (MQ=255)
                                                                       ttgttTAAATTGCCGAAGGCCGGATAAGGCGTTCACGCAGCATACGGCAGATGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCATACAGGTTTGAGCCGAACCGTAGGCCGGATAAGGCGATCACGCCGCATCCGGCa                                                                                >  2:2716624/1‑150 (MQ=14)
                                                                        tgttTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTCCAGGTATGTGCCGAACCGAAGGCCGGATACGGCGTTCACGCCGCATCCGGCAg                                                                               >  1:4374611/1‑150 (MQ=255)
                                                                        tgttTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCTCGCCGCATCCGGCAg                                                                               >  1:5240992/1‑150 (MQ=255)
                                                                        tgttTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAg                                                                               <  2:1711115/150‑1 (MQ=255)
                                                                         gttTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTCACGCGTCTTATCATGCCTACAGGGAAGAGCCGCACCGTAGGCCGCATAAGGCGTTCGCGCCGCATCCGGCCGt                                                                              >  1:2600971/1‑150 (MQ=14)
                                                                           ttAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTg                                                                            <  1:3033136/150‑1 (MQ=255)
                                                                             aaCGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgc                                                                          <  1:1004391/150‑1 (MQ=255)
                                                                             aaCGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgc                                                                          <  1:5503854/150‑1 (MQ=255)
                                                                              aCGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                         <  1:4661141/150‑1 (MQ=255)
                                                                                gTGCCGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTCACGCGTCTTATCATGCCTACAGGTGTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTGGGCGca                                                                       >  1:4546422/1‑150 (MQ=255)
                                                                                gTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGca                                                                       >  2:1228712/1‑150 (MQ=255)
                                                                                   ccGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACaa                                                                    >  2:3414303/1‑150 (MQ=255)
                                                                                     ggaGGCCGGAGCAGGCGTTCCCGCCGCCTCCGGCAGGTGGCGCACACTGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATg                                                                  <  2:165309/148‑1 (MQ=14)
                                                                                      tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGc                                                                 <  2:1035173/150‑1 (MQ=255)
                                                                                       aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGcc                                                                <  1:5316715/150‑1 (MQ=255)
                                                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCt                                                               <  1:2473894/150‑1 (MQ=255)
                                                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTGCCCGCAGCATCCGGCAGGTGCAGGACACTGCCt                                                               >  2:3358966/1‑150 (MQ=11)
                                                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTg                                                              >  2:2200918/1‑150 (MQ=255)
                                                                                          ccGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCtt                                                                                                                                              <  1:5378134/69‑1 (MQ=40)
                                                                                           cGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATCAGGCGTGCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGAt                                                            >  2:3697639/1‑150 (MQ=255)
                                                                                            ggATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATg                                                           <  2:1934936/150‑1 (MQ=255)
                                                                                             gATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGc                                                          <  1:2439709/150‑1 (MQ=255)
                                                                                               tAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGa                                                        <  2:3243966/150‑1 (MQ=255)
                                                                                                   gCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTCACGCGTCTTATCATGCCTACAGGTGTGTGCCGAACCGTAGGCCGGATAAGGCGTGCCCGCCGCATCCCGCCGGAGGCGCACACTGCATCAGGCGACGCt                                                    >  2:393011/1‑150 (MQ=11)
                                                                                                   gCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCt                                                    >  2:1345525/1‑150 (MQ=255)
                                                                                                    cGTTCACGCCGCCTCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCtt                                                   <  2:4710677/150‑1 (MQ=255)
                                                                                                    cGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCtt                                                   <  2:1737075/150‑1 (MQ=255)
                                                                                                      ttCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGt                                                 >  1:1869547/1‑150 (MQ=255)
                                                                                                           gccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGTCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGt                                            <  2:4966852/150‑1 (MQ=255)
                                                                                                            ccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTc                                           >  1:4611358/1‑150 (MQ=255)
                                                                                                                 tCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATc                                      >  2:5278503/1‑150 (MQ=21)
                                                                                                                    ggCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATg                                   <  2:3084406/150‑1 (MQ=255)
                                                                                                                    ggCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATg                                   >  1:4764991/1‑150 (MQ=255)
                                                                                                                      cAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGcc                                 <  1:1948697/150‑1 (MQ=255)
                                                                                                                       aGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCt                                >  2:4171635/1‑150 (MQ=12)
                                                                                                                           ggCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACaa                            <  1:3340841/150‑1 (MQ=255)
                                                                                                                           ggCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACaa                            <  1:1720211/150‑1 (MQ=255)
                                                                                                                            gcgcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAAc                                                                                                                    <  1:3887926/61‑1 (MQ=255)
                                                                                                                               cacaATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCt                        <  1:2113282/150‑1 (MQ=255)
                                                                                                                                acacggCCTGCTGCGACGCTTAACGCGTCTTATCATGCCTACCGGTTTGTGCCGACCCGTAGGCCGGCTAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCtg                       <  2:2147957/145‑1 (MQ=14)
                                                                                                                                  aaTGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCtgtg                     >  2:4427002/1‑150 (MQ=255)
                                                                                                                                  aaTGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCtgtg                     <  2:4082068/150‑1 (MQ=255)
                                                                                                                                       cTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACACGTCTGTGCCGaa                >  2:1793796/1‑150 (MQ=255)
                                                                                                                                        tGATGCGAAGCATAACGCGTCTTATAATGCATACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATACGGCAGTAGGCGCACAATGCCTGATGCGACGCTAGTCGCGTCTTAACATGCCTACAAGTATGTGCCGAAc               >  2:5058917/1‑150 (MQ=14)
                                                                                                                                              gACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgc                                                                                           >  1:1240290/1‑68 (MQ=255)
                                                                                                                                              gACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAgg         >  2:4340185/1‑150 (MQ=255)
                                                                                                                                              gACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCGTATCATGCCTACAAGTCTGTGCCGAACCGTAgg         >  2:4221042/1‑150 (MQ=255)
                                                                                                                                               aCGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGc        >  1:3221685/1‑150 (MQ=255)
                                                                                                                                                     aaCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATa  <  1:5965704/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATA  >  NC_000913/4614172‑4614470

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: