Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,083 G→C 17.1% 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,0830GC17.1% 536.9 / 98.6 225intergenic (‑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 (93/90);  new base C (17/21);  total (111/114)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-02

Only 100 of 229 total aligned reads displayed.

GTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC  >  NC_000913/4311987‑4312182
                                                                                                |                                                                                                   
gTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCg                                                                                                  >  1:3840294/1‑100 (MQ=255)
 ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                                 <  2:2821157/100‑1 (MQ=255)
 ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                                 <  1:1770646/100‑1 (MQ=255)
  tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGa                                                                                                <  1:369353/100‑1 (MQ=255)
      aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACATGTAACTCCGCGATGGCGGAATTa                                                                                            <  2:614101/100‑1 (MQ=255)
      aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTa                                                                                            <  2:691797/100‑1 (MQ=255)
       acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                           >  2:438534/1‑100 (MQ=255)
       acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                           >  2:1537818/1‑100 (MQ=255)
        caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                          <  2:2639001/100‑1 (MQ=38)
        caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                          >  2:1186257/1‑100 (MQ=255)
         aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCt                                                                                         <  2:3393598/100‑1 (MQ=38)
         aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCt                                                                                         >  1:2988125/1‑100 (MQ=255)
         aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCCATCGCGGAGTTACCt                                                                                         <  2:858605/100‑1 (MQ=255)
          ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                        <  1:2023857/100‑2 (MQ=38)
          ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                        <  2:3514365/100‑2 (MQ=38)
            aaGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCtttt                                                                                      <  2:1291299/100‑1 (MQ=255)
             aGCCCACGAGTGTGCCGCTGAAATATTTTATATTTTTTTTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtc                                                                                      <  1:3253344/99‑4 (MQ=25)
             aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttt                                                                                     <  1:1609820/100‑1 (MQ=255)
             aGCCCACGAGTGTGCCGCTGAAATATATCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttt                                                                                     >  2:195222/1‑100 (MQ=255)
              gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                                    >  1:1822185/1‑100 (MQ=255)
              gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                                    <  1:2701075/100‑1 (MQ=255)
               cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                                   >  1:1476740/1‑100 (MQ=255)
                 cACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                                 <  2:2135139/100‑1 (MQ=255)
                  atgaatgtgACGCTGAAATAATTCATATAGTTATTCATAATAATGTTCTCATCTTAAGGTAACACATAACTCCGCGATGTCAGAATTACCTTTTTc                                                                                    <  2:1592609/91‑1 (MQ=255)
                  aCGAGTGTGCCGCTTAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGTTGACAGGTAACTCCTCGATGGCGGAATTACCTTTTTCtt                                                                                  <  2:1833719/98‑1 (MQ=255)
                   cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtcccctta                                                                                <  2:655960/99‑10 (MQ=25)
                   cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                               <  2:2988098/100‑1 (MQ=255)
                       tgtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAg                                                                           <  2:1201871/100‑1 (MQ=255)
                        gtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCAGAATTACCTTTTTCTTATTTCAgc                                                                          >  1:2736291/1‑100 (MQ=255)
                          gCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgc                                                                        <  1:2813612/100‑1 (MQ=255)
                           ccGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcg                                                                       >  2:3855600/1‑100 (MQ=255)
                            cGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAATGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcgc                                                                      <  2:2588130/100‑1 (MQ=255)
                            cGCTGAAATATCTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcgc                                                                      >  2:3105918/1‑100 (MQ=255)
                                gAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCAACGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTgg                                                                  >  1:1220633/1‑100 (MQ=255)
                                 aaaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGt                                                                 <  2:1886211/100‑1 (MQ=255)
                                 aaaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTAAGCGCGCCTGGt                                                                 >  2:3209688/1‑100 (MQ=255)
                                   atatTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTga                                                               <  1:1914579/100‑1 (MQ=255)
                                     atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                             <  2:1306751/100‑1 (MQ=255)
                                     atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                             >  1:458883/1‑100 (MQ=255)
                                     atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAATACCTTTTTCATATTTCAGCGCGCCTGGTgaga                                                             >  1:2295040/1‑100 (MQ=255)
                                     atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCAATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                             >  1:25052/1‑100 (MQ=255)
                                      tttCAGATATTTATTCATAATAATGTTATCACGGTAAGGGGACAGGTAAATCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgag                                                              >  2:2271368/1‑98 (MQ=255)
                                        tCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaaga                                                          >  1:1029610/1‑100 (MQ=255)
                                         cAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTCTTCTTATTTAAGCGCGCCTGGTGAgaataa                                                         >  2:2323071/1‑100 (MQ=255)
                                             tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGATATca                                                     <  1:2595837/100‑1 (MQ=255)
                                             tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATca                                                     <  2:2029879/100‑1 (MQ=255)
                                             tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTAAGCGCGCCTGGTGAGAAGaaa                                                        >  2:2254585/1‑97 (MQ=255)
                                               ttttttTATTATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcat                                                    <  2:2056366/99‑2 (MQ=255)
                                                ttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAGCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaaca                                                  <  2:3510033/100‑1 (MQ=255)
                                                ttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaaca                                                  >  2:2954081/1‑100 (MQ=255)
                                                    tCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTg                                              <  1:3629556/100‑1 (MQ=255)
                                                    tCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCAGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAtt                                               >  2:2491867/1‑99 (MQ=255)
                                                     cATAATAATGTTCGCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                             >  2:2422207/1‑100 (MQ=255)
                                                      ataataATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcg                                            <  1:142928/100‑1 (MQ=255)
                                                       taataaTGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGAGCACCTGGTGAGAAGAAATCAACAATTgcgc                                           >  1:1190105/1‑100 (MQ=255)
                                                        aataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGcc                                          >  1:2009393/1‑100 (MQ=255)
                                                        aataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGcc                                          >  1:3444959/1‑100 (MQ=255)
                                                        aataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGcc                                          <  1:1707409/100‑1 (MQ=255)
                                                        aataatGTTCTCACGGTAAGAGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGcc                                          >  2:2388866/1‑100 (MQ=255)
                                                         ataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCg                                         >  2:1426945/1‑100 (MQ=255)
                                                          taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                                        <  1:2950795/100‑1 (MQ=255)
                                                          taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGATGAATTACTTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAt                                                >  2:1873506/1‑92 (MQ=255)
                                                            atGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtt                                      <  1:2555406/100‑1 (MQ=255)
                                                             tGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAACGCGCCTGTTGAGAAGAAATCAACAATTACGCCGGttt                                     >  1:776109/1‑100 (MQ=255)
                                                                tctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                                  >  1:902870/1‑100 (MQ=255)
                                                                tctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTGCAGCGCGCCTGGGGAGAAGAAATCAACAATTGCGCCGGtttttt                                  >  2:3125298/1‑100 (MQ=255)
                                                                  tcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                                 <  1:594888/99‑1 (MQ=255)
                                                                   cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCa                               <  2:3451209/100‑1 (MQ=255)
                                                                   cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAACGCGCCTGGTGAGAAGAAATCAACAATTACGCCGATTTTTTGCa                               >  1:670816/1‑100 (MQ=255)
                                                                    aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCat                              <  1:999884/100‑1 (MQ=255)
                                                                     cGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCAGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCata                             >  1:2073191/1‑100 (MQ=255)
                                                                      ggTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatat                            <  1:696609/100‑1 (MQ=255)
                                                                       gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatt                           >  1:2202868/1‑100 (MQ=255)
                                                                       gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatt                           <  1:1683340/100‑1 (MQ=255)
                                                                         aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                         <  2:3481205/100‑1 (MQ=255)
                                                                         aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                         <  1:3120451/100‑1 (MQ=255)
                                                                           ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata                       >  1:866945/1‑100 (MQ=255)
                                                                           ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata                       >  2:1106860/1‑100 (MQ=255)
                                                                           agaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGCGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata                       >  1:2871454/7‑100 (MQ=25)
                                                                             ggAAAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatt                     <  2:532439/100‑1 (MQ=255)
                                                                              aaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGTatttttatatt                     >  2:3020105/4‑99 (MQ=25)
                                                                                aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc                  >  1:523795/2‑100 (MQ=38)
                                                                                aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc                  >  2:2765472/2‑100 (MQ=38)
                                                                                  ggTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt                >  2:743751/1‑100 (MQ=38)
                                                                                  ggTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt                >  1:2132737/1‑100 (MQ=38)
                                                                                  ggTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt                >  1:188617/1‑100 (MQ=38)
                                                                                  ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt                >  1:3169948/1‑100 (MQ=255)
                                                                                  ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt                 >  1:1656589/1‑99 (MQ=255)
                                                                                  ggTAACTCCGCGATAGCAGCATTACCTTTTACTTATTTCGGCGCACATGGTGAGACGCAATCAACCATTTCACCGGTCTTTTGCATATTCATATTAGCtt                >  2:3543795/1‑100 (MQ=255)
                                                                                    tAACTCCGCGATGGTGGAAATACATTAATCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCAGGTTTTTTGCatattaatattaacgtat              >  2:1337129/1‑94 (MQ=255)
                                                                                    tAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGt              <  1:2462268/100‑1 (MQ=255)
                                                                                    tAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGt              >  2:3701946/1‑100 (MQ=255)
                                                                                       cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCa           >  2:220857/1‑100 (MQ=255)
                                                                                         ccGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt         >  1:3750659/1‑100 (MQ=255)
                                                                                         ccGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt         >  1:3020894/1‑100 (MQ=255)
                                                                                          cgcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATg        >  2:447544/1‑100 (MQ=255)
                                                                                          cgcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGAATATTAATATTAGCTTGTTCAAt         >  1:2342818/1‑99 (MQ=255)
                                                                                           gcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGt       <  2:3750620/100‑1 (MQ=255)
                                                                                             gATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTg     <  1:1202143/100‑1 (MQ=255)
                                                                                                ggCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGc  <  2:1706153/100‑1 (MQ=255)
                                                                                                |                                                                                                   
GTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC  >  NC_000913/4311987‑4312182

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