Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,839,302 T→C intergenic (‑38/‑222) ascG ← / → ascF asc operon transcriptional repressor; prpBC operon repressor/cellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,839,3020TC98.8% 309.1 / ‑4.4 84intergenic (‑38/‑222)ascG/ascFasc operon transcriptional repressor; prpBC operon repressor/cellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (33/50);  total (33/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

CGCGATCTTTAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGT  >  NC_000913/2839159‑2839448
                                                                                                                                               |                                                                                                                                                  
cgcgATCTTTAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGt                                                                                                                                             <  2:1045145/151‑1 (MQ=255)
  cgATCTTTAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGttt                                                                                                                                           >  2:872926/1‑151 (MQ=255)
         tAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGAtt                                                                                                                                    >  1:289641/1‑151 (MQ=255)
           gTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGt                                                                                                                                  <  1:130675/151‑1 (MQ=255)
              tcCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCt                                                                                                                               >  2:1113794/1‑151 (MQ=255)
                cTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTg                                                                                                                             >  2:893441/1‑151 (MQ=255)
                  ggCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTg                                                                                                                             <  1:889093/149‑1 (MQ=255)
                  ggCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATAGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGc                                                                                                                            >  1:1153507/1‑150 (MQ=255)
                    cTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGcac                                                                                                                          >  2:384009/1‑150 (MQ=255)
                     tGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGcaca                                                                                                                         >  2:642961/1‑150 (MQ=255)
                      gACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGcaca                                                                                                                         <  2:871758/149‑1 (MQ=255)
                        cGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATcg                                                                                                                     <  1:381948/151‑1 (MQ=255)
                        cGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATcg                                                                                                                     >  1:271993/1‑151 (MQ=255)
                           aGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATcggcg                                                                                                                  <  1:434229/151‑1 (MQ=255)
                           aGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATcggc                                                                                                                   <  1:1019533/150‑1 (MQ=255)
                             ccATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCggg                                                                                                                >  2:1223700/1‑151 (MQ=255)
                              cATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTGCACAATCGGCggg                                                                                                                >  2:554514/1‑149 (MQ=255)
                               aTTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACTCCGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCggg                                                                                                                >  1:983294/1‑149 (MQ=255)
                                ccaCCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACAGTTTGCCGATTGTCCTTGCACAATcggc                                                                                                                   >  1:148508/3‑145 (MQ=255)
                                 tACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGaaa                                                                                                             <  1:113654/150‑1 (MQ=255)
                                 tACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGaaa                                                                                                             <  2:657572/150‑1 (MQ=255)
                                      gAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATAtt                                                                                                       <  1:654238/151‑1 (MQ=255)
                                          gCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAgg                                                                                                   <  1:362856/151‑1 (MQ=255)
                                           cACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGt                                                                                                  <  2:486776/151‑1 (MQ=255)
                                            aCGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTg                                                                                                 <  2:775602/151‑1 (MQ=255)
                                             cgcgGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGa                                                                                                <  2:1079006/151‑1 (MQ=255)
                                              gcgGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTg                                                                                                 <  1:439962/149‑1 (MQ=255)
                                               cgGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGAcc                                                                                              <  2:1122041/151‑1 (MQ=255)
                                               cgGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGAcc                                                                                              <  2:672814/151‑1 (MQ=255)
                                                 ggAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTTAGGTGACCgg                                                                                            >  1:780125/1‑151 (MQ=255)
                                                  gAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGt                                                                                           <  2:1196318/151‑1 (MQ=255)
                                                  gAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGt                                                                                           <  2:83435/151‑1 (MQ=255)
                                                    aaCGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAATATTCAGGGGACCgg                                                                                            >  1:1130652/1‑148 (MQ=255)
                                                    aaCGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGttt                                                                                         >  2:725723/1‑151 (MQ=255)
                                                        gTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTcaca                                                                                     <  2:391471/151‑1 (MQ=255)
                                                          cGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACaa                                                                                    <  1:893396/150‑1 (MQ=255)
                                                            cTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGGCCGGGGTCACaaa                                                                                   >  1:1150426/1‑149 (MQ=255)
                                                              tttGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAatat                                                                                <  2:355554/150‑1 (MQ=255)
                                                              tttGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAata                                                                                 <  2:1016819/149‑1 (MQ=255)
                                                               ttGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAAtata                                                                               <  1:58349/150‑1 (MQ=255)
                                                                tGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATaa                                                                              <  1:1243278/150‑1 (MQ=255)
                                                                 gAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATaaaa                                                                            <  2:1146881/151‑1 (MQ=255)
                                                                   atccccGGCGCGCTTCGCCACTTCCAGCAGCGTCGTCATCATTTTCACCCCTTTACTCGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATaaaaaa                                                                          <  2:351129/149‑1 (MQ=255)
                                                                      cccGGCGCGGTTCGCCACTTCCAGCATCGTCGTCACCATTTTCATCCCTATACACGCAATCAACGCATTGTACCGCCCCGGTTGCCGAATGACCATGCACAAGCGGGGGGGAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGa                                                                       <  1:679637/151‑1 (MQ=255)
                                                                      cccGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGa                                                                       <  1:1067862/151‑1 (MQ=255)
                                                                        cGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGaa                                                                      >  2:174481/1‑150 (MQ=255)
                                                                        cGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAAc                                                                     <  1:708479/151‑1 (MQ=255)
                                                                        cGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGGTTCACCAAAATAAAAAATGaa                                                                      >  1:895148/1‑150 (MQ=255)
                                                                         ggCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACa                                                                    <  1:69091/151‑1 (MQ=255)
                                                                         ggCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACa                                                                    <  2:266963/151‑1 (MQ=255)
                                                                           cgcgcTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAt                                                                  <  2:812485/151‑1 (MQ=255)
                                                                            gcgcTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTTAGGTGGCCGGGTTCACAAATATAAAAAATGAACAAtt                                                                 >  2:925110/1‑151 (MQ=255)
                                                                            gcgcTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAtt                                                                 <  2:1262085/151‑1 (MQ=255)
                                                                             cgcTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTc                                                                <  1:857787/151‑1 (MQ=255)
                                                                              gcTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCa                                                               <  1:506049/151‑1 (MQ=255)
                                                                               cTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCa                                                               <  1:604873/150‑1 (MQ=255)
                                                                                ttCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAATAATTAGGTGAACGGGTTCACAAATATAAAAAATGAAAAAtt                                                                 >  1:486595/1‑147 (MQ=255)
                                                                                 tCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCAct                                                             <  1:414033/150‑1 (MQ=255)
                                                                                   gCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActc                                                            >  1:472043/1‑149 (MQ=255)
                                                                                     cACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGTGTCACAAATATAAAAAATGAACAATTCACTCTCtt                                                        >  2:644791/1‑151 (MQ=255)
                                                                                      aCTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCtt                                                        >  2:948719/1‑150 (MQ=255)
                                                                                      aCTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTg                                                       <  1:485266/151‑1 (MQ=255)
                                                                                      aCTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTTACTCTCtt                                                        >  1:614275/1‑150 (MQ=255)
                                                                                         tCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCt                                                     <  1:960633/150‑1 (MQ=255)
                                                                                         tCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCt                                                     <  1:156770/150‑1 (MQ=255)
                                                                                             gCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTAttt                                                <  1:199979/151‑1 (MQ=255)
                                                                                                        atcatTTTCATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTa                                     <  2:239891/151‑1 (MQ=255)
                                                                                                                cATCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGa                              <  2:833129/150‑1 (MQ=255)
                                                                                                                 aTCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGAGCGGCTTCACAAAAATAAAAAATGAACAATTCACTCTCTTGCTTATTTTGTGACAACTATTCATGAtt                            >  1:471531/1‑151 (MQ=255)
                                                                                                                 aTCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGAtt                            <  1:243824/151‑1 (MQ=255)
                                                                                                                 aTCCCTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCAGGTGACCGGGTTCACAAAAATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGAtt                            >  1:95945/1‑151 (MQ=255)
                                                                                                                    ccTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTtg                         <  1:1095560/151‑1 (MQ=255)
                                                                                                                    ccTTTACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGCTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTTGTGACAACTATTCATGAtttt                          >  2:154572/1‑150 (MQ=255)
                                                                                                                      tttACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTTAGGTGACCGGGTTCACAAATATAAAAAATGAAAAATTCACTCTCTTGCTTATTTTGTGACAACTATTCATGAtttttt                        >  1:1248481/1‑148 (MQ=255)
                                                                                                                        tACACGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAATATTTAGGGGACCGGTGACACAAATATAAAAAATGAACAATTCAATCTCTTGCTTATTTTGTGACAAATATTCATGATTTtgt                        >  1:489958/1‑148 (MQ=255)
                                                                                                                         acacGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGaaa                    <  1:136353/151‑1 (MQ=255)
                                                                                                                         acacGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTTTTGACAAATATTCATGATTTTGTTaa                     >  1:706303/1‑150 (MQ=255)
                                                                                                                          cacGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAAc                   <  1:1223636/151‑1 (MQ=255)
                                                                                                                           acGCAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAAcc                  >  2:999584/1‑151 (MQ=255)
                                                                                                                            cGCAATCAACGCAGTGTACCACACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATAGTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATCCGCTCGCTTGCTTATTTAGTGACAACTATTGATGATTTTGTGAAACCg                 <  1:302020/151‑1 (MQ=255)
                                                                                                                              cAATCAACGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCgg                >  1:474346/1‑150 (MQ=255)
                                                                                                                                   aaCGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCtt          <  1:495053/151‑1 (MQ=255)
                                                                                                                                   aaCGCAGTGTACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCtt          >  2:960490/1‑151 (MQ=255)
                                                                                                                                            tACCGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGt  <  1:725778/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CGCGATCTTTAGTCTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGT  >  NC_000913/2839159‑2839448

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