Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,074 C→T 18.8% 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,0740CT18.8% 306.9 / 53.8 137intergenic (‑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 (58/53);  new base T (20/6);  total (78/59)
Fisher's exact test for biased strand distribution p-value = 2.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 138 total aligned reads displayed.

CAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCA  >  NC_000913/4311976‑4312173
                                                                                                  |                                                                                                   
caAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACt                                                                                                    <  2:1801047/100‑1 (MQ=255)
 aaaaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTc                                                                                                   <  2:2373260/100‑1 (MQ=255)
 aaaaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTCCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTc                                                                                                   >  1:2248682/1‑100 (MQ=255)
       gcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCCCGa                                                                                              >  1:506151/1‑99 (MQ=255)
         tGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATgg                                                                                           >  2:2334086/1‑100 (MQ=255)
          ggTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGc                                                                                          >  1:751619/1‑100 (MQ=255)
          ggTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGc                                                                                          >  2:2131558/1‑100 (MQ=255)
          ggTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGc                                                                                          >  1:13442/1‑100 (MQ=255)
           gTTGACAACATTAAGCCCACGAGTGTGCCTCTGAAATATTTCAGATATTTATTCATAATAATGTTCCCACGGTAAGGGGACAGGTAACTCCGCGATGGCg                                                                                         <  1:1748254/100‑1 (MQ=255)
           gTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCg                                                                                         <  2:793139/100‑1 (MQ=255)
           gTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCCATcgcg                                                                                         <  1:309067/100‑1 (MQ=255)
            ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                        >  2:2815235/1‑100 (MQ=255)
            ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                        >  2:1198850/1‑100 (MQ=255)
            ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                        >  1:1866617/1‑100 (MQ=255)
             tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGa                                                                                       <  2:1408873/100‑1 (MQ=255)
               acaacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAt                                                                                     >  1:256157/1‑100 (MQ=255)
               acaacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAt                                                                                     >  1:1580263/1‑100 (MQ=255)
               aaaacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTGACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAt                                                                                     <  1:456709/98‑1 (MQ=255)
                caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAt                                                                                     >  2:2411530/1‑99 (MQ=255)
                  acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                  >  1:820334/1‑100 (MQ=255)
                   caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                 <  2:1284366/100‑1 (MQ=39)
                      tAAGGCCACGGGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgt                                                                              <  1:956555/100‑3 (MQ=37)
                      tAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgt                                                                              <  1:1713262/100‑3 (MQ=38)
                      tAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgt                                                                              <  1:2591696/100‑3 (MQ=38)
                        aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtcc                                                                            <  1:1205906/100‑5 (MQ=37)
                         gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                           >  2:2776707/1‑100 (MQ=255)
                         gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                           >  1:1514307/1‑100 (MQ=255)
                         gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                           <  2:251508/100‑1 (MQ=255)
                          cccACGAGTGTGCCGCTGAAATCTTTCAGATATTTATTCATCATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgt                                                                              <  1:725437/96‑3 (MQ=37)
                          cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                          >  2:2298488/1‑100 (MQ=255)
                            cACGAGTGTGCCGCTGAAATATTTCAGATATATATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAATACCTTTTTCTTa                                                                        >  1:188971/1‑100 (MQ=255)
                             aCGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAt                                                                       <  2:2177638/100‑1 (MQ=255)
                              cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACTGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                      <  1:1556754/100‑1 (MQ=255)
                              cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                      <  2:1207689/100‑1 (MQ=255)
                                aGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTc                                                                    >  2:160418/1‑100 (MQ=255)
                                aGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTc                                                                    <  1:2729129/100‑1 (MQ=255)
                                 gtgtGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCa                                                                   <  1:722434/100‑1 (MQ=255)
                                      ccGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcg                                                              >  1:33940/1‑100 (MQ=255)
                                          tGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTg                                                          <  1:2054299/100‑1 (MQ=255)
                                            aaaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGt                                                        <  2:256159/100‑1 (MQ=255)
                                             aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                       <  2:1442576/100‑1 (MQ=255)
                                             aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                       >  2:1656779/1‑100 (MQ=255)
                                             aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                       <  2:282802/100‑1 (MQ=255)
                                                 tttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCTCGATGGCGGAATTACATTTTTCTTATTTCAGCGCGCCTGGTGagaa                                                   >  1:1278632/1‑100 (MQ=255)
                                                   tCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaaga                                                 <  2:307300/100‑1 (MQ=255)
                                                       atatTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATc                                             <  2:2567707/100‑1 (MQ=255)
                                                        tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgagatcggaaga                                            >  1:2207921/1‑92 (MQ=255)
                                                        tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATca                                            >  2:1688980/1‑100 (MQ=255)
                                                        tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATca                                            <  2:1050566/100‑1 (MQ=255)
                                                         cgatcttCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaa                                           <  2:1099380/95‑1 (MQ=255)
                                                          tttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGTCTGGTGAGAAGAAATcaac                                          <  1:368976/100‑1 (MQ=255)
                                                             attCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAt                                       >  2:2170876/1‑100 (MQ=255)
                                                              ttCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAtt                                      >  1:1399805/1‑100 (MQ=255)
                                                              ttCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTGATTTCAGCGCGCCAGGTGAGAAGAAATCAACAAtt                                      >  1:2525367/1‑100 (MQ=255)
                                                               tCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaaagatcg                                     >  1:1099378/1‑94 (MQ=255)
                                                               tCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTg                                     <  2:359609/100‑1 (MQ=255)
                                                               tCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTg                                     >  2:1736400/1‑100 (MQ=255)
                                                                cATAATAATGTTCTCACGGTAAGGGGACAGGTCACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTg                                     >  1:1783798/1‑99 (MQ=255)
                                                                cATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                    >  2:2110018/1‑100 (MQ=255)
                                                                cATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                    >  2:30420/1‑100 (MQ=255)
                                                                cATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                    >  2:569521/1‑100 (MQ=255)
                                                                  taataaTGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                  >  1:993093/1‑100 (MQ=255)
                                                                    ataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCg                                >  2:440684/1‑100 (MQ=255)
                                                                     taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                               >  2:2324503/1‑100 (MQ=255)
                                                                       atGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtt                             <  2:2429805/100‑1 (MQ=255)
                                                                       atGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtt                             <  2:2842252/100‑1 (MQ=255)
                                                                        tGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAATTCAGCGCGCCTGGTGCGAAGAAATCAACAATTGCGCCGGatt                            >  1:1485253/1‑100 (MQ=255)
                                                                           tctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                         <  2:1944440/100‑1 (MQ=255)
                                                                            ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                        >  2:599781/1‑100 (MQ=255)
                                                                            ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                        <  2:56389/100‑1 (MQ=255)
                                                                            ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                        <  1:1477942/100‑1 (MQ=255)
                                                                            ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                        <  2:1008311/100‑1 (MQ=255)
                                                                              cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgagatcggaag                      >  2:1751978/1‑90 (MQ=255)
                                                                              cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCa                      >  2:1899287/1‑100 (MQ=255)
                                                                              cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCa                      >  2:449925/1‑100 (MQ=255)
                                                                              cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCa                      >  2:2056207/1‑100 (MQ=255)
                                                                                cGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCata                    >  2:342393/1‑100 (MQ=255)
                                                                                  gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatt                  <  1:2012890/100‑1 (MQ=255)
                                                                                  aataagaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatt                  >  2:2088289/11‑100 (MQ=25)
                                                                                    aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                <  1:1459105/100‑1 (MQ=255)
                                                                                     aagaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaat               >  2:939860/8‑100 (MQ=25)
                                                                                     aGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaat               <  2:2149561/100‑1 (MQ=255)
                                                                                      agaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              >  1:923336/7‑100 (MQ=37)
                                                                                        ggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatt            <  2:37784/100‑1 (MQ=255)
                                                                                        aaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatt            >  2:151079/5‑100 (MQ=37)
                                                                                         gACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatact            >  2:1065562/1‑97 (MQ=255)
                                                                                          aaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAg          >  1:987338/3‑100 (MQ=38)
                                                                                          aaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAg          >  1:1373582/3‑100 (MQ=38)
                                                                                           aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         >  2:2893383/2‑100 (MQ=38)
                                                                                           aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         >  1:2410731/2‑100 (MQ=38)
                                                                                           aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         >  1:1874584/2‑100 (MQ=38)
                                                                                           aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAACAAATCAACAATTGCTCGGGTTTTTTGCATATTAATATTAg          >  1:1415496/2‑99 (MQ=37)
                                                                                            tggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt        <  1:1362680/99‑1 (MQ=255)
                                                                                            aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt        >  2:1197177/1‑100 (MQ=38)
                                                                                             ggTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt       >  1:988492/1‑100 (MQ=39)
                                                                                             ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt       >  1:1950158/1‑100 (MQ=255)
                                                                                             ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt       >  1:1926285/1‑100 (MQ=255)
                                                                                             ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCta       >  2:1362682/1‑99 (MQ=255)
                                                                                                 aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc   <  2:277782/100‑1 (MQ=255)
                                                                                                  cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCa  >  1:2763057/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCA  >  NC_000913/4311976‑4312173

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