Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,083 G→C 16.9% intergenic (‑41/+18) 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,0830GC16.9% 299.4 / 52.5 113intergenic (‑41/+18)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 G (46/47);  new base C (11/8);  total (57/56)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 113 total aligned reads displayed.

GTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGA  >  NC_000913/4311987‑4312180
                                                                                                |                                                                                                 
gTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCg                                                                                                >  1:1553402/1‑100 (MQ=255)
 ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGcag                                                                                               <  1:1643110/100‑3 (MQ=255)
 ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                               <  1:439424/100‑1 (MQ=255)
  tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGa                                                                                              >  1:1123423/1‑100 (MQ=255)
  tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGa                                                                                              >  2:1455846/1‑100 (MQ=255)
  tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGa                                                                                              >  1:257205/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)
         aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTATGGGGACAGGTAACTCCGCGATGGCGGAATTACCt                                                                                       >  2:1257814/1‑100 (MQ=255)
          ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                      <  2:559020/100‑2 (MQ=38)
          ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                      <  2:2118873/100‑2 (MQ=38)
          ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                      <  2:1931486/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)
                ccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTATCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCtt                                                                                >  2:1038165/1‑100 (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:1779495/100‑1 (MQ=255)
                                  aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                              <  1:1324563/100‑1 (MQ=255)
                                  aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                              <  2:522284/100‑1 (MQ=255)
                                       ttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                         >  2:1629002/1‑100 (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)
                                                       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)
                                                                tctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                                >  1:157621/1‑100 (MQ=255)
                                                                 ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                               >  1:1518187/1‑100 (MQ=255)
                                                                   cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCa                             <  2:1618755/100‑1 (MQ=255)
                                                                    aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCat                            >  2:144003/1‑100 (MQ=255)
                                                                    aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCat                            <  1:1275848/100‑1 (MQ=255)
                                                                    aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCat                            <  1:1260928/100‑1 (MQ=255)
                                                                      ggTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatat                          <  1:1702098/100‑1 (MQ=255)
                                                                      ggTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatat                          >  1:1180402/1‑100 (MQ=255)
                                                                      ggTAAGCGGACAGGTAACTCCGCGATTGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTCGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatat                          <  2:370070/100‑1 (MQ=255)
                                                                       gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatt                         >  1:1335763/1‑100 (MQ=255)
                                                                       gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatt                         <  1:2275896/100‑1 (MQ=255)
                                                                         aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                       <  2:1553391/100‑1 (MQ=255)
                                                                         aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                       <  2:1865493/100‑1 (MQ=255)
                                                                         aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                       <  1:1909660/100‑1 (MQ=255)
                                                                         aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                       <  1:109458/100‑1 (MQ=255)
                                                                          aagaaaaaGGTAATTCCACGATCGAGGAGATACCTTTTTCTTATTTCAGCGCGCCTGGGGAGAAGAATTTATGAATTTAGCAGGTTTTTTGCatatt                         >  2:468778/8‑97 (MQ=255)
                                                                           ggagACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatta                        <  2:2385505/94‑1 (MQ=255)
                                                                           agaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata                     >  1:2314436/7‑100 (MQ=25)
                                                                           agaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata                     >  1:953140/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:1078982/2‑100 (MQ=38)
                                                                                aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc                >  2:2102562/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:987805/1‑100 (MQ=255)
                                                                                  ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt              >  1:365150/1‑100 (MQ=255)
                                                                                   atAACTCCGCGATGGCGGAATTACCCTTTTCTTATTTCAACACGCCTAGTAAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATAGTAGc                <  2:2232952/96‑1 (MQ=255)
                                                                                    tAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGt            >  1:516534/1‑100 (MQ=255)
                                                                                      tcTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc          <  1:1603409/99‑1 (MQ=255)
                                                                                      attCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc          <  2:1614296/98‑1 (MQ=255)
                                                                                      aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc          <  1:914456/100‑1 (MQ=255)
                                                                                      aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc          <  1:1248779/100‑1 (MQ=255)
                                                                                      aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc          <  1:280163/100‑1 (MQ=255)
                                                                                       cTCCGCGATGGCGGAATTACCTTTTTTTTTTTTTTTCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc          >  2:1603397/1‑99 (MQ=255)
                                                                                       cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCa         >  1:1514619/1‑100 (MQ=255)
                                                                                         ccGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt       <  1:1038169/100‑1 (MQ=255)
                                                                                         ccGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt       >  2:689223/1‑100 (MQ=255)
                                                                                         ccGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt       >  2:1420958/1‑100 (MQ=255)
                                                                                         ccGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt       >  1:1251156/1‑100 (MQ=255)
                                                                                          cgcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCaa        <  1:817275/98‑1 (MQ=255)
                                                                                          cgcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATg      <  2:2196642/100‑1 (MQ=255)
                                                                                          cgcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATg      <  2:122375/100‑1 (MQ=255)
                                                                                          cgcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATg      <  1:25265/100‑1 (MQ=255)
                                                                                                |                                                                                                 
GTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGA  >  NC_000913/4311987‑4312180

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: