Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,985,112 C→G 16.7% intergenic (‑43/+27) araG ← / ← araF L‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,985,1120CG16.7% 364.4 / 66.2 138intergenic (‑43/+27)araG/araFL‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (54/60);  new base G (12/11);  total (67/71)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-01

Only 100 of 151 total aligned reads displayed.

GCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTA  >  NC_000913/1985013‑1985211
                                                                                                   |                                                                                                   
gCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGcatc                                                                                                     >  2:319609/1‑100 (MQ=255)
 ccGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGcatca                                                                                                    >  1:287716/1‑100 (MQ=255)
  cGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGcatgat                                                                                                   <  1:214385/100‑1 (MQ=255)
  cGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGCTCACACCCCCTCTGcatgat                                                                                                   <  1:210800/100‑1 (MQ=255)
   gggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGcatc                                                                                                     <  2:3707607/97‑1 (MQ=255)
   gggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATg                                                                                                  <  2:967053/100‑1 (MQ=255)
    ggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc                                                                                                 <  1:1372594/100‑1 (MQ=255)
    ggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc                                                                                                 <  1:3704340/100‑1 (MQ=255)
     gAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCa                                                                                                <  2:2484856/100‑1 (MQ=255)
       aaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga                                                                                              >  2:3569852/1‑100 (MQ=255)
         cGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAgg                                                                                            <  1:3204402/100‑1 (MQ=255)
         cGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAgg                                                                                            <  1:3738834/100‑1 (MQ=255)
          gTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggg                                                                                           <  1:1905461/100‑1 (MQ=255)
          gTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggg                                                                                           <  1:1936267/100‑1 (MQ=255)
          gTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAggg                                                                                           <  1:2275347/100‑1 (MQ=255)
           ttttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgggg                                                                                          <  1:259897/100‑1 (MQ=255)
           ttttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgggg                                                                                          <  1:2235333/100‑1 (MQ=255)
            tttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggggg                                                                                         <  2:526428/100‑2 (MQ=255)
              tACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgggggtn                                                                                       <  1:3874676/100‑4 (MQ=255)
              tACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAAt                                                                                       >  2:2637519/1‑100 (MQ=255)
                ccGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgggggtgtg                                                                                     <  2:94557/100‑6 (MQ=255)
                ccGATGCCGCGAAATGAGAGATACGAGATAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAAttt                                                                                     >  2:719480/1‑100 (MQ=255)
                  gATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAAttttt                                                                                   >  2:1367023/1‑100 (MQ=255)
                     gCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgggggtgtgaacga                                                                                <  1:3595768/100‑11 (MQ=38)
                     gCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCg                                                                                <  2:674862/100‑1 (MQ=255)
                      ccGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCgg                                                                               <  2:602615/100‑1 (MQ=255)
                      ccGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCgg                                                                               >  1:816094/1‑100 (MQ=255)
                      ccGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCgg                                                                               <  2:1886905/100‑1 (MQ=255)
                       cgcgAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGc                                                                              <  1:839300/100‑1 (MQ=255)
                        gcgAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCa                                                                             >  1:1481008/1‑100 (MQ=255)
                         cgAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACTCCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCaa                                                                            >  2:2141799/1‑100 (MQ=255)
                            aaTGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAtt                                                                         <  2:1199309/100‑1 (MQ=255)
                            aaTGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAtt                                                                         >  1:2783290/1‑100 (MQ=255)
                             aTGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAtta                                                                        >  1:3309631/1‑100 (MQ=255)
                               gagagaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGGTCACACCCCCTCTGCATAATGCAGAGGGGAATTTTTCCGGCAAAttact                                                                      >  2:304481/1‑100 (MQ=255)
                                agagaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAttactt                                                                     <  2:828772/100‑1 (MQ=255)
                                 gagaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAttactta                                                                    >  2:2011190/1‑100 (MQ=255)
                                    aTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTACATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCg                                                                 <  1:3935256/100‑1 (MQ=255)
                                         gggTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTaa                                                            <  1:110223/100‑1 (MQ=255)
                                          ggTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTaaa                                                           >  1:2517767/1‑100 (MQ=255)
                                            tAGACTGTTGCATAACGTCTCCGTGAATCAATGCTCGTTCACACCCCCTCTGCATCATCCAGAGGGGAATCTTTCCAGCAAATTACTTAACGCCTaaa                                                           >  2:2595514/1‑98 (MQ=255)
                                              gACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGCGTGGAATTTTTCCGACCAATTACTTACCGCCTAAACCt                                                        >  2:904874/1‑99 (MQ=255)
                                               aCTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCttt                                                      <  1:1007608/100‑1 (MQ=255)
                                                cTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCtttt                                                     <  2:205802/100‑1 (MQ=255)
                                                cTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCtttt                                                     >  1:894361/1‑100 (MQ=255)
                                                cTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACACCCTCTGCATCATGCCGAGGGGACTTTTTCCGGCAAATTACTTACCGCCTAAACCtttt                                                     >  1:3074481/1‑100 (MQ=255)
                                                  gTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCtttttt                                                   >  2:9423/1‑100 (MQ=255)
                                                   ttGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtc                                                  <  1:366606/100‑1 (MQ=255)
                                                     gCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtctc                                                >  2:13285/1‑100 (MQ=255)
                                                      cATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTcc                                               <  2:1458265/100‑1 (MQ=255)
                                                      cATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTcc                                               <  2:1297871/100‑1 (MQ=255)
                                                       aTAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCa                                              >  1:3836219/1‑100 (MQ=255)
                                                           cGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttc                                          >  1:830328/1‑100 (MQ=255)
                                                             tctcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttctt                                        <  1:2542635/100‑1 (MQ=255)
                                                              ctcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttcttc                                       <  2:3135445/100‑1 (MQ=255)
                                                                ccGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTcttctt                                     <  1:1971076/100‑1 (MQ=255)
                                                                 cGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCttt                                    <  2:1631858/100‑1 (MQ=255)
                                                                 cGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCttt                                    >  2:1181106/1‑100 (MQ=255)
                                                                 cGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTACGGCAAATTACTAACCGCCTAAACCTTTTTTCTCCAGTTcttctt                                     >  2:1420708/1‑99 (MQ=255)
                                                                   tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaa                                  <  2:3719051/100‑1 (MQ=255)
                                                                   ccgatctcTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaa                                  <  2:3832447/93‑1 (MQ=255)
                                                                    gAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaaa                                 >  1:2906144/1‑100 (MQ=255)
                                                                     atctaCTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAg                                <  2:621868/96‑1 (MQ=255)
                                                                     aaTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCTTTAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAg                                >  2:456246/1‑100 (MQ=255)
                                                                     aaTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAg                                <  2:208587/100‑1 (MQ=255)
                                                                        cACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTg                             <  1:24592/100‑1 (MQ=255)
                                                                        cACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTg                             <  2:1243905/100‑1 (MQ=255)
                                                                         aCTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGt                            <  2:274918/100‑1 (MQ=255)
                                                                          cTGGTCGGTCACCCCCACTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCGAAACCTTTTTTCTCCAGTTCGTCTTTAAAGTTGt                            <  2:1937338/99‑1 (MQ=255)
                                                                           tGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCa                          <  2:1461162/100‑1 (MQ=255)
                                                                           tGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCa                          <  2:3018180/100‑1 (MQ=255)
                                                                            ggTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAc                         >  1:2746755/1‑100 (MQ=255)
                                                                            ggTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAc                         <  1:264225/100‑1 (MQ=255)
                                                                            ggTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAc                         >  2:1720377/1‑100 (MQ=255)
                                                                            ggTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAc                         >  2:2357251/1‑100 (MQ=255)
                                                                             gTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcg                        >  1:1270417/1‑100 (MQ=255)
                                                                               ggaaaaattccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcgcg                      >  1:1344606/10‑100 (MQ=38)
                                                                               cGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcgcg                      <  1:1546344/100‑1 (MQ=255)
                                                                               cGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcgcg                      >  2:3363217/1‑100 (MQ=255)
                                                                                gTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGt                     <  1:625904/100‑1 (MQ=255)
                                                                                   cacaCCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGAt                  >  2:3249316/1‑100 (MQ=255)
                                                                                    aattccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATc                 >  1:2985297/5‑100 (MQ=255)
                                                                                     cacCCCCTCTGCATCATGCAGAGGGGAATTTTTTTTTTTTTTTACTTACAGCATAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTTTAACGTGTGAt                  >  2:160847/1‑98 (MQ=255)
                                                                                     cacCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCa                >  1:478002/1‑100 (MQ=255)
                                                                                     attccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATc                 >  2:3245449/4‑99 (MQ=255)
                                                                                      ttccccTCTGCATGATGCAGAGGGGAATTTTTTTTTCTTTTTTTTTATCGCCTAACCCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTAACGCGTGATc                 >  2:1140752/3‑98 (MQ=19)
                                                                                       cccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCTTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGt              <  2:1702700/100‑1 (MQ=255)
                                                                                        ccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGGTCAGTa             >  1:3304653/1‑100 (MQ=255)
                                                                                        ccccTCTGCATGACGCAGAGGGGAATTTTTCCGGCAAAATACTTACCGCCTAACCCTTTTTTCTCCAGTTCTTCTTTAAACTTGTCACGCGTGATCAGTa             >  2:2432145/1‑100 (MQ=255)
                                                                                         cccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTAc            >  1:3187951/1‑100 (MQ=255)
                                                                                          ccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAgatcg           >  1:1219453/1‑96 (MQ=255)
                                                                                          ccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTAcc           >  1:348/1‑100 (MQ=255)
                                                                                           ctctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACTTTTTTTCTCCAGTTCTTCTTTACAATTTCCACGAGTGATCAGTACCa          >  2:3695059/1‑100 (MQ=255)
                                                                                           ctctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCa          <  2:3558337/100‑1 (MQ=255)
                                                                                                 atcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCgg    <  2:268747/100‑1 (MQ=255)
                                                                                                 atcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCgg    >  1:3671932/1‑100 (MQ=255)
                                                                                                 atcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCgg    <  2:3387959/100‑1 (MQ=255)
                                                                                                 atcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCgg    <  2:280518/100‑1 (MQ=255)
                                                                                                   catGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTa  <  1:1978240/100‑1 (MQ=255)
                                                                                                   catGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTa  <  1:2357263/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTA  >  NC_000913/1985013‑1985211

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