Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,074 C→T 16.0% intergenic (‑32/+27) alsB ← / ← alsR D‑allose ABC transporter periplasmic binding protein/d‑allose‑inducible als operon transcriptional repressor; autorepressor; repressor of rpiR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,312,0740CT16.0% 520.4 / 72.4 218intergenic (‑32/+27)alsB/alsRD‑allose ABC transporter periplasmic binding protein/d‑allose‑inducible als operon transcriptional repressor; autorepressor; repressor of rpiR
Reads supporting (aligned to +/- strand):  ref base C (90/92);  new base T (18/17);  total (109/109)
Fisher's exact test for biased strand distribution p-value = 8.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

Only 100 of 228 total aligned reads displayed.

GCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCA  >  NC_000913/4311975‑4312173
                                                                                                   |                                                                                                   
gcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAAc                                                                                                     >  1:3557967/1‑100 (MQ=255)
gcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAAc                                                                                                     >  2:2446610/1‑100 (MQ=255)
gcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAAc                                                                                                     <  2:3120340/100‑1 (MQ=255)
 caAAAGCGCTGGTTGACAACATTAATCCCATGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAAc                                                                                                     <  1:601674/99‑1 (MQ=255)
 caAAAGCGCTGGTTGACAACATTAAGCCCACGATTGTGTCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACTGTAAGGGGACAGGTAAc                                                                                                     <  2:3111192/99‑1 (MQ=255)
      tcgcTGTTTGACAACATTAAGCACACGAGTGTGCCGCTGAAATATTACAAATATTTATTCATAATTATGTTATCATGATAAGGGGTCAGGTAACt                                                                                                    >  2:286821/2‑95 (MQ=255)
       cgcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGa                                                                                              >  2:2086594/1‑100 (MQ=255)
            gTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCg                                                                                         >  1:3840294/1‑100 (MQ=255)
            gTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCg                                                                                         <  2:130267/100‑1 (MQ=255)
            gTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCg                                                                                         <  1:2156991/100‑1 (MQ=255)
             ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                        <  2:2821157/100‑1 (MQ=255)
             ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                        <  1:2446633/100‑1 (MQ=255)
                  aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTa                                                                                   <  2:691797/100‑1 (MQ=255)
                   acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                  >  2:438534/1‑100 (MQ=255)
                   acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                  >  2:1537818/1‑100 (MQ=255)
                    caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                 <  2:2639001/100‑1 (MQ=38)
                    caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                 >  2:1186257/1‑100 (MQ=255)
                     aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCCATCGCGGAGTTACCt                                                                                <  2:858605/100‑1 (MQ=255)
                      ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                               <  2:3512326/100‑2 (MQ=38)
                      ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                               <  1:1910791/100‑2 (MQ=38)
                       tAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgt                                                                              <  1:2074927/100‑3 (MQ=38)
                       tAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCtt                                                                               <  2:2259415/99‑1 (MQ=255)
                         aGCCCACGAGTGTGTCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtcc                                                                            <  1:918219/100‑5 (MQ=25)
                         aGCCCACGAGTGTGCCGCTGAAATATTTTATATTTTTTTTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtc                                                                             <  1:3253344/99‑4 (MQ=25)
                         aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtcc                                                                            <  1:679081/100‑5 (MQ=37)
                         aGCCCACGAGTGTGCCGCTGAAATATATCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttt                                                                            >  2:195222/1‑100 (MQ=255)
                          gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                           <  1:2701075/100‑1 (MQ=255)
                           cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                          >  1:1476740/1‑100 (MQ=255)
                           cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                          >  1:1636922/1‑100 (MQ=255)
                            ccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCtt                                                                         >  2:2152175/1‑100 (MQ=255)
                            ccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCtt                                                                         >  2:2215106/1‑100 (MQ=255)
                            ccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                          <  2:341254/99‑1 (MQ=255)
                             cACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                        >  2:552742/1‑100 (MQ=255)
                             cACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                        <  2:229976/100‑1 (MQ=255)
                              atgaatgtgACGCTGAAATAATTCATATAGTTATTCATAATAATGTTCTCATCTTAAGGTAACACATAACTCCGCGATGTCAGAATTACCTTTTTc                                                                           <  2:1592609/91‑1 (MQ=255)
                                 aGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTc                                                                    <  1:2266446/100‑1 (MQ=255)
                                   tgtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                      <  2:500433/96‑1 (MQ=255)
                                   tgtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAg                                                                  <  2:1201871/100‑1 (MQ=255)
                                       ccGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcg                                                              >  2:3855600/1‑100 (MQ=255)
                                       ccGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcg                                                              >  2:2358706/1‑100 (MQ=255)
                                        cGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAATGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcgc                                                             <  2:2588130/100‑1 (MQ=255)
                                         gCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGcc                                                            <  1:1885005/100‑1 (MQ=255)
                                            gAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTgg                                                         <  2:659730/100‑1 (MQ=255)
                                            gAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCAACGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTgg                                                         >  1:1220633/1‑100 (MQ=255)
                                             aaaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGt                                                        <  2:1886211/100‑1 (MQ=255)
                                                tatTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTga                                                      <  2:3141756/99‑1 (MQ=255)
                                                 atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                    >  1:458883/1‑100 (MQ=255)
                                                 atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                    >  2:641193/1‑100 (MQ=255)
                                                 atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCAATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                    >  1:25052/1‑100 (MQ=255)
                                                   ttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAGCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                  <  1:3720817/100‑1 (MQ=255)
                                                    taatatatTTATTCATAATAATGTTTTCACGGTAATGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaa                                                   <  2:3479340/94‑1 (MQ=255)
                                                    tCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaaga                                                 >  1:1029610/1‑100 (MQ=255)
                                                      aGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGaaa                                               <  1:1701325/100‑1 (MQ=255)
                                                        atatTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATc                                             <  2:3071165/100‑1 (MQ=255)
                                                         tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATca                                            >  1:1741845/1‑100 (MQ=255)
                                                         tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATca                                            <  2:2029879/100‑1 (MQ=255)
                                                           ttttttTATTATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcat                                           <  2:2056366/99‑2 (MQ=255)
                                                            ttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAGCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaaca                                         <  2:3510033/100‑1 (MQ=255)
                                                                tCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCAGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAtt                                      >  2:2491867/1‑99 (MQ=255)
                                                                   taataaTGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGAGCACCTGGTGAGAAGAAATCAACAATTgcgc                                  >  1:1190105/1‑100 (MQ=255)
                                                                    aataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGcc                                 >  1:3444959/1‑100 (MQ=255)
                                                                    aataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGcc                                 >  1:2009393/1‑100 (MQ=255)
                                                                     ataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCg                                <  2:2969194/100‑1 (MQ=255)
                                                                      taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                               <  2:1007317/100‑1 (MQ=255)
                                                                      taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                               <  1:2950795/100‑1 (MQ=255)
                                                                      taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGATGAATTACTTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAt                                       >  2:1873506/1‑92 (MQ=255)
                                                                         tGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAACGCGCCTGTTGAGAAGAAATCAACAATTACGCCGGttt                            >  1:776109/1‑100 (MQ=255)
                                                                           ttCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGttttt                          <  1:3048154/100‑1 (MQ=255)
                                                                             ctcACGGTAATGGTTTAGTTATCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGATAACATCTCAACCATTTCGCCTGttttttt                        <  2:3539910/100‑2 (MQ=255)
                                                                             ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                        <  1:835171/100‑1 (MQ=255)
                                                                              tcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                        <  1:594888/99‑1 (MQ=255)
                                                                              tcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGc                       >  1:1778458/1‑100 (MQ=255)
                                                                               cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCa                      >  2:3011046/1‑100 (MQ=255)
                                                                                aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCat                     <  1:2582191/100‑1 (MQ=255)
                                                                                aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTCTTTTGCat                     >  1:3635861/1‑100 (MQ=255)
                                                                                 cGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCata                    >  2:3618861/1‑100 (MQ=255)
                                                                                     aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTTGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                <  2:1220619/100‑1 (MQ=255)
                                                                                     aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                <  1:3120451/100‑1 (MQ=255)
                                                                                      aagaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaat               >  2:3749701/8‑100 (MQ=25)
                                                                                      aGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                >  1:2596970/1‑99 (MQ=255)
                                                                                       ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              <  1:2816179/100‑1 (MQ=255)
                                                                                       ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              >  2:1106860/1‑100 (MQ=255)
                                                                                       agaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              >  1:1928364/7‑100 (MQ=37)
                                                                                       agaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGCGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              >  1:2871454/7‑100 (MQ=25)
                                                                                           aaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatta           >  1:2338650/3‑99 (MQ=38)
                                                                                           aaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAg          >  1:1508693/3‑100 (MQ=38)
                                                                                            aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         >  1:523795/2‑100 (MQ=38)
                                                                                            aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         >  2:2765472/2‑100 (MQ=38)
                                                                                            aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTAAGAAGAAAACAAAAATTACGCCGGTTTTTTGCATATTAATATTAGc         >  1:2029670/2‑100 (MQ=25)
                                                                                            aaGGTAATTCCGACATCGCGGAGTAACATTTTTCTTATTTAAGCGCGCCTGGTGAGAAGACATTAACAATTGCGCCGGGTTTTTGCATATAAATATTAGc         >  2:3354849/2‑100 (MQ=255)
                                                                                             aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         >  2:1477593/1‑99 (MQ=38)
                                                                                             aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATATCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatagta           >  2:241232/1‑97 (MQ=38)
                                                                                             aGGTAACTCCGCGATTTCTTTATTATCTTTTTCTTAGTTCAGCGTGCCTGGGGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                <  1:867481/92‑1 (MQ=255)
                                                                                             aGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt        <  2:1711443/100‑1 (MQ=255)
                                                                                              ggTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt       >  2:743751/1‑100 (MQ=38)
                                                                                              ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt       >  1:3169948/1‑100 (MQ=255)
                                                                                              ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt        >  1:1656589/1‑99 (MQ=255)
                                                                                                tAACTCCGCGATGGTGGAAATACATTAATCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCAGGTTTTTTGCatattaatattaacgtat     >  2:1337129/1‑94 (MQ=255)
                                                                                                tAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGt     <  1:1426960/100‑1 (MQ=255)
                                                                                                  aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGtt    >  1:2736307/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
GCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCA  >  NC_000913/4311975‑4312173

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: