Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,486,771 T→C 18.0% intergenic (+19/‑20) crp → / → yhfK cAMP‑activated global transcription factor, mediator of catabolite repression/putative transporter, FUSC superfamily inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,486,7710TC18.0% 432.8 / 87.6 206intergenic (+19/‑20)crp/yhfKcAMP‑activated global transcription factor, mediator of catabolite repression/putative transporter, FUSC superfamily inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (81/86);  new base C (15/22);  total (96/110)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-02

Only 100 of 207 total aligned reads displayed.

GGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATG  >  NC_000913/3486672‑3486867
                                                                                                   |                                                                                                
ggACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGt                                                                                                  >  1:3552556/1‑100 (MQ=255)
  aCGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTa                                                                                                <  2:1010155/100‑1 (MQ=255)
    gCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTAcc                                                                                              >  1:2600770/1‑100 (MQ=255)
     cATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCt                                                                                             <  1:3841452/100‑1 (MQ=255)
      aTTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTg                                                                                            >  2:1082504/1‑100 (MQ=255)
      aTTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTg                                                                                            <  2:2370289/100‑1 (MQ=255)
        tCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGt                                                                                          <  2:14301/100‑1 (MQ=255)
          tGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAg                                                                                        <  2:3584939/100‑1 (MQ=255)
           gAAGCTGCTCGAAGATCAGAACCTGATCTCCGCACACGGTAATACCATCGTCGTTTACGGCACTCGGTAATTCCATCGGAGTGGCGCGCTACCTGGTa                                                                                         >  1:586658/1‑98 (MQ=255)
             aGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAgcgc                                                                                     <  2:3211834/100‑1 (MQ=255)
              gATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAgcgcg                                                                                    <  2:1409886/100‑1 (MQ=255)
               gtgtgcggaGATCAGGTTCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAgcgcgc                                                                                   >  2:3732019/9‑100 (MQ=255)
               aTGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAgcgcgc                                                                                   >  2:658239/1‑100 (MQ=255)
               aTGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAgcgcgc                                                                                   >  2:1710535/1‑100 (MQ=255)
               aTGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAacgcgc                                                                                   <  2:246516/100‑1 (MQ=255)
               aTGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAATGGCGCGCTACCAGGTAacgcgc                                                                                   <  2:3283270/100‑1 (MQ=255)
                tGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGcc                                                                                  <  2:4284230/100‑1 (MQ=255)
                tGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGcc                                                                                  <  2:4369824/100‑1 (MQ=255)
                tGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGcc                                                                                  <  2:2595258/100‑1 (MQ=255)
                tGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGcc                                                                                  <  2:3122739/100‑1 (MQ=255)
                 gcgggAAGATCAGAGCCTGATCGCCGCAGACGGTAAAACCATCGTCGGTTACGGCACTCGGTAATCCCGTCGGAGTGGCGCGCTACCAGGTAacgcggc                                                                                  <  1:1772785/96‑3 (MQ=255)
                 gCTGGAAGATCAGCACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACGTGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGcc                                                                                  <  1:830258/99‑1 (MQ=255)
                 gCTGGAAGATCAGAACCTGGTCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGCCa                                                                                 <  2:1164078/100‑1 (MQ=255)
                 gCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGCCa                                                                                 <  2:1998979/100‑1 (MQ=255)
                  cTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAt                                                                                >  2:2723419/1‑100 (MQ=255)
                  cTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGCCac                                                                                <  2:1896184/100‑2 (MQ=39)
                  cTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGCCac                                                                                <  2:1320352/100‑2 (MQ=39)
                  cTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAACGCGCCac                                                                                <  2:1769362/100‑2 (MQ=39)
                    ggAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAttt                                                                              >  2:2459329/1‑100 (MQ=255)
                    ggAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAttt                                                                              >  1:886871/1‑100 (MQ=255)
                    ggAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGATAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAttt                                                                              >  1:2618781/1‑100 (MQ=255)
                      aaGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGACACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTg                                                                            <  1:3379994/100‑1 (MQ=255)
                       aGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGt                                                                           >  2:3772816/1‑100 (MQ=255)
                        gATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGtt                                                                          <  2:2509414/100‑1 (MQ=255)
                        gATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGtt                                                                          >  1:411970/1‑100 (MQ=255)
                        gATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTAGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGtt                                                                          >  2:685718/1‑100 (MQ=255)
                         aTCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGttt                                                                         >  2:4076450/1‑100 (MQ=255)
                         aTCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGtt                                                                          <  2:3974243/99‑1 (MQ=255)
                            aGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTccc                                                                      <  1:2157819/100‑1 (MQ=255)
                            aGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTcc                                                                       <  2:271096/99‑1 (MQ=255)
                             gAACCTGATCTCCGCACACGGTAAAAACATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTCGCGCGCCATTTTGTTTctcc                                                                     >  1:2770386/1‑100 (MQ=255)
                                 cTGATCGCCGCGCACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGAt                                                                 <  1:854165/100‑1 (MQ=255)
                                  tGATCTCCGCCCACGGTAAAACCATCGTCGTTTCCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTcc                                                                       >  1:3264739/1‑93 (MQ=255)
                                    aTCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGAtgtg                                                              >  2:902732/1‑100 (MQ=255)
                                       tacGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGgcg                                                           <  2:2397074/98‑1 (MQ=255)
                                          gCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAg                                                        <  2:333794/100‑1 (MQ=255)
                                            acacGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGAc                                                      <  2:354253/100‑1 (MQ=255)
                                                ggTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGAt                                                  <  2:55086/100‑1 (MQ=255)
                                                ggTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGAt                                                  <  1:2568237/100‑1 (MQ=255)
                                                     aaCCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATc                                             <  1:1806836/100‑1 (MQ=255)
                                                       ccATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGCTTCCCCCGATGTGGCGCAGACTGATTTATCAc                                           >  1:235845/1‑100 (MQ=255)
                                                       ccATCGACGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTCGTCGCGCTCCCTTTCGTCTCCCCCGATGTCGCGCAGACTGATTTAt                                              >  1:2743618/1‑97 (MQ=255)
                                                          tcgtcgTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCAcccc                                        <  1:2830749/100‑1 (MQ=255)
                                                           cgtcgtTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCg                                       >  1:554714/1‑100 (MQ=255)
                                                           cgtcgtTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCg                                       <  1:797431/100‑1 (MQ=255)
                                                            gtggggtACGGCACTCGTTAAGCGCGTCGGAGTGGCTCGGGACGGAGCAGCGCGCCATTTTGTGTCGACCGATGTGGCGCAGACTTATTGATCAcccc                                        <  1:4214974/92‑1 (MQ=255)
                                                             ttgtTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGa                                      <  2:3685938/97‑1 (MQ=255)
                                                               gtTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGatat                                   >  1:1710796/1‑100 (MQ=255)
                                                                   aCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAAc                               >  2:3328762/1‑100 (MQ=255)
                                                                    cGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACt                              <  2:4239612/100‑1 (MQ=255)
                                                                     ggCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCaa                                <  1:3200788/97‑1 (MQ=255)
                                                                     ggCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTa                             >  2:860777/1‑100 (MQ=255)
                                                                     ggCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTa                             <  1:115884/100‑1 (MQ=255)
                                                                      gCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTAt                            <  1:2515831/100‑1 (MQ=255)
                                                                      gCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTAt                            <  2:3568112/100‑1 (MQ=255)
                                                                      gCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCACCGATGTGGCGCAGACTGATTTATCATCCCGATATCAACTAt                            >  1:3083233/1‑100 (MQ=255)
                                                                      gCAATCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTAt                            <  1:1976133/100‑1 (MQ=255)
                                                                       cACTCGTTCCTCCCGTCGGAGTGGCGCGTTACCTGATAGCGCGTCATTTTGTTTCCCCCGATGTGTCGCAGACTGATTTATCACCCCGATATCAACTa                             <  1:1040146/98‑1 (MQ=255)
                                                                       cACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGTGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATg                           >  1:286234/1‑100 (MQ=255)
                                                                         cTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCa                         >  2:2519096/1‑100 (MQ=255)
                                                                          tCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCAc                        <  2:247672/100‑1 (MQ=255)
                                                                          tCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCAc                        >  2:1541163/1‑100 (MQ=255)
                                                                             ttAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTc                     <  1:4009065/100‑1 (MQ=255)
                                                                              tAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCg                    >  2:829769/1‑100 (MQ=255)
                                                                              tAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCg                    <  1:3940847/100‑1 (MQ=255)
                                                                               aaTCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGa                   <  1:2188374/100‑1 (MQ=255)
                                                                                aTCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGAc                  >  2:360162/1‑100 (MQ=255)
                                                                                 tCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACa                 <  1:1343333/100‑1 (MQ=255)
                                                                                   ccGTCTGAGTGGCGCGGTACCTGCGAGCACGCCAGTGTGGTACGCCCGGTGTGGCGCAGACTGATTAATGACCCCGATATCAACTATGCACTTCGACaa                <  1:3730814/99‑1 (MQ=255)
                                                                                    cGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAAc              <  1:3669803/100‑1 (MQ=255)
                                                                                       caaaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCt           >  2:2791293/6‑100 (MQ=38)
                                                                                        ggAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTg          >  2:2260897/1‑100 (MQ=255)
                                                                                         aaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACTCCGATATCAACTATGCACTTCGACAAACGCTg          >  1:1881776/4‑99 (MQ=38)
                                                                                         aaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTgg         >  2:464863/4‑100 (MQ=38)
                                                                                         aaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTgg         >  2:913514/4‑100 (MQ=38)
                                                                                         aaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTgg         >  2:298131/4‑100 (MQ=38)
                                                                                         aaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTg          >  2:1512057/4‑99 (MQ=38)
                                                                                           gTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGtn       <  1:3346349/100‑2 (MQ=255)
                                                                                           gTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTg       >  1:1251971/1‑100 (MQ=255)
                                                                                           atGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGGCAAACGCTGGTg       >  2:3464691/2‑100 (MQ=38)
                                                                                            tGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGc      >  2:4162378/1‑100 (MQ=255)
                                                                                            tGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGc      <  2:1051173/100‑1 (MQ=255)
                                                                                             ggCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGca     >  1:1051169/1‑99 (MQ=255)
                                                                                             ggCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCt     >  2:3007508/1‑100 (MQ=255)
                                                                                             ggCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCt     >  2:119633/1‑100 (MQ=255)
                                                                                             ggCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCt     >  2:538052/1‑100 (MQ=255)
                                                                                             ggCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTCTCACCCCGATATCAACTATGCACTTCGCCAAACGCTGGTGCt     >  1:1530658/1‑100 (MQ=255)
                                                                                             ggCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCt     >  1:4387417/1‑100 (MQ=255)
                                                                                               cgcgTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCCCCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTAt   >  1:330343/1‑100 (MQ=255)
                                                                                                gcgTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATg  <  1:3757320/100‑1 (MQ=255)
                                                                                                   |                                                                                                
GGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATG  >  NC_000913/3486672‑3486867

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: