Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,074 C→T 17.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,0740CT17.8% 247.4 / 36.1 104intergenic (‑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 (43/42);  new base T (11/8);  total (54/50)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

Only 100 of 111 total aligned reads displayed.

GCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCA  >  NC_000913/4311975‑4312173
                                                                                                   |                                                                                                   
tctaaaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAAc                                                                                                     <  1:1640809/97‑1 (MQ=255)
gcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAAc                                                                                                     <  2:10461/100‑1 (MQ=255)
 caAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCCTATTAATGTTCTCACGGTAAGGGGACAGGTAAc                                                                                                     <  2:1563580/99‑1 (MQ=255)
 caAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACt                                                                                                    >  2:2285226/1‑100 (MQ=255)
   aaaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAAcaga                                                                                                  >  2:1640796/1‑97 (MQ=255)
       cgttGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATTTTCTCACGGTAAGGGGACAAGTAACTCcgcg                                                                                               <  2:2268363/96‑1 (MQ=255)
       cgcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCcgcg                                                                                               >  2:1643940/1‑99 (MQ=255)
          tGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCCAtcg                                                                                           <  2:2043632/100‑3 (MQ=255)
            gTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGc                                                                                          <  2:1706568/99‑1 (MQ=255)
            gTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCg                                                                                         <  1:905792/100‑1 (MQ=255)
            gTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCg                                                                                         >  1:1553402/1‑100 (MQ=255)
             ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCgg                                                                                        <  1:1433073/100‑1 (MQ=255)
             ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGcag                                                                                        <  1:1643110/100‑3 (MQ=255)
             ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                        <  1:439424/100‑1 (MQ=255)
              tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGa                                                                                       >  1:257205/1‑100 (MQ=255)
              tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGa                                                                                       >  1:1123423/1‑100 (MQ=255)
              tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGa                                                                                       >  2:1455846/1‑100 (MQ=255)
               gACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGaa                                                                                      >  1:1090134/1‑100 (MQ=255)
                  aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTACAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTa                                                                                   <  1:1293841/100‑1 (MQ=255)
                    caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                 <  1:1174408/100‑1 (MQ=38)
                    caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                 <  1:2197516/100‑1 (MQ=38)
                    caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                 >  1:2199780/1‑100 (MQ=255)
                      ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                               <  2:1931486/100‑2 (MQ=38)
                      ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                               <  2:2118873/100‑2 (MQ=38)
                      ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                               <  2:559020/100‑2 (MQ=38)
                      ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCtt                                                                               >  2:922470/1‑100 (MQ=255)
                       tAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgt                                                                              <  1:2095995/100‑3 (MQ=38)
                          gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                           >  2:2232659/1‑100 (MQ=255)
                            ccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                           <  2:2198333/98‑1 (MQ=255)
                            ccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCtt                                                                         >  1:2196662/1‑100 (MQ=255)
                            ccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTATCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCtt                                                                         >  2:1038165/1‑100 (MQ=255)
                             cACGAGTGTGGCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGTCATGTAATTCCGTGATGGCGGAATTACCttttac                                                                           <  2:239589/97‑3 (MQ=255)
                               ctaGTGTGCCGCTAAAATATTTCAGATATTTATTCATAATAAAGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                      >  1:1615086/3‑100 (MQ=255)
                                 aGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                      <  2:240125/98‑1 (MQ=255)
                                  gtgtGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCa                                                                   <  2:612072/100‑1 (MQ=255)
                                    gtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgc                                                                 >  2:1931234/1‑100 (MQ=255)
                                     tgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTAACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcg                                                                >  2:456632/1‑100 (MQ=255)
                                      gCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCAGCGATGGCGGAATTACCTTTTTCTTATTTCAgcg                                                                >  2:1460714/1‑99 (MQ=255)
                                      gCCGCTGAAATATTTCAGATATTTATTCATAATAATGATCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgc                                                               >  2:1275841/1‑100 (MQ=255)
                                       gcGCTGAAGTATTTCAGATATTTATTCATAATAATGTTAGTACGGTTAGGGGACAAGTAACTCCGCGTTGGCGGAGTTAAGTTATGATTATTTCAgcgcg                                                              >  2:564088/2‑100 (MQ=255)
                                       ccGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcg                                                              <  1:556143/100‑1 (MQ=255)
                                        cGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcgc                                                             >  1:1205421/1‑100 (MQ=255)
                                            gAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTgg                                                         >  1:2376182/1‑100 (MQ=255)
                                            gAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTgg                                                         >  1:1477255/1‑100 (MQ=255)
                                              aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTGCCTTTTTCTTATTTCAGCGCGCCTGGTg                                                       <  2:1720402/100‑1 (MQ=255)
                                              aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                       <  2:522284/100‑1 (MQ=255)
                                              aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                       <  2:1779495/100‑1 (MQ=255)
                                              aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                       <  1:1324563/100‑1 (MQ=255)
                                                   ttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                  >  2:1629002/1‑100 (MQ=255)
                                                      aGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGTTAACTCCGCGATGGCGGCATTACCTTTTTCTTATTTCCGCGCGCCTGGTGagaa                                                   >  2:1507109/1‑96 (MQ=255)
                                                       gATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAAt                                              <  2:1615285/100‑1 (MQ=255)
                                                        atatTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATc                                             <  1:43246/100‑1 (MQ=255)
                                                         tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCGGCGCGCCTGGTGAGAAGAAATca                                            <  1:79217/100‑1 (MQ=255)
                                                         tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATca                                            >  2:1598103/1‑100 (MQ=255)
                                                          atttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaa                                           <  2:754897/100‑1 (MQ=255)
                                                          atttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaa                                           <  1:1877431/100‑1 (MQ=255)
                                                                tCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTg                                     >  2:888643/1‑100 (MQ=255)
                                                                  ataataATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGTGCCTGGTGAGAAGAAATCAACAATTgcg                                   >  2:754235/1‑100 (MQ=255)
                                                                  ataataATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcg                                   >  1:1293594/1‑100 (MQ=255)
                                                                   taataaTGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                  >  1:1223055/1‑100 (MQ=255)
                                                                   taataaTGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTGCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                  >  2:1230614/1‑100 (MQ=255)
                                                                     ataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCg                                <  1:950224/100‑1 (MQ=255)
                                                                      taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                               >  2:432899/1‑100 (MQ=255)
                                                                          gTTCTCACGGTAAGGGGACAGGTAACGCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttt                           <  2:1865026/100‑1 (MQ=255)
                                                                           ttCTCACGGTAAGGGGACAGGTAACTACGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGttttt                          <  1:1660903/100‑1 (MQ=255)
                                                                            tctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                         <  2:1123418/100‑1 (MQ=255)
                                                                             ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                        >  1:1518187/1‑100 (MQ=255)
                                                                             ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                        <  1:305604/100‑1 (MQ=255)
                                                                               cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCa                      <  2:1618755/100‑1 (MQ=255)
                                                                                aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCat                     <  1:1260928/100‑1 (MQ=255)
                                                                                aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCat                     >  2:144003/1‑100 (MQ=255)
                                                                                  ggTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatat                   <  1:1702098/100‑1 (MQ=255)
                                                                                  ggTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatat                   >  1:1180402/1‑100 (MQ=255)
                                                                                  ggTAAGCGGACAGGTAACTCCGCGATTGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTCGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatat                   <  2:370070/100‑1 (MQ=255)
                                                                                   gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGTTTTTTTGCatatt                  <  1:1620353/100‑1 (MQ=255)
                                                                                   gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatt                  <  1:2275896/100‑1 (MQ=255)
                                                                                   gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatt                  >  1:1335763/1‑100 (MQ=255)
                                                                                     aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                >  2:563172/1‑100 (MQ=255)
                                                                                     aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                <  2:1553391/100‑1 (MQ=255)
                                                                                     aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                <  1:109458/100‑1 (MQ=255)
                                                                                       ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              <  2:1205414/100‑1 (MQ=255)
                                                                                       ggagACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatta                 <  2:2385505/94‑1 (MQ=255)
                                                                                       agaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              >  1:953140/7‑100 (MQ=25)
                                                                                       agaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              >  1:2314436/7‑100 (MQ=25)
                                                                                       agaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCGGGTGAGAAGAAGTCAACAATTGCGCCGGTTTTTTGCatattaata              >  2:2071530/7‑100 (MQ=25)
                                                                                        gggACAGGTCACTCCGCGATGGGGGAATTACCGTGTTCGTAGGTCAGCACGCCTGGGGAGAAGAATTCAACAATTGCGCCTGTTGTGGGCatattaa                >  2:195988/1‑97 (MQ=255)
                                                                                        gggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatat             >  1:2120718/1‑100 (MQ=255)
                                                                                         aaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatt            >  1:1102679/5‑100 (MQ=37)
                                                                                           aaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAg          >  1:411355/3‑100 (MQ=38)
                                                                                            aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         >  2:2102562/2‑100 (MQ=38)
                                                                                            aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         >  2:1078982/2‑100 (MQ=38)
                                                                                             aagTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt        >  1:1614308/3‑100 (MQ=38)
                                                                                             aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGACGAAATCAACAATTGCGCCGGTTTTTTTCATATTAATATTAg          >  2:1588103/1‑98 (MQ=37)
                                                                                             aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt        >  2:457081/1‑100 (MQ=38)
                                                                                              ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt       >  1:365150/1‑100 (MQ=255)
                                                                                              ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt       >  1:987805/1‑100 (MQ=255)
                                                                                                tAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGt     >  1:516534/1‑100 (MQ=255)
                                                                                                  tcTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc   <  1:1603409/99‑1 (MQ=255)
                                                                                                  aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc   <  1:1248779/100‑1 (MQ=255)
                                                                                                   cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCa  >  1:1514619/1‑100 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: