Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,728,118 (C)8→9 coding (1254/1455 nt) xylB ← xylulokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,728,1101.C97.6% 498.1 / 4.7 124coding (1262/1455 nt)xylBxylulokinase
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (23/98);  total (25/99)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 227 total aligned reads displayed.

CGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAAC  >  NC_000913/3728010‑3728204
                                                                                                     |                                                                                              
 gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATc                                                                                               <  1:6328908/101‑1 (MQ=255)
  aCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcc                                                                                              >  2:6329446/1‑101 (MQ=255)
   cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccc                                                                                             <  1:6328822/101‑1 (MQ=255)
   cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccc                                                                                             >  1:6329447/1‑101 (MQ=255)
    tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccc                                                                                            <  1:6329065/101‑1 (MQ=255)
    tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccc                                                                                            <  1:6328970/101‑1 (MQ=255)
     gTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccc                                                                                           >  2:6329449/1‑101 (MQ=255)
      ttCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccc                                                                                          <  2:6328962/101‑1 (MQ=255)
       tCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccc                                                                                         <  1:6328612/101‑1 (MQ=255)
         tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       <  2:6329263/101‑2 (MQ=255)
         tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       <  2:6329024/101‑2 (MQ=255)
         tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329452/1‑100 (MQ=255)
         tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       <  1:6328980/101‑2 (MQ=255)
          aaCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccccg                                                                                      >  2:6329453/1‑99 (MQ=255)
           aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccccgt                                                                                     <  1:6328687/101‑4 (MQ=255)
           aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccccgt                                                                                     <  2:6328934/101‑4 (MQ=255)
           aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccccgt                                                                                     <  1:6328816/101‑4 (MQ=255)
           aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  2:6329457/1‑98 (MQ=255)
            cGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccccgtc                                                                                    >  2:6329461/1‑97 (MQ=255)
            cGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccccgta                                                                                    >  1:6329459/1‑97 (MQ=255)
              gTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACg                                                                                   >  2:6329464/1‑101 (MQ=255)
              gTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACg                                                                                   <  2:6328698/101‑1 (MQ=255)
              gTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACg                                                                                   <  1:6329051/101‑1 (MQ=255)
               tAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACgg                                                                                  <  1:6328788/101‑1 (MQ=255)
                aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccc                                                                                             >  1:6329467/1‑88 (MQ=255)
                 gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTa                                                                                <  1:6329108/101‑1 (MQ=255)
                  ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  2:6329469/1‑91 (MQ=255)
                   tGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAAt                                                                              <  2:6329130/101‑1 (MQ=255)
                   tGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAAt                                                                              <  1:6328893/101‑1 (MQ=255)
                   tGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329471/1‑90 (MQ=255)
                    gCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACg                                                                                   <  2:6328920/95‑1 (MQ=255)
                     cGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCg                                                                            <  1:6328986/101‑1 (MQ=255)
                       gCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAg                                                                          <  2:6328892/101‑1 (MQ=255)
                       gCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGCCCCCACATCCCCCCCCCGTACGGTaa                                                                               <  1:6328719/96‑1 (MQ=255)
                        caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329475/1‑85 (MQ=255)
                         aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  2:6329480/1‑84 (MQ=255)
                          acaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329482/1‑83 (MQ=255)
                            aaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  2:6329484/1‑81 (MQ=255)
                             aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGctgctg                                                                    <  1:6328926/101‑1 (MQ=255)
                             aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329486/1‑80 (MQ=255)
                             aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329488/1‑80 (MQ=255)
                              ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGa                                                                   <  2:6328789/101‑1 (MQ=255)
                              ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGa                                                                   <  1:6328873/101‑1 (MQ=255)
                              ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccc                                                                                         >  1:6329491/1‑78 (MQ=255)
                               tCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGAc                                                                  >  2:6329492/1‑101 (MQ=255)
                                   tGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCt                                                              <  2:6328991/101‑1 (MQ=255)
                                   tGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccc                                                                                        >  2:6329496/1‑74 (MQ=255)
                                        gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329500/1‑69 (MQ=255)
                                        gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  2:6329501/1‑69 (MQ=255)
                                         aTTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329503/1‑68 (MQ=255)
                                           ttCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329506/1‑66 (MQ=255)
                                           ttCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccc                                                                                        >  2:6329507/1‑66 (MQ=255)
                                            tctcTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGc                                                     <  2:6328900/101‑1 (MQ=255)
                                              tctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                   <  1:6329041/101‑1 (MQ=255)
                                                tGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  2:6329511/1‑61 (MQ=255)
                                                 ggATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCa                                                <  1:6329202/101‑1 (MQ=255)
                                                 ggATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329513/1‑60 (MQ=255)
                                                  gATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCAt                                               <  1:6329223/101‑1 (MQ=255)
                                                  gATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329516/1‑59 (MQ=255)
                                                   aTTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATc                                              >  1:6329518/1‑101 (MQ=255)
                                                    ttCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                             <  1:6328839/101‑1 (MQ=255)
                                                    ttCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329520/1‑57 (MQ=255)
                                                     tCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTg                                            <  1:6329139/101‑1 (MQ=255)
                                                     tCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  2:6329521/1‑56 (MQ=255)
                                                       gccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  2:6329525/1‑54 (MQ=255)
                                                        ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  2:6329527/1‑53 (MQ=255)
                                                         cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                        >  2:6329529/1‑101 (MQ=255)
                                                          gcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  2:6329532/1‑51 (MQ=255)
                                                           cgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGAc                                          <  2:6328829/97‑1 (MQ=255)
                                                            gATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGAc                                          <  2:6329225/96‑1 (MQ=255)
                                                            gATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAg                                     <  1:6328942/101‑1 (MQ=255)
                                                            gATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  2:6329534/1‑49 (MQ=255)
                                                             aTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329536/1‑48 (MQ=255)
                                                              tCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATc                                                                             <  2:6329140/59‑1 (MQ=255)
                                                               cTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTAc                                  <  2:6328495/101‑1 (MQ=255)
                                                               cTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329538/1‑46 (MQ=255)
                                                                tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  2:6329540/1‑45 (MQ=255)
                                                                 gcgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329542/1‑44 (MQ=255)
                                                                  cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCa                               <  2:6328971/101‑1 (MQ=255)
                                                                  cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  2:6329544/1‑43 (MQ=255)
                                                                   gcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccccc                                                                                       >  1:6329545/1‑42 (MQ=255)
                                                                      aGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAc                           <  1:6329288/101‑1 (MQ=255)
                                                                      aGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccc                                                                                        >  2:6329549/1‑39 (MQ=255)
                                                                         cTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAcgcg                        <  2:6329047/101‑1 (MQ=255)
                                                                           tgctgcGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGcc                                       >  2:6329554/1‑84 (MQ=255)
                                                                           tgctgcGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgcc                      <  2:6328874/101‑1 (MQ=255)
                                                                             ctgcGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgcccc                    <  2:6329158/101‑1 (MQ=255)
                                                                              tgcGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccg                   <  1:6329011/101‑1 (MQ=255)
                                                                                    cAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCa             <  2:6329198/101‑1 (MQ=255)
                                                                                    cAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCa             <  2:6329009/101‑1 (MQ=255)
                                                                                     aGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAg                                                                          >  2:6329557/1‑39 (MQ=38)
                                                                                     aGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCaa            <  1:6328976/101‑1 (MQ=255)
                                                                                       tGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATc          <  1:6328933/101‑1 (MQ=255)
                                                                                          tGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc       <  1:6329259/101‑1 (MQ=255)
                                                                                           ggCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGt                                    <  1:6328949/71‑1 (MQ=255)
                                                                                           ggCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACg      <  2:6329020/101‑1 (MQ=255)
                                                                                             ccgcacaTCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAcg                          <  2:6329121/76‑1 (MQ=255)
                                                                                             ccccACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                      <  2:6328739/67‑1 (MQ=255)
                                                                                              cccACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATccgcccgc                                                 >  1:6329559/1‑55 (MQ=255)
                                                                                               ccACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAt           <  2:6329216/92‑1 (MQ=255)
                                                                                                     |                                                                                              
CGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAAC  >  NC_000913/3728010‑3728204

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