Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,985,112 C→G 17.0% 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,1120CG17.0% 316.0 / 58.2 124intergenic (‑43/+27)araG/araFL‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (59/44);  new base G (9/12);  total (68/56)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.58e-01

Only 100 of 132 total aligned reads displayed.

GCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTA  >  NC_000913/1985013‑1985211
                                                                                                   |                                                                                                   
gCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGcatc                                                                                                     <  1:2332800/100‑1 (MQ=255)
 ccGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGcatca                                                                                                    >  1:3134508/1‑100 (MQ=255)
  cGGGAGACGCTTTACCGATGGCGCGAGATGAGAAATACGGGGTGGACTGTTGCAGAAGGTCTCCGTGAACATCAGGTCGTTCACACCCCCTCTGgatgat                                                                                                   >  1:1390426/1‑100 (MQ=255)
   gggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATg                                                                                                  >  2:304991/1‑100 (MQ=255)
    ggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc                                                                                                 <  2:1478093/100‑1 (MQ=255)
    ggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc                                                                                                 <  2:2863082/100‑1 (MQ=255)
    ggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTACGTGACTCACTGATCGATCACCCCCCCTCTGCATCATGc                                                                                                 >  1:867704/1‑100 (MQ=255)
       aaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga                                                                                              <  1:973810/100‑1 (MQ=255)
         cGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAgg                                                                                            <  2:189380/100‑1 (MQ=255)
         cGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAgg                                                                                            <  2:1951945/100‑1 (MQ=255)
          gTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggg                                                                                           <  2:1542183/100‑1 (MQ=255)
          gTTTTACCGATGCCGCGAAATGAGAAAAACGGGGTAGACTGGTGCATAACGTCTCCGTGAATCACTGGTCGTTCACCCCCCCTCTGCATGATGCAGAggg                                                                                           <  1:3862739/100‑1 (MQ=255)
           ttttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAgggg                                                                                          <  2:678948/100‑1 (MQ=255)
            tttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggggg                                                                                         <  2:3540798/100‑2 (MQ=255)
            tttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggggg                                                                                         <  2:468214/100‑2 (MQ=255)
            tttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggggg                                                                                         <  1:2442858/100‑2 (MQ=255)
             ttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGaa                                                                                        <  2:3662253/100‑1 (MQ=255)
               aCCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAAtt                                                                                      <  1:47734/100‑1 (MQ=255)
                 cGATGCCGCGAAATGAGAGATACGTGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgggtg                                                                                         <  1:419187/95‑3 (MQ=255)
                  gATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAAttttt                                                                                   >  2:690119/1‑100 (MQ=255)
                    tGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgggggtgtgaacg                                                                                 <  2:1835151/100‑10 (MQ=38)
                       cgcgAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGc                                                                              <  2:3241979/100‑1 (MQ=255)
                       cgcgAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGc                                                                              >  1:1952956/1‑100 (MQ=255)
                       cgcgAAATGAGAGATACGGGCTATACTGTTCCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgggggtgtgaa                                                                                   <  1:1773848/95‑8 (MQ=37)
                        gcgAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCa                                                                             >  2:2131114/1‑100 (MQ=255)
                        gcgAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCa                                                                             <  1:4224282/100‑1 (MQ=255)
                             aTGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAtta                                                                        >  2:2859780/1‑100 (MQ=255)
                             aTGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAtta                                                                        >  2:683823/1‑100 (MQ=255)
                               gagagaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAttact                                                                      <  1:3566313/100‑1 (MQ=255)
                               gagagaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTCCATCATGCAGAGGGGAATTTCTCCGGCAAAttact                                                                      >  2:811884/1‑100 (MQ=255)
                                agagaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAttactt                                                                     <  1:1798739/100‑1 (MQ=255)
                                agagaTACGGGGTAGACTGTTGAATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAttactt                                                                     >  1:2553610/1‑100 (MQ=255)
                                     tACGGGGTAGTCTATTGCATAACGTCTCCTTGATTCACTGGTCTTTTACACCCCCTCTGCATCATGCAGAGGTGTATTTTTACGGCAAATTACTAACCg                                                                 <  1:103659/99‑1 (MQ=255)
                                      aCGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGcc                                                               <  1:1047417/100‑1 (MQ=255)
                                          ggTAGACTGTTGCATAACGTCTCCGTGAATCACTGGGCGTTCACACCCCCTCTGCATCATGCAGAGGGGCATTTTTCCGGCAAATTACTTACCGCCaaa                                                            >  2:2022670/1‑99 (MQ=255)
                                           gTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAAc                                                          >  2:3469096/1‑100 (MQ=255)
                                                 tGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCttttt                                                    <  1:3439770/100‑1 (MQ=255)
                                                  gTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCtttttt                                                   >  1:2877079/1‑100 (MQ=255)
                                                   ttGCATCACGTCTCCGTGAATCACTGGTCGTTCACACCCCCACTGCATCATGCAGAGGGGACTTTTTCCGGCAAATTACTTACCGCCTAAACCttttt                                                    <  1:2136215/98‑1 (MQ=255)
                                                    tGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACAGCCTAAAACTATTTtc                                                  >  1:3658343/1‑99 (MQ=255)
                                                     gCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtctc                                                <  1:908413/100‑1 (MQ=255)
                                                      cATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTcc                                               <  1:1823983/100‑1 (MQ=255)
                                                      cATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTcc                                               <  1:2136335/100‑1 (MQ=255)
                                                       aTAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTGCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCa                                              >  1:1228777/1‑100 (MQ=255)
                                                        tAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAg                                             >  1:3224980/1‑100 (MQ=255)
                                                        tAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAg                                             >  2:2948975/1‑100 (MQ=255)
                                                        tAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAg                                             >  2:532651/1‑100 (MQ=255)
                                                         aaCGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAg                                             >  1:2138027/1‑99 (MQ=255)
                                                           cGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttc                                          <  2:111451/100‑1 (MQ=255)
                                                            gTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttct                                         <  1:2859764/100‑1 (MQ=255)
                                                               tcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTtcttct                                      >  1:23333/1‑100 (MQ=255)
                                                               tcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTtcttct                                      >  2:336548/1‑100 (MQ=255)
                                                               tcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTtcttct                                      >  1:1314711/1‑100 (MQ=255)
                                                               tcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTtcttct                                      >  2:1938747/1‑100 (MQ=255)
                                                               tcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTtcttct                                      <  2:895752/100‑1 (MQ=255)
                                                                ccGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTcttctt                                     >  2:23087/1‑100 (MQ=255)
                                                                ccGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTcttctt                                     <  1:394035/100‑1 (MQ=255)
                                                                   tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaa                                  <  1:2448665/100‑1 (MQ=255)
                                                                   tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGAAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaa                                  >  1:4320965/1‑100 (MQ=255)
                                                                     aaTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAg                                <  2:2777157/100‑1 (MQ=255)
                                                                     aaTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAg                                <  2:2413963/100‑1 (MQ=255)
                                                                      aTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCTTCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTagatt                               >  1:536039/1‑96 (MQ=255)
                                                                      aTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGt                               >  2:3759880/1‑100 (MQ=255)
                                                                      aTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGt                               <  1:2919543/100‑1 (MQ=255)
                                                                        cACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTg                             >  2:3978105/1‑100 (MQ=255)
                                                                           tGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCa                          >  2:3901718/1‑100 (MQ=255)
                                                                           tGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCa                          <  2:883452/100‑1 (MQ=255)
                                                                           tGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCa                          <  2:3228563/100‑1 (MQ=255)
                                                                           tGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCAAGTTCTTCTTTAAAGTTGtaa                          >  1:2439839/1‑98 (MQ=255)
                                                                            ggTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCa                          >  1:597084/1‑99 (MQ=255)
                                                                            ggTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAc                         <  1:3652630/100‑1 (MQ=255)
                                                                               cGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcgcg                      >  1:442869/1‑100 (MQ=255)
                                                                               cGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcgcg                      <  1:532649/100‑1 (MQ=255)
                                                                               cGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcgcg                      >  1:1971995/1‑100 (MQ=255)
                                                                                 ttCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTCACGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAATGTTGTAACGAGt                     >  1:3434340/1‑99 (MQ=255)
                                                                                 ttCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTg                    <  1:3112718/100‑1 (MQ=255)
                                                                                   cacaCCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGAt                  >  2:843959/1‑100 (MQ=255)
                                                                                   cacaCCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGAt                  >  1:1664333/1‑100 (MQ=255)
                                                                                   cacaCCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGAt                  >  2:1208851/1‑100 (MQ=255)
                                                                                    acacCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATc                 >  2:1611095/1‑100 (MQ=255)
                                                                                     cacCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCa                <  1:1616891/100‑1 (MQ=255)
                                                                                     attccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCa                >  2:832154/4‑100 (MQ=255)
                                                                                      ttccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTATTTAAAGCTGTCACGCGTGATCAg               >  2:4111830/3‑100 (MQ=255)
                                                                                       tccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGt              >  2:3723056/2‑100 (MQ=255)
                                                                                       cccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAg               <  1:3217100/99‑1 (MQ=255)
                                                                                        ccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTa             >  2:909312/1‑100 (MQ=255)
                                                                                        ccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTa             >  2:31686/1‑100 (MQ=255)
                                                                                        ccccTCTGCATGAAGCAGAGGGGAATTTTTCGGGCAAATTACGTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTa             >  2:2974596/1‑100 (MQ=255)
                                                                                        ccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTa             <  1:1051317/100‑1 (MQ=255)
                                                                                         cccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTAc            <  2:3240465/100‑1 (MQ=255)
                                                                                          ccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTAcc           >  2:1991858/1‑100 (MQ=255)
                                                                                          ccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTAcc           <  1:909308/100‑1 (MQ=255)
                                                                                          ccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTAcc           >  1:3996658/1‑100 (MQ=255)
                                                                                               gCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTc      >  2:4249605/1‑100 (MQ=255)
                                                                                               gCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTc      >  1:1496859/1‑100 (MQ=255)
                                                                                               gCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTc      >  1:3766306/1‑100 (MQ=255)
                                                                                                 atcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCgg    <  2:2416870/100‑1 (MQ=255)
                                                                                                 atcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCgg    >  1:2952832/1‑100 (MQ=255)
                                                                                                   catGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTa  <  1:2797948/100‑1 (MQ=255)
                                                                                                   catGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTa  <  1:3613543/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: