Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,080 G→C 17.4% intergenic (‑38/+21) 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,0800GC17.4% 281.5 / 51.7 109intergenic (‑38/+21)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 (47/43);  new base C (10/9);  total (57/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

Only 100 of 115 total aligned reads displayed.

CGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATG  >  NC_000913/4311982‑4312176
                                                                                                  |                                                                                                
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:1553402/1‑100 (MQ=255)
     gTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCg                                                                                            <  1:905792/100‑1 (MQ=255)
      ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCgg                                                                                           <  1:1433073/100‑1 (MQ=255)
      ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGcag                                                                                           <  1:1643110/100‑3 (MQ=255)
      ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                           <  1:439424/100‑1 (MQ=255)
       tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGa                                                                                          >  2:1455846/1‑100 (MQ=255)
       tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGa                                                                                          >  1:257205/1‑100 (MQ=255)
       tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGa                                                                                          >  1:1123423/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: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)
                              tgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTAACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcg                                                                   >  2:456632/1‑100 (MQ=255)
                               gCCGCTGAAATATTTCAGATATTTATTCATAATAATGATCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgc                                                                  >  2:1275841/1‑100 (MQ=255)
                                ccGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcg                                                                 <  1:556143/100‑1 (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)
                                            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)
                                                         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:305604/100‑1 (MQ=255)
                                                                      ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                           >  1:1518187/1‑100 (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:1180402/1‑100 (MQ=255)
                                                                           ggTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatat                      <  1:1702098/100‑1 (MQ=255)
                                                                            gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGTTTTTTTGCatatt                     <  1:1620353/100‑1 (MQ=255)
                                                                            gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatt                     >  1:1335763/1‑100 (MQ=255)
                                                                            gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatt                     <  1:2275896/100‑1 (MQ=255)
                                                                              aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                   <  2:1865493/100‑1 (MQ=255)
                                                                              aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                   <  1:109458/100‑1 (MQ=255)
                                                                              aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                   >  2:563172/1‑100 (MQ=255)
                                                                              aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                   <  2:1553391/100‑1 (MQ=255)
                                                                              aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                   <  1:1909660/100‑1 (MQ=255)
                                                                               aagaaaaaGGTAATTCCACGATCGAGGAGATACCTTTTTCTTATTTCAGCGCGCCTGGGGAGAAGAATTTATGAATTTAGCAGGTTTTTTGCatatt                     >  2:468778/8‑97 (MQ=255)
                                                                                ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata                 <  2:1205414/100‑1 (MQ=255)
                                                                                ggagACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatta                    <  2:2385505/94‑1 (MQ=255)
                                                                                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)
                                                                                     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: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:280163/100‑1 (MQ=255)
                                                                                           aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc      <  1:1248779/100‑1 (MQ=255)
                                                                                            cTCCGCGATGGCGGAATTACCTTTTTTTTTTTTTTTCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc      >  2:1603397/1‑99 (MQ=255)
                                                                                            cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCa     >  1:1514619/1‑100 (MQ=255)
                                                                                              ccGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt   >  2:1420958/1‑100 (MQ=255)
                                                                                              ccGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt   >  2:689223/1‑100 (MQ=255)
                                                                                              ccGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt   >  1:1251156/1‑100 (MQ=255)
                                                                                              ccGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt   <  1:1038169/100‑1 (MQ=255)
                                                                                               cgcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCaa    <  1:817275/98‑1 (MQ=255)
                                                                                               cgcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATg  <  1:25265/100‑1 (MQ=255)
                                                                                                  |                                                                                                
CGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATG  >  NC_000913/4311982‑4312176

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: