Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,985,112 C→G 15.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,1120CG15.0% 226.6 / 24.8 85intergenic (‑43/+27)araG/araFL‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base G (4/9);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGT  >  NC_000913/1985014‑1985210
                                                                                                  |                                                                                                  
ccGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATCACGTCTACGTGCCTCACTGGTCGTTCACACCCCCTCTGcatca                                                                                                   <  2:248525/100‑1 (MQ=255)
ccGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGcatca                                                                                                   <  2:1711111/100‑1 (MQ=255)
ccGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGcatca                                                                                                   <  2:2554122/100‑1 (MQ=255)
 cGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGcatcat                                                                                                  <  2:832661/100‑1 (MQ=255)
  gggAGACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGcatcat                                                                                                  <  2:2472462/99‑1 (MQ=255)
  gggAAACGTTTTACCGATGCCGCCAAAAGAGCGATACGGGATAGACGGTTGCATCACATCTCCGTGAATCACTGGTCGTCCACACCCCCTCTGCATCATg                                                                                                 >  2:1955132/1‑100 (MQ=255)
   ggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc                                                                                                <  1:2803262/100‑1 (MQ=255)
   ggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc                                                                                                >  1:1792988/1‑100 (MQ=255)
   ggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc                                                                                                <  1:1924536/100‑1 (MQ=255)
     acccGTTTTACCGATGGCGAGTAATGGGAGATACTGTGTAGACTGTTGCATAACGCCTCCCTGTATCACTGGTCGTTCGCACCTCCTCTGcgtca                                                                                                   <  2:1216185/92‑1 (MQ=255)
     aaaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCag                                                                                              <  2:2182961/100‑1 (MQ=255)
       aCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCagag                                                                                            >  1:343417/1‑100 (MQ=255)
        cGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAgg                                                                                           <  1:1696695/100‑1 (MQ=255)
        cGTTGTACCTATGCCCCGAAAGGAGAGATACGGGGTTGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTTACACCCCCTCTGCATGATGTAGAgg                                                                                           <  2:2115928/100‑1 (MQ=255)
         gtcttACCGATTCCGCGAAATGAGAGATACGGGGTACACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCTCACCCCCTCTGCATGATGCAGAggg                                                                                          <  2:368044/97‑1 (MQ=255)
         gTTTTACCGATGCCGCGAAATGAGCGATACGGGGTAGACTGTTGCATAACGTCTCCGTTAATCACTGGTCGTTCACGCCCCCTCTGCATGATGCAGAggg                                                                                          <  2:873441/100‑1 (MQ=255)
         gTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTGCGTGAATCACTGGTCGTTCACACCCCCTTTGCATCATGCAGCggg                                                                                          >  2:2104531/1‑100 (MQ=255)
         gTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAggg                                                                                          <  1:3042970/100‑1 (MQ=255)
         gTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAgg                                                                                           >  1:152346/1‑99 (MQ=255)
          ttttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgggg                                                                                         <  2:1745733/100‑1 (MQ=255)
           tttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAATGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggggg                                                                                        <  1:2652668/100‑2 (MQ=255)
           tttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggggg                                                                                        <  2:2328420/100‑2 (MQ=255)
           tttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggggg                                                                                        <  1:2263353/100‑2 (MQ=255)
             tACCGATGCTGCTAAATGAGAGATACGGGGTATACTGTTTCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggggg                                                                                        <  1:247928/98‑2 (MQ=255)
             tACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgggggtg                                                                                      <  2:1600288/100‑4 (MQ=255)
             tACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAAt                                                                                      >  1:899188/1‑100 (MQ=255)
                cGATGCCTCGAAATGAGATATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAAtt                                                                                     <  1:1643934/98‑1 (MQ=255)
                   tGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTcc                                                                                <  2:716010/100‑1 (MQ=255)
                    tccGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCg                                                                               <  2:1456681/99‑1 (MQ=255)
                    gCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCg                                                                               <  1:2697509/100‑1 (MQ=255)
                    gCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCg                                                                               <  2:1792967/100‑1 (MQ=255)
                    gCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCg                                                                               <  1:1930517/100‑1 (MQ=255)
                     ccGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCgg                                                                              >  2:2839033/1‑100 (MQ=255)
                     ccGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCga                                                                              >  1:1456695/1‑99 (MQ=255)
                      cgcgAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGc                                                                             >  1:1883325/1‑100 (MQ=255)
                          aaaTGAGTGATAAGGGGTTGACTGTTGCATAACGTCTCCGTTAATCACTGGTCGTTCACCCCCCCTCTGCATCATGCAGAGGGGAATTTTTCCgg                                                                              <  2:2254855/95‑1 (MQ=255)
                              gagagaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAttact                                                                     >  2:1285035/1‑100 (MQ=255)
                               agagaTACGGGGTAGACTGTTGCATAACGTCTACGTGAAGCACTGGTCGTTCACACCAACTCTGCATCATGCAGAGGGGAATTGTTCCGGCAAAtta                                                                       <  2:1437017/97‑1 (MQ=255)
                                    tACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGc                                                               >  1:1890038/1‑100 (MQ=255)
                                     aCGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGcc                                                              >  2:349070/1‑100 (MQ=255)
                                     aCGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCACTGCATCATGCAGAGGGGAAATTTTCCGGCAAATTACTTACCGcc                                                              >  1:767056/1‑100 (MQ=255)
                                      cGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCt                                                             >  1:2785075/1‑100 (MQ=255)
                                       ggggTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTa                                                            <  2:676252/100‑1 (MQ=255)
                                        gggTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTaa                                                           >  1:2722053/1‑100 (MQ=255)
                                         ggTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGAGAATTTTTCCAGCAAATTACTTACCGCCTaaa                                                          >  1:380480/1‑100 (MQ=255)
                                          gTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTTACACCTCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAAc                                                         <  1:402520/100‑1 (MQ=255)
                                            aGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCt                                                       >  2:1176712/1‑100 (MQ=255)
                                                tGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCttttt                                                   <  1:2069835/100‑1 (MQ=255)
                                                 gTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCtttttt                                                  >  1:2209915/1‑100 (MQ=255)
                                                  ttGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtc                                                 <  2:2686270/100‑1 (MQ=255)
                                                  ttGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtc                                                 <  2:1662349/100‑1 (MQ=255)
                                                   tGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCTGGCAAATTACTTACCGCCTAAACCTTTTTtct                                                <  2:701392/100‑1 (MQ=255)
                                                   tGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtct                                                >  2:3008049/1‑100 (MQ=255)
                                                    gCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCGCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTGTTTtctc                                               >  2:1621956/1‑100 (MQ=255)
                                                    gCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtctc                                               <  2:2629541/100‑1 (MQ=255)
                                                    gCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtctc                                               <  1:1097690/100‑1 (MQ=255)
                                                     cATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTcc                                              <  1:2692642/100‑1 (MQ=255)
                                                     cATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTcc                                              >  1:500549/1‑100 (MQ=255)
                                                      aTAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCa                                             >  2:1881678/1‑100 (MQ=255)
                                                       tAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCTGTTTTTTTACTTACCGTCTAAACCTTTTTTTTCCAg                                            >  2:2416020/1‑100 (MQ=255)
                                                       tAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAg                                            <  2:1883303/100‑1 (MQ=255)
                                                       tAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAg                                            <  2:2711589/100‑1 (MQ=255)
                                                          cGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCCAAACCTTTTTTCTCCAGttc                                         <  1:2522889/100‑1 (MQ=255)
                                                            tctcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttctt                                       >  2:1330135/1‑100 (MQ=255)
                                                             ctcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttcttc                                      <  2:1981000/100‑1 (MQ=255)
                                                              tcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTtcttct                                     <  1:612012/100‑1 (MQ=255)
                                                               ccGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTcttctt                                    <  2:2772046/100‑1 (MQ=255)
                                                               ccGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTcttctt                                    >  1:144350/1‑100 (MQ=255)
                                                               ccGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTcttctt                                    >  1:2509059/1‑100 (MQ=255)
                                                                   gAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaaa                                <  2:3018141/100‑1 (MQ=255)
                                                                   gAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaaa                                <  2:1843987/100‑1 (MQ=255)
                                                                     aTCACTGTTCGTTCACTCTCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAg                               <  2:861878/99‑1 (MQ=255)
                                                                     aTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTCCTTTAAAGt                              >  1:1657297/1‑100 (MQ=255)
                                                                        aCTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCATAAACCTTTTTTCTCCAGTTCTTCTTTAAAGtt                             >  1:1675350/1‑98 (MQ=255)
                                                                         cTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTc                          >  2:1377537/1‑100 (MQ=255)
                                                                          tGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCa                         <  2:2244752/100‑1 (MQ=255)
                                                                            gTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcg                       >  2:2122336/1‑100 (MQ=255)
                                                                              cGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcgcg                     <  1:1913120/100‑1 (MQ=255)
                                                                                    cacCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCa               >  2:19633/1‑100 (MQ=255)
                                                                                    cacCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCa               >  2:831833/1‑100 (MQ=255)
                                                                                    cacCCCCTCTGCATCATCCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATc                <  2:338934/99‑1 (MQ=255)
                                                                                    attccccTCTGCATGATACAGAGGGGAATTTTTCCCGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCa               >  2:1453597/4‑100 (MQ=255)
                                                                                      tccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAAGTACTTACCGAATAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCCTGATCAGt             >  2:1314642/2‑100 (MQ=255)
                                                                                      cccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGt             >  2:1249316/1‑100 (MQ=255)
                                                                                      cccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGt             <  2:783743/100‑1 (MQ=255)
                                                                                       ccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCCAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTa            >  1:1919772/1‑100 (MQ=255)
                                                                                        cccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTAc           >  2:2515649/1‑100 (MQ=255)
                                                                                        cccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTAc           >  1:2778097/1‑100 (MQ=255)
                                                                                         ccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTAcc          <  1:878355/100‑1 (MQ=255)
                                                                                            ctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACg       >  2:567563/1‑100 (MQ=255)
                                                                                            ctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACg       >  2:1631252/1‑100 (MQ=255)
                                                                                               catcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCg    >  2:1427510/1‑100 (MQ=255)
                                                                                                 tcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATAAGTACCACGTCGGt  >  1:1141094/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
CCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGT  >  NC_000913/1985014‑1985210

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