Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,172,350 A→T 16.6% L175Q (CTG→CAG)  gatD ← galactitol‑1‑phosphate dehydrogenase, Zn‑dependent and NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,172,3500AT16.6% 274.7 / 44.4 114L175Q (CTG→CAG) gatDgalactitol‑1‑phosphate dehydrogenase, Zn‑dependent and NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (53/41);  new base T (8/11);  total (61/53)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

Only 100 of 122 total aligned reads displayed.

TGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAA  >  NC_000913/2172251‑2172449
                                                                                                   |                                                                                                   
tGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCca                                                                                                     <  2:2339387/100‑1 (MQ=255)
  cATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCcagc                                                                                                   >  1:697837/1‑100 (MQ=255)
  cATCGCACCGAAAGATTATGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCcagc                                                                                                   <  1:1097241/100‑1 (MQ=255)
    tCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACTGCGCACTGAATGGCcagcag                                                                                                 <  1:1588009/100‑1 (MQ=255)
    tCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCcagcag                                                                                                 <  2:412599/100‑1 (MQ=255)
      gCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGc                                                                                               <  1:1737110/100‑1 (MQ=255)
        aacGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCa                                                                                             >  2:1746578/3‑100 (MQ=255)
         ccGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCaa                                                                                            >  2:2340762/1‑100 (MQ=255)
           gAAAGGTTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACTCGCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattc                                                                                          <  1:503163/100‑4 (MQ=255)
           gAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattc                                                                                          <  2:511312/100‑4 (MQ=255)
           gAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattc                                                                                          <  1:752883/100‑4 (MQ=255)
           gAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATg                                                                                          >  2:2011429/1‑100 (MQ=255)
           gAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAAGGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattc                                                                                          <  2:611390/100‑4 (MQ=255)
            aaaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattca                                                                                         <  2:1049259/100‑5 (MQ=255)
            aaaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattca                                                                                         <  2:2536988/100‑5 (MQ=255)
            aaaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATgg                                                                                         >  1:1376064/1‑100 (MQ=255)
             aaGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattcag                                                                                        <  1:1827121/100‑6 (MQ=255)
              aGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattcagt                                                                                       <  1:410866/100‑7 (MQ=255)
               gATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattcagtg                                                                                      <  2:811308/100‑8 (MQ=255)
               gATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTc                                                                                      >  2:1093532/1‑100 (MQ=255)
                aTTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCATCACACTCTTTGCTCCCAGAGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTcc                                                                                     >  1:1411347/1‑100 (MQ=255)
                 ttttGCCAGTGCCAGTTTTTTTTTTTTTTTTTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCg                                                                                    >  1:486758/1‑100 (MQ=255)
                 ttttGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCg                                                                                    <  2:541889/100‑1 (MQ=255)
                  tttGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCgg                                                                                   >  2:2146841/1‑100 (MQ=255)
                     gccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCAc                                                                                >  1:2266134/1‑100 (MQ=255)
                      ccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCAcc                                                                               >  2:2298603/1‑100 (MQ=255)
                      ccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCAcc                                                                               >  1:1522948/1‑100 (MQ=255)
                      ccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCAcc                                                                               <  2:2522135/100‑1 (MQ=255)
                      ccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCAcc                                                                               >  2:626756/1‑100 (MQ=255)
                        agtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaa                                                                             <  1:1943263/100‑1 (MQ=255)
                         gtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaat                                                                            >  2:1800301/1‑100 (MQ=255)
                          tgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaata                                                                           <  1:2112464/100‑1 (MQ=255)
                           gccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaataa                                                                          <  2:1797280/100‑1 (MQ=255)
                              agtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaataataa                                                                       <  1:54592/100‑1 (MQ=255)
                               gtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGTAGGCCAATGGTTCCGGCACCaataataat                                                                      >  1:1562211/1‑100 (MQ=255)
                                  tttCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAAc                                                                   >  1:2177318/1‑100 (MQ=255)
                                    tCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGt                                                                 >  2:2445474/1‑100 (MQ=255)
                                    tCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGt                                                                 >  2:2446500/1‑100 (MQ=255)
                                    tCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGt                                                                 >  1:317613/1‑100 (MQ=255)
                                     ccgAACTAATTTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACg                                                                  <  1:216297/96‑1 (MQ=255)
                                       gAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGtttt                                                              <  2:754746/100‑1 (MQ=255)
                                       gAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGtttt                                                              <  2:2266149/100‑1 (MQ=255)
                                          cTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCGGGCAACAATAATAATAACGttttt                                                             >  1:972797/1‑98 (MQ=255)
                                            actGTCGCTCGCCGTCACACTCTTAGAACCCACCTCAACAGCGCACAGAATAGACATCCGTCTAATGGTTCCGCCACCAATAATAATAACGTTTTTATtc                                                         <  1:1599205/98‑1 (MQ=255)
                                            aaTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtc                                                         >  2:1546904/1‑100 (MQ=255)
                                            aaTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtc                                                         >  1:1832110/1‑100 (MQ=255)
                                             aTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtct                                                        <  1:370783/100‑1 (MQ=255)
                                               gTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTca                                                      <  1:417960/100‑1 (MQ=255)
                                                tCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTcac                                                     <  2:2280603/100‑1 (MQ=255)
                                                 cGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTcaca                                                    <  2:815904/100‑1 (MQ=255)
                                                  gATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                   >  2:1711576/1‑100 (MQ=255)
                                                  gATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                   <  2:1886954/100‑1 (MQ=255)
                                                  gATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                   >  2:1534282/1‑100 (MQ=255)
                                                  gATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                   >  2:381033/1‑100 (MQ=255)
                                                    tCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAAcc                                                 <  2:1045035/100‑1 (MQ=255)
                                                        cGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTgc                                             >  1:1256887/1‑100 (MQ=255)
                                                        cGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTgc                                             <  1:534582/100‑1 (MQ=255)
                                                           cacaCTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCCCCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCt                                          >  1:1826215/1‑100 (MQ=255)
                                                             cacTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTaa                                        >  1:1721645/1‑100 (MQ=255)
                                                               ctctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTaaa                                       <  1:879333/99‑1 (MQ=255)
                                                               ctctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAAt                                      >  2:1155642/1‑100 (MQ=255)
                                                               ctctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAAt                                      <  2:661241/100‑1 (MQ=255)
                                                                 ctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGa                                    <  1:440837/100‑1 (MQ=255)
                                                                      cTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGc                               >  2:330079/1‑100 (MQ=255)
                                                                        cccAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCAt                             >  2:626517/1‑100 (MQ=255)
                                                                        cccAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCAt                             >  1:1430567/1‑100 (MQ=255)
                                                                        cccAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCAt                             >  1:1368324/1‑100 (MQ=255)
                                                                         ccAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATg                            >  1:1795757/1‑100 (MQ=255)
                                                                          cAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGc                           <  1:475261/100‑1 (MQ=255)
                                                                          cAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGc                           <  2:2177332/100‑1 (MQ=255)
                                                                          cAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGCTCCGGCACCAATAATAATACCGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGc                           <  1:1412135/100‑1 (MQ=255)
                                                                              gcgACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGAc                       <  1:1590033/100‑1 (MQ=255)
                                                                               cgACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGAcc                      <  2:2141547/100‑1 (MQ=255)
                                                                                    gcgcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACgg                 >  2:1157027/1‑100 (MQ=255)
                                                                                    ccggaaccattGGCCTGCTGGCCATTGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACgg                 >  2:2464978/11‑100 (MQ=38)
                                                                                     cggaaccattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGt                >  2:467971/10‑100 (MQ=255)
                                                                                     cgcACTGAATGGCCAGCATGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTACTCTCATACCCTTGCGCTAAATGAAAAGCATGCAGACCAACgg                 <  1:271031/99‑1 (MQ=255)
                                                                                     cgcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCACTAATAATAACGTTTTTATTCTCACAACCTTGCACTAAATGAAAAGCATGCAGACCAACgg                 >  1:555038/1‑99 (MQ=255)
                                                                                     cgcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACgg                 >  2:2593934/1‑99 (MQ=255)
                                                                                      gcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTa               <  1:57199/100‑1 (MQ=255)
                                                                                       cACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTaa              <  1:559758/100‑1 (MQ=255)
                                                                                        aaccattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACATTGCGCTAAATGAAAAGCATGCAGACCAACGGTAAt             >  1:2255521/7‑100 (MQ=255)
                                                                                         cTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATc            >  2:782760/1‑100 (MQ=255)
                                                                                         cTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATc            <  2:1218101/100‑1 (MQ=255)
                                                                                         cTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCg           <  1:2268214/100‑1 (MQ=255)
                                                                                         accattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATc            >  2:630458/6‑100 (MQ=255)
                                                                                          tGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCg           <  2:1376072/100‑1 (MQ=255)
                                                                                           gAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCgg          >  1:1077672/1‑100 (MQ=255)
                                                                                           gAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCgg          >  2:2450527/1‑100 (MQ=255)
                                                                                           gAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCgg          >  1:1014482/1‑100 (MQ=255)
                                                                                           cattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCgg          >  2:1827176/4‑100 (MQ=255)
                                                                                           cattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCgg          >  2:1252636/4‑100 (MQ=255)
                                                                                            aaTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGc         <  1:2588767/100‑1 (MQ=255)
                                                                                              tGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCt        <  1:1997997/99‑1 (MQ=255)
                                                                                               ggCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACGTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCa      >  1:1662200/1‑100 (MQ=255)
                                                                                                gCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCaa     >  1:412439/1‑100 (MQ=255)
                                                                                                gCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCaa     >  1:1179133/1‑100 (MQ=255)
                                                                                                 ccAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAAt    >  2:1557990/1‑100 (MQ=255)
                                                                                                  cagcagGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATa   <  2:1014485/100‑1 (MQ=255)
                                                                                                   agcagGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaa  <  2:1042490/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAA  >  NC_000913/2172251‑2172449

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