Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,728,125 T→G 7.7% Y416S (TAC→TCC)  xylB ← xylulokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,728,1250TG 7.7% 2369.1 / 41.5 874Y416S (TAC→TCC) xylBxylulokinase
Reads supporting (aligned to +/- strand):  ref base T (259/524);  new base G (84/7);  total (343/531)
Fisher's exact test for biased strand distribution p-value = 3.61e-29
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 890 total aligned reads displayed.

GATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727976‑3728274
                                                                                                                                                     |                                                                                                                                                     
   aaGCGGCATATCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAAt                                                                                                                                                    <  2:6392413/150‑1 (MQ=255)
      cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGa                                                                                                                                                 <  1:3447418/150‑1 (MQ=255)
           tAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATACAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgctgc                                                                                                                                            <  2:6734650/150‑1 (MQ=255)
             aCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGGCGGGCGTCGGGctgctg                                                                                                                                          >  2:3377318/1‑150 (MQ=255)
              cGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGa                                                                                                                                         <  1:8012865/150‑1 (MQ=255)
              cGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTACGGTAATCGCGctgctgg                                                                                                                                         >  2:8158755/1‑149 (MQ=255)
               gCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGAc                                                                                                                                        <  2:79199/150‑1 (MQ=255)
                  gcgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGGCCGGGGGTCGAGCTGCTGGCCGc                                                                                                                                     >  2:5430680/1‑150 (MQ=255)
                  gcgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGAGTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGc                                                                                                                                     <  2:8021258/150‑1 (MQ=255)
                   cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGTTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCt                                                                                                                                    <  1:2533207/150‑1 (MQ=255)
                   cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACACCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCt                                                                                                                                    <  2:5249491/150‑1 (MQ=255)
                    gcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAAGGAGCGAGTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTg                                                                                                                                   <  2:2207426/150‑1 (MQ=255)
                     cATCTTGTAGATGCGTCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGTTTTCTCTGGATTCGCCGCGTTCTGCGCCAGCCTTTCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGa                                                                                                                                  <  1:1116937/150‑1 (MQ=255)
                     cATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTGCGGTAATCGAGCTGCGGACCGCTGa                                                                                                                                  >  2:4092094/1‑150 (MQ=255)
                         tGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCGCCGTACGGTACTCGAGCTGGTGACCGCTGATATc                                                                                                                              >  2:1146042/1‑150 (MQ=255)
                          ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                             <  1:3111181/150‑1 (MQ=255)
                          ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                             <  2:2285439/150‑1 (MQ=255)
                            tAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGc                                                                                                                           <  1:7833038/150‑1 (MQ=255)
                                     cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTg                                                                                                                  <  2:3642217/150‑1 (MQ=255)
                                        ttCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACg                                                                                                               <  1:5894204/150‑1 (MQ=255)
                                         tCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                              <  2:6053007/150‑1 (MQ=255)
                                         tCTAACGGTAGTTGCGGCAAAAATTCAATGAGCGTTTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                              <  1:4027408/150‑1 (MQ=255)
                                           tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                            <  1:4397282/150‑1 (MQ=255)
                                            aaCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAg                                                                                                           <  1:5253450/150‑1 (MQ=255)
                                            aaCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAg                                                                                                           <  1:6938447/150‑1 (MQ=255)
                                               ggTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTGCGGGGCTCGTGCGGCTGAGCGCGGGGAGCCGGCCGCAGCTGGCGcccggtg                                                                                                        >  1:2011781/1‑145 (MQ=255)
                                                gTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACt                                                                                                       <  1:8156519/150‑1 (MQ=255)
                                                   gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCAc                                                                                                    <  2:6255941/150‑1 (MQ=255)
                                                   gTTGCGGCAACAATTCAATGAGCGAGTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCAc                                                                                                    <  2:8367619/150‑1 (MQ=255)
                                                    ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGGGAGCGGGCGGCTGGCCGCGGCTCGCCGCCGGCGTCTGGCGCCAGTGCGCACt                                                                                                   >  2:3948097/1‑150 (MQ=255)
                                                    ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGGCCTCGAGCGGCTGCCCGCGGAGCTCCGCCCGCAGCTGACGCCAGTCCTCAcg                                                                                                   >  2:2593892/1‑149 (MQ=255)
                                                     tGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTa                                                                                                  <  2:7698874/150‑1 (MQ=255)
                                                      gCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAc                                                                                                 <  1:226233/150‑1 (MQ=255)
                                                       cGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAcg                                                                                                <  2:1710585/150‑1 (MQ=255)
                                                       cGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAcg                                                                                                <  2:4514737/150‑1 (MQ=255)
                                                           aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGGAGTCGAGCTGCTGACGGCTGTGGTCCGCCGGGGTGTGACGCCAGTACTCGCGACGcgcc                                                                                            >  2:74856/1‑150 (MQ=255)
                                                           aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTGAGCGAGCTGCTGACCGCTGATAGCCGCCCGCATCTGACGCCAGTTCTCGCTACGcgcc                                                                                            >  2:7088480/1‑150 (MQ=255)
                                                           aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgcc                                                                                            <  2:2392151/150‑1 (MQ=255)
                                                               aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgc                                                                                        <  1:9100645/150‑1 (MQ=255)
                                                                 tCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgccc                                                                                      <  2:3298325/150‑1 (MQ=255)
                                                                  cAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCgccccgcccc                                                                                     <  2:5600734/150‑1 (MQ=255)
                                                                   aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGGACGGTAATCGAGCTGCTGGCCGCTGATATCCGCCAGCAGCTGACGCCAGTACTCACTGCGCGCCCCGccccc                                                                                    >  2:576586/1‑150 (MQ=255)
                                                                    aTGAGCGATTTCTGTGGATTCGCGGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCa                                                                                   <  2:5011497/150‑1 (MQ=255)
                                                                    aTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCa                                                                                   <  2:3038877/150‑1 (MQ=255)
                                                                      gAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAt                                                                                 <  1:8272154/150‑1 (MQ=255)
                                                                       aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGGCGGGACTCGCGCTGCGGGCCGCTGGGATCCGCGGGCCGCTGACGGCCGTCCGCGCGGGGCGCCCCGcgcccgctc                                                                                >  1:1616594/1‑146 (MQ=255)
                                                                       aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTAGGGTCATCGAGCTGCTGAGCGCTGATATCCGCCCGCAGCTGACGCCCGGCCTCACTCCGCGCCCCGCCCCCCATc                                                                                >  2:7831877/1‑150 (MQ=255)
                                                                       aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGTGCTGCGGTCCGCGGATATCCGCCAGCGTCTGACGGCGGGACGCTCGACGCGCCCCGcccccgcgc                                                                                >  1:952059/1‑146 (MQ=255)
                                                                        gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGGCCGCGGGCTGCTGACCGCTGCTGTCCGCCAGCCTCTGGCGCCAGTCCTCTCTACGCGCCCCGcccccgcgcc                                                                               >  2:906640/1‑145 (MQ=255)
                                                                        gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGGCGCCCGTACTCGCTACGCGCCCCGCCCCCTATCa                                                                               >  2:211563/1‑150 (MQ=255)
                                                                         cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCCGACGc                                                                                                              >  1:5570251/1‑118 (MQ=255)
                                                                          gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                                             <  1:8670427/150‑1 (MQ=255)
                                                                             ttCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgccc                                                                                      >  2:4771607/1‑138 (MQ=255)
                                                                                tctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCATCGAGCTGCTGACCGCTGCTATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAGGCCACGTAacg                                                                       >  2:6877650/1‑149 (MQ=255)
                                                                                     aTTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTg                                                                  <  1:4862959/150‑1 (MQ=255)
                                                                                       tCGCCGCGATCTGCGCCAGCCTTTCTGCGCCGAGTGCTGGCCCCACATCCCCCCCCGTACGGTTATCGAGCTGCTGACCGCTGAGATCCGCCAGCATCTTACGCGAGTACTGACTACGCGCCCCGCCCCCTATCAACGTAACACTTTGCg                                                                <  2:7352029/150‑1 (MQ=255)
                                                                                          ccgcGATCTGCGCCCGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                             <  1:8858016/150‑1 (MQ=255)
                                                                                          ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGCCGGGAATCGCGCTGGTGACCGCTGAGATCCGCCCGCCTCTGACGCCAGTCCTCACTACGCGCCCCGCCCCCCATCAGCGGACCACGTTGCGGtt                                                             >  2:5896520/1‑150 (MQ=255)
                                                                                             cgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTaa                                                          <  1:6780758/150‑1 (MQ=255)
                                                                                                tCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATAc                                                       <  2:8494254/150‑1 (MQ=255)
                                                                                                  tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg                                                     <  2:7224766/150‑1 (MQ=255)
                                                                                                  tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGCATCGAGCTGCTGACCGCTGGTATCCGCCAGCGTCTGACGCCAGTACTCGCTACGCGCCCCGCCCCCCATCAACGGAACACTTTGCGGGTTAATACCg                                                     >  2:2922750/1‑150 (MQ=255)
                                                                                                  tGCGCCAGCATTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg                                                     <  2:1559741/150‑1 (MQ=255)
                                                                                                   gcgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGc                                                    <  2:4893529/150‑1 (MQ=255)
                                                                                                    cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCa                                                   <  2:1088728/150‑1 (MQ=255)
                                                                                                     gcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGGTTACTACCGcgg                                                  >  2:34988/1‑148 (MQ=255)
                                                                                                      ccAGCCTTGCTGCGCCCAGTGCTGTCCCCACATCCCCCCCCGTGCGGTAATCGGGCTGCTGGCCGCGGATATCCGCCCGCGTCTGGCGCCCGTGCTCCCTGCGCGCCCCGCCCCCTATCAACGTGACACTTTGCGGTTTGATGCCGCAgg                                                 >  2:2751556/1‑150 (MQ=255)
                                                                                                      ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTAGTGGAGCTGGTGACCGCTGATATCCGCCAGCATCTGACGCCGGTACTCGCTGCGCGCCCCGCCCCCAAGCGCCGGAACACTTTGGGGTTGAATACCGCAgg                                                 >  1:7249118/1‑150 (MQ=255)
                                                                                                      ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCAGCTGGCGCCAGTACGCACGACGCGCCCCGCCCCCAAGCAACGGAACGCTGTGCGGGGTAATACCGCAgg                                                 >  2:4592820/1‑150 (MQ=255)
                                                                                                       cAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGc                                                <  1:8396844/150‑1 (MQ=255)
                                                                                                           cTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGc                                            <  1:6474912/150‑1 (MQ=255)
                                                                                                            ttGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCa                                           <  1:1031036/150‑1 (MQ=255)
                                                                                                            ttGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGAACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGTCGCCAGTACTCACTACGCGCCCCGCCCCCAAACCACGTAACACTTTGCGGTTTAATACCGCAGGCATGcc                                           <  1:2981104/150‑2 (MQ=255)
                                                                                                            ttGCTGAGCCCAGTTCTGGACCCACATCCCACCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATATGACGCCAGTACTCACTACGCGCCCCGCCCCCCATCAACGTAACACTTTGCGGTTTAATACCTCAGGCATGCa                                           <  2:4381721/150‑1 (MQ=255)
                                                                                                               ctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGa                                        <  1:5422620/150‑1 (MQ=255)
                                                                                                                  cgcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACAt                                     <  1:8558461/150‑1 (MQ=255)
                                                                                                                   gcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATc                                    <  1:2742178/150‑1 (MQ=255)
                                                                                                                   gcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATc                                    <  1:61243/150‑1 (MQ=255)
                                                                                                                    cccAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATGTCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGGACCCCTTTGCGGGTTAATACCGCAGGCATGCACGACATcc                                   >  1:5146935/1‑150 (MQ=255)
                                                                                                                      cAGTGCTGGCCCCACATCCCCCCCCGTACGGTACTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                                 <  2:8674833/150‑1 (MQ=255)
                                                                                                                      cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                                 <  1:7651610/150‑1 (MQ=255)
                                                                                                                      cAGTGCTGGCCCCACATCCCCCCCCGTACCGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCGTCTGACGCCCGTACTCACGACGCGCCCCGCCCCCCATCAACGTAACGCGTTGCGGTTTAATACCGCAGGCATGCCCGACCTCCAt                                 >  2:3570003/1‑150 (MQ=255)
                                                                                                                       aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATGGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTTGGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                                <  2:889866/150‑1 (MQ=255)
                                                                                                                           cTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCAt                            >  2:3745093/1‑150 (MQ=255)
                                                                                                                           cTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCAGACGCCAGTACGCACTACGCGCCCCGCCCCCAATCAACCTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCaa                            >  1:3314131/1‑149 (MQ=255)
                                                                                                                              gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg                         <  1:4198177/150‑1 (MQ=255)
                                                                                                                              gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTGTGCGGTTTAATACCGCAGGCATGCACGACAGCCATGCCATCTg                         >  2:248646/1‑150 (MQ=255)
                                                                                                                               ccccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCCATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCCCGACATCCATGCCATCTGc                        >  2:6580582/1‑150 (MQ=255)
                                                                                                                               ccccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCAGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGc                        <  1:1448836/150‑1 (MQ=255)
                                                                                                                                cccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGcc                       <  1:3843775/150‑1 (MQ=255)
                                                                                                                                 ccACATCCCCCCCCGTACGGTACTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa                      <  2:1886078/150‑1 (MQ=255)
                                                                                                                                  cacaTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCACCGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAg                     <  2:1382982/150‑1 (MQ=255)
                                                                                                                                   acaTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgc                    <  1:7077099/150‑1 (MQ=255)
                                                                                                                                     aTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCGGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCAACTGCCAgcgc                  >  1:738117/1‑150 (MQ=255)
                                                                                                                                        cccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATa               >  1:6355684/1‑150 (MQ=255)
                                                                                                                                         ccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATaa              <  1:4681586/150‑1 (MQ=255)
                                                                                                                                          cccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAc             >  2:4265155/1‑150 (MQ=255)
                                                                                                                                              cGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCAc         <  2:4176825/150‑1 (MQ=255)
                                                                                                                                                     tAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTcc  >  1:1859003/1‑150 (MQ=255)
                                                                                                                                                     tAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTcc  <  1:8946954/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727976‑3728274

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