Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,172,353 A→T 16.1% L174Q (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,3530AT16.1% 279.5 / 43.1 114L174Q (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/10);  total (61/53)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 120 total aligned reads displayed.

TCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAA  >  NC_000913/2172255‑2172451
                                                                                                  |                                                                                                  
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)
       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)
         aaGAGTTTGCCAGTGCCAGTTTTTCTGAACTAAGGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCGGGCCattcag                                                                                          <  1:530744/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)
           gATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTc                                                                                        <  2:2200547/100‑1 (MQ=255)
            aTTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattcagtgc                                                                                       <  1:277114/100‑9 (MQ=255)
            aTTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCATCACACTCTTTGCTCCCAGAGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTcc                                                                                       >  1:1411347/1‑100 (MQ=255)
             ttttGCCAGTGCCAGTTTTTTTTTTTTTTTTTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCg                                                                                      >  1:486758/1‑100 (MQ=255)
              tttGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCgg                                                                                     >  2:2146841/1‑100 (MQ=255)
                 gccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCAc                                                                                  >  1:2266134/1‑100 (MQ=255)
                  ccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCAcc                                                                                 <  2:697840/100‑1 (MQ=255)
                  ccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCAcc                                                                                 >  2:626756/1‑100 (MQ=255)
                  ccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCAcc                                                                                 >  2:2298603/1‑100 (MQ=255)
                  ccagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCAcc                                                                                 <  2:2522135/100‑1 (MQ=255)
                    agtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaa                                                                               <  1:1943263/100‑1 (MQ=255)
                      tgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaata                                                                             <  1:2112464/100‑1 (MQ=255)
                       gccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaataa                                                                            <  2:1797280/100‑1 (MQ=255)
                         cagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaataata                                                                          <  2:38595/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                                                                   >  1:317613/1‑100 (MQ=255)
                                tCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGt                                                                   >  2:2445474/1‑100 (MQ=255)
                                tCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGt                                                                   >  2:2446500/1‑100 (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                                                           >  1:1832110/1‑100 (MQ=255)
                                        aaTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtc                                                           >  2:1546904/1‑100 (MQ=255)
                                         aTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtct                                                          <  1:370783/100‑1 (MQ=255)
                                           gTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTca                                                        <  1:417960/100‑1 (MQ=255)
                                            tCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTcac                                                       >  2:2394542/1‑100 (MQ=255)
                                             cGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTcaca                                                      <  2:815904/100‑1 (MQ=255)
                                              gATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                     >  2:381033/1‑100 (MQ=255)
                                              gATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                     <  2:1886954/100‑1 (MQ=255)
                                              gATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                     >  2:1711576/1‑100 (MQ=255)
                                                tCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAAcc                                                   <  2:1045035/100‑1 (MQ=255)
                                                   ccGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCGCCAATAATAATAACGTTTTTATTCTCACAACCTTg                                                >  1:1966518/1‑100 (MQ=255)
                                                    cGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTgc                                               >  1:1256887/1‑100 (MQ=255)
                                                    cGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTgc                                               <  1:534582/100‑1 (MQ=255)
                                                     gTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTgcg                                              >  1:1590096/1‑100 (MQ=255)
                                                     gTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTgcg                                              <  2:1095638/100‑1 (MQ=255)
                                                      tCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTgcgc                                             <  1:1271599/100‑1 (MQ=255)
                                                       cacaCTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCCCCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCt                                            >  1:1826215/1‑100 (MQ=255)
                                                       cacaCTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCt                                            <  1:439813/100‑1 (MQ=255)
                                                        acacTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTa                                           >  2:1905954/1‑100 (MQ=255)
                                                         cacTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTaa                                          >  1:1721645/1‑100 (MQ=255)
                                                           ctctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTaaa                                         <  1:879333/99‑1 (MQ=255)
                                                            tctTTGCTCCCAGCGCGACAGCGCACTGCATGGCCCGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAAt                                        <  2:1411355/99‑1 (MQ=255)
                                                             ctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGa                                      <  1:440837/100‑1 (MQ=255)
                                                                    cccAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCAt                               >  1:1430567/1‑100 (MQ=255)
                                                                    cccAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCAt                               >  2:626517/1‑100 (MQ=255)
                                                                     ccAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATg                              >  1:1795757/1‑100 (MQ=255)
                                                                      cAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGc                             <  2:2177332/100‑1 (MQ=255)
                                                                      cAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGc                             <  1:475261/100‑1 (MQ=255)
                                                                      cAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGCTCCGGCACCAATAATAATACCGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGc                             <  1:1412135/100‑1 (MQ=255)
                                                                        gcgcGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAg                           <  2:1003770/100‑1 (MQ=255)
                                                                          gcgACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGAc                         <  1:1590033/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)
                                                                                 cgcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTTTTTTTTTTTCCTTGCGCTAAATGAAAAGCATGCAGACCAACGGt                  >  1:2263281/1‑100 (MQ=255)
                                                                                 cgcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACgg                   >  2:2593934/1‑99 (MQ=255)
                                                                                   cACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTaa                <  1:559758/100‑1 (MQ=255)
                                                                                    aaccattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAAt               >  1:1800487/7‑100 (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            >  1:1014482/1‑100 (MQ=255)
                                                                                       cattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCgg            >  2:1827176/4‑100 (MQ=255)
                                                                                       cattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCgg            >  2:1252636/4‑100 (MQ=255)
                                                                                       cattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCgg            >  1:581650/4‑100 (MQ=255)
                                                                                        aaTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGc           >  2:2484334/1‑100 (MQ=255)
                                                                                          tGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCt          <  1:1997997/99‑1 (MQ=255)
                                                                                            gCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCaa       >  2:568666/1‑100 (MQ=255)
                                                                                            gCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCaa       >  1:1179133/1‑100 (MQ=255)
                                                                                            gCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCACTAAATAAAAAGAATGCAGACAAACAGTAATCATCTCaa       >  1:1451766/1‑100 (MQ=255)
                                                                                             ccAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAAt      >  2:1557990/1‑100 (MQ=255)
                                                                                               agcagGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaa    <  2:1042490/100‑1 (MQ=255)
                                                                                                gcagGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaaa   >  2:1956897/1‑100 (MQ=255)
                                                                                                 cagGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaaaa  >  2:2203795/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
TCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAA  >  NC_000913/2172255‑2172451

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: