Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,094,415 T→G 5.1% intergenic (‑143/+308) rhaD ← / ← rhaA rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,094,4150TG 5.1% 1076.6 / 30.7 354intergenic (‑143/+308)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base T (143/188);  new base G (12/9);  total (157/197)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 359 total aligned reads displayed.

TTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCT  >  NC_000913/4094266‑4094564
                                                                                                                                                     |                                                                                                                                                     
tttGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCt                                                                                                                                                       <  2:2681465/150‑1 (MQ=255)
tttGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCt                                                                                                                                                       <  2:2610461/150‑1 (MQ=255)
     aTGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATaa                                                                                                                                                  <  2:846873/150‑1 (MQ=255)
     aTGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATaa                                                                                                                                                  <  2:1108641/150‑1 (MQ=255)
       gTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGTCGACCTGTAGGCCTGATAAgg                                                                                                                                                <  1:2451310/150‑1 (MQ=255)
       gTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAgg                                                                                                                                                <  2:2559069/150‑1 (MQ=255)
         tCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCa                                                                                                                                              <  2:2041865/150‑1 (MQ=255)
           tgtgTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTAATCAGGCCTACAACTCGTGCCGACCAGTAGGCCTGATAAGGCCCt                                                                                                                                            >  1:4122994/1‑150 (MQ=21)
               ttCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTtgtg                                                                                                                                        >  1:3289189/1‑150 (MQ=255)
                tCCTGTAAATGCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGc                                                                                                                                       <  2:2272032/150‑1 (MQ=255)
                 gcggTACATGCGGGGTTGTCGGATGCCCGACCGGTAGGCCGGATAAGGCTCTTGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTCTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                                                                                                                      <  2:4007183/147‑1 (MQ=11)
                 ccTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGcc                                                                                                                                      <  2:275616/150‑1 (MQ=255)
                 ccTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGcc                                                                                                                                      <  1:3953276/150‑1 (MQ=255)
                 ccTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCATGATAAGGCACTTGTGcc                                                                                                                                      <  1:884215/150‑1 (MQ=255)
                         ttCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCgg                                                                                                                              <  1:559289/150‑1 (MQ=255)
                         ttCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAACGCGCTAGCGCCGCATACGGCAATAAATGCCTGAAGCGACTATGTCGCGTCTTATCAGGCCTACAACAATTGCCGACCTGTAGGCCTGAGAGGGCACGTGTGCCGCGGCCgg                                                                                                                              >  2:4093098/1‑150 (MQ=17)
                           cGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGCGCCGCATCCGGca                                                                                                                            >  1:3993558/1‑150 (MQ=255)
                            ggTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATTCCAGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGAAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaa                                                                                                                           >  1:1039/1‑150 (MQ=255)
                              tgttgtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaatc                                                                                                                         <  2:2131329/150‑1 (MQ=255)
                              tggtgtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaatc                                                                                                                         <  1:281977/147‑1 (MQ=255)
                               gttgtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaatca                                                                                                                        >  1:2654467/1‑150 (MQ=255)
                                 tgtCGGGTGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaatcaat                                                                                                                      <  1:1696682/150‑1 (MQ=255)
                                 tgtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaatcaat                                                                                                                      <  1:646531/150‑1 (MQ=255)
                                 tgtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaatcaat                                                                                                                      <  1:97043/150‑1 (MQ=255)
                                 tgtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGcaatcaat                                                                                                                      <  1:1141271/150‑1 (MQ=255)
                                  gtCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATg                                                                                                                     <  2:1718106/150‑1 (MQ=255)
                                   tCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGc                                                                                                                    <  2:1855063/150‑1 (MQ=255)
                                       aTGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGa                                                                                                                <  1:140677/150‑1 (MQ=255)
                                           aCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCAACCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCCGGCCTACAACTATTGCCGACCAGTAGGCCGGAAAAGGCACTTGTGCCGCATCCGGCACTCAAGGCCTGCTGCg                                                                                                            >  1:3491835/1‑150 (MQ=11)
                                            cGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGa                                                                                                           <  1:7112/150‑1 (MQ=255)
                                            cGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCATTGCCTGATGCGa                                                                                                           >  1:1609354/1‑150 (MQ=255)
                                             gACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGAc                                                                                                          <  2:2107189/150‑1 (MQ=255)
                                              aCCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCAGATGCGACGCTCTCTATGATTATCATGACACCAACAAATGCCGACCCGTGGTCCGGATATGGCCCTAGTGCCGCAGCCGCCCCACACTGCCTGAGGCGACg                                                                                                         >  1:1962386/1‑150 (MQ=2)
                                               cccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGc                                                                                                        <  2:133809/150‑1 (MQ=255)
                                                  gTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                                                                                     <  2:236127/150‑1 (MQ=39)
                                                  gTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                                                                                     <  1:2288585/150‑1 (MQ=39)
                                                     ggCCGGATAAGGCTCTAGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGTCGCTGTcgc                                                                                                  <  2:3267123/150‑1 (MQ=255)
                                                      gCCGGATAGGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGTTGCGCCGCTGTAGCGTCTTATAAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAATGCCCTTGTGCCGCTGACGTCCAACACTGACTGATGCGACGCAGTcgct                                                                                                 >  1:3396192/1‑149 (MQ=11)
                                                      gCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg                                                                                                 >  2:4206779/1‑150 (MQ=21)
                                                       ccGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGt                                                                                                <  1:1147306/150‑1 (MQ=39)
                                                        cGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                                               >  1:794757/1‑150 (MQ=39)
                                                        cGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                                               >  2:1477136/1‑150 (MQ=39)
                                                        cGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                                               >  1:3759607/1‑150 (MQ=39)
                                                              aGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCa                                                                                         <  1:3604572/150‑1 (MQ=39)
                                                               ggCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAg                                                                                        <  2:309873/150‑1 (MQ=39)
                                                               ggCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAGTGCCGACCTGGAGGCCTGATAATGCACTGGGGCCGCCTCCGGCCATCTATGCCTGATGCGACGCTGGCGCGTcgtctccg                                                                                        >  2:3447469/1‑143 (MQ=14)
                                                                gcgcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAgg                                                                                       >  2:1316431/1‑150 (MQ=38)
                                                                gcgcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAGTCCCCACCTGTAGGCCTGATACGGCACTTGTGCCGCATCCGGCAATCAAAGCCTGATGCGACGCTGTCGCGTCTTATCAgg                                                                                       >  2:3216282/1‑150 (MQ=11)
                                                                 cgcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                                                      <  1:2366731/150‑1 (MQ=38)
                                                                 cgcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                                                      >  1:3267432/1‑150 (MQ=38)
                                                                 cgcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                                                      <  1:210153/150‑1 (MQ=38)
                                                                  gcTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCAACCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCATATCAGGcc                                                                                     >  1:3159955/1‑150 (MQ=255)
                                                                      gcgcCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACACCTATTGCCGACCTGTAGGCCTGATAAGGCCCTGGTGCCGCATCCGGCCATCAATGCCTGTTGCGACGCGGTCGCGTCTTCTCAGGCCTGCa                                                                                 >  2:535278/1‑150 (MQ=11)
                                                                         ccgcATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACt                                                                              <  1:3681044/150‑1 (MQ=32)
                                                                          cgcATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTg                                                                             >  2:3953621/1‑150 (MQ=32)
                                                                           gcATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGt                                                                            <  1:3101205/150‑1 (MQ=32)
                                                                              tCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGc                                                                         <  1:2752606/150‑1 (MQ=32)
                                                                                 ggCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGa                                                                      <  1:2872076/150‑1 (MQ=32)
                                                                                  gCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGcc                                                                                     >  2:639299/1‑134 (MQ=32)
                                                                                  gCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCACATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGAc                                                                     <  2:2943946/150‑1 (MQ=17)
                                                                                          tGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGc                                                             <  1:4035126/150‑1 (MQ=34)
                                                                                           gCCTGCTGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGcc                                                            <  2:763554/150‑1 (MQ=21)
                                                                                              tGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGa                                                         <  1:4201728/150‑1 (MQ=34)
                                                                                              tGATGCGACGCTGTCGCGTCTTATCAGGCCAACAACGATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGCCGCTGTCGCGTCTTATCAGGGCTACAACTGTTGCCGACCTGTAGGCCGGa                                                         >  1:3337784/1‑150 (MQ=255)
                                                                                               gATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGAt                                                        <  1:1014131/150‑1 (MQ=34)
                                                                                                ctGCGACGCTGTCGCGTCTTATCCGGCCTACAACTAGTGCCGCCCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATa                                                       <  2:476658/149‑1 (MQ=255)
                                                                                                 tGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATaa                                                      <  2:2490148/150‑1 (MQ=34)
                                                                                                    gACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGc                                                   <  2:675595/150‑1 (MQ=34)
                                                                                                     aCGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCAGTAGGCCGGATAAGGcc                                                  >  2:3803324/1‑149 (MQ=17)
                                                                                                          gTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTtgt                                             <  1:3688163/150‑1 (MQ=37)
                                                                                                            cgcgTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGc                                           <  1:2952480/150‑1 (MQ=37)
                                                                                                             gcgTCTTCTCAGGCCCACAACTACTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                          <  1:1326664/150‑1 (MQ=255)
                                                                                                             gcgTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                          <  1:3712067/150‑1 (MQ=37)
                                                                                                              cgTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCg                                         <  1:2611314/150‑1 (MQ=37)
                                                                                                                     tCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAAGGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCgg                                  <  1:2184074/150‑1 (MQ=37)
                                                                                                                        ggCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaa                               <  2:3324279/150‑1 (MQ=37)
                                                                                                                          ccTACAACTATCGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatc                             <  1:609525/150‑1 (MQ=37)
                                                                                                                             aCAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatcaat                          <  2:1128882/150‑1 (MQ=37)
                                                                                                                                  tATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTg                     >  2:161238/1‑150 (MQ=35)
                                                                                                                                  tATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTg                     >  1:4168053/1‑150 (MQ=35)
                                                                                                                                  tATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTg                     >  2:1841260/1‑150 (MQ=35)
                                                                                                                                  tATTGCCGACCTGTAGGCCTGATAAGGCACTAGGGCCGCATCCGGCATTCAATCCCAGATGCGACGCTGTAGCTAATTATCATTCCTACAAAGGATCCCCACCAGTACGCAGGATAATGCCCTAGAGACGAATCCCGCACTCATTGccag                     >  2:4051743/1‑148 (MQ=2)
                                                                                                                                     tGCCGACCTGTAGGCCTGATAAGGCTCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATg                  <  1:1182163/150‑1 (MQ=21)
                                                                                                                                     tGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATg                  <  2:3405276/150‑1 (MQ=34)
                                                                                                                                       ccGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCg                >  2:1925506/1‑150 (MQ=34)
                                                                                                                                        cGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGa               <  2:2001600/150‑1 (MQ=34)
                                                                                                                                        cGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGa               <  1:3840625/150‑1 (MQ=34)
                                                                                                                                         gACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGAc              >  2:369237/1‑150 (MQ=34)
                                                                                                                                         gACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGAc              >  1:1590379/1‑150 (MQ=34)
                                                                                                                                             tGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCGGTAGGCCGGATAAGGCACTGGCGCCGCATCCGGCACTCAATCCCAGATGCGACGCTg          >  2:1246665/1‑150 (MQ=11)
                                                                                                                                              gTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCCATCAATGCCTGATGCGACGCTGt         >  2:3653296/1‑150 (MQ=21)
                                                                                                                                              gTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt         >  2:3736912/1‑150 (MQ=34)
                                                                                                                                               tAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCACTCAATGCCTGATGCGACGCGGTc        >  1:3244234/1‑150 (MQ=21)
                                                                                                                                               tAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCACTCAATCCCAGCTGCGACGCtcta        >  1:94387/1‑147 (MQ=11)
                                                                                                                                               tAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTc        >  2:1577126/1‑150 (MQ=34)
                                                                                                                                               tAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTc        <  1:1755242/150‑1 (MQ=34)
                                                                                                                                               tAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTc        >  2:3362948/1‑150 (MQ=34)
                                                                                                                                                  gCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgcg     >  1:2930885/1‑150 (MQ=34)
                                                                                                                                                   ccTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGt    <  1:986374/150‑1 (MQ=34)
                                                                                                                                                    cTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc   >  1:3189295/1‑150 (MQ=33)
                                                                                                                                                     |                                                                                                                                                     
TTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCT  >  NC_000913/4094266‑4094564

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: