Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,652,855:1 (G)6→7 100% coding (362/846 nt) sseA → 3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,652,8491.G86.5% 318.9 / 45.3 104coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (7/7);  new base G (41/49);  total (48/56)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 115 total aligned reads displayed.

AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652750‑2652943
                                                                                                    |                                                                                                 
aGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCt                                                                                                    >  1:1037864/1‑100 (MQ=255)
aGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCt                                                                                                    >  1:1307764/1‑100 (MQ=255)
     aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑‑ggg‑g‑g                                                                                           <  1:907223/100‑1 (MQ=255)
     aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑‑ggg‑g‑g                                                                                           >  2:1208850/1‑100 (MQ=255)
       gccccTGATTGTATATGACGAAGGTAATGTTTTCTCAGCCCAACGAGCATGGTTGATGCTGCGAATCTTTGGTGTAGAGAAAGTGTCGATTCtn                                                                                                   <  1:808376/91‑2 (MQ=255)
       gaaCCTGATTGTCTATGACGCAGGTAATCTTTTCTCCGCCCCACGAAAATTGTGGATGCTGCGCATCTTTGGTGTAGAGAAAGTTTCGTTTCT‑‑TGG‑G‑GGt                                                                                         <  1:965505/98‑1 (MQ=255)
       gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGGCCGACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑‑ggg‑g‑ggg                                                                                         <  1:2547479/100‑2 (MQ=255)
       gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑‑ggg‑g‑ggg                                                                                         >  1:490655/1‑99 (MQ=255)
        cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑‑ggg‑g‑gggt                                                                                        >  1:1130037/1‑98 (MQ=255)
        cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑‑GGG‑G‑GGTg                                                                                        >  1:1803881/1‑100 (MQ=255)
          ccTGATTGGCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑‑ggg‑g‑gggtgg                                                                                      <  1:2153182/100‑5 (MQ=255)
            atcttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGTTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGAc                                                                                     <  1:1616477/97‑1 (MQ=255)
             gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACt                                                                                    <  1:1407227/100‑1 (MQ=255)
              aTTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGTGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑‑GGG‑G‑GGTGGACTTg                                                                                  <  2:801109/100‑1 (MQ=255)
              aTTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑‑GGG‑G‑GGTGGACTTg                                                                                  <  2:482072/100‑1 (MQ=255)
                tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGAcagat                                                                                 >  2:1616481/1‑96 (MQ=255)
                tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGc                                                                                 <  1:2474862/100‑1 (MQ=255)
                 gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCa                                                                                <  2:2673460/100‑1 (MQ=255)
                  tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAg                                                                               >  1:100761/1‑100 (MQ=255)
                   cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAgg                                                                              <  1:33052/100‑1 (MQ=255)
                   cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAg                                                                               <  1:200800/99‑1 (MQ=255)
                   cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTAGATTCTG‑GGG‑G‑GGTGGACTTGCAgg                                                                              >  2:453316/1‑100 (MQ=255)
                   cTATAACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGATCACCTCTGGTGTAGAAAAAGTGTCGATTATG‑GGG‑G‑GGTTGCCTTGCAg                                                                               >  1:659239/1‑99 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTg                                                                           <  2:584141/100‑1 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACCTGCAGGCTg                                                                           >  2:2446127/1‑100 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTgg                                                                          <  2:1608469/100‑1 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTgg                                                                          >  1:1347533/1‑100 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTgg                                                                          >  2:1498355/1‑100 (MQ=255)
                         cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCa                                                                        <  2:1945965/100‑1 (MQ=255)
                         cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCa                                                                        <  2:490681/100‑1 (MQ=255)
                          gAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGAGGCTGCGCACCTTTGGTGTAGAGAAAGTGTGGATTCTGGGGG‑G‑GGTGGACTTGCAGGCTGGc                                                                         <  1:1136995/99‑1 (MQ=255)
                           aaGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGCAAGTGTCGATCCG‑‑GGG‑GAGGGAGACTTGCAGGCTGGCAgc                                                                      >  1:1859436/1‑100 (MQ=255)
                            aGGTAATCTTTTCTCAGCCCCACGGGCATGGTGGATGCTGCGCACCTTTGGTGTAGGCAAAGTGTCGAGTCTG‑GGG‑G‑GGTGGACTTTCAGGCTGGCAgag                                                                     <  2:1965130/100‑3 (MQ=255)
                            aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGATAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAggg                                                                     <  2:1967854/100‑3 (MQ=255)
                            aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAgcg                                                                     <  1:101011/100‑1 (MQ=255)
                             ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAgcgc                                                                    <  1:2766123/100‑1 (MQ=255)
                             ggTAATCTGTTCTCAGCCCCACGGTCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTTTG‑GGG‑G‑GGTGGACTTGCAGGCTGGAAgcgc                                                                    <  1:4171/100‑1 (MQ=255)
                              gTAATCTTTGCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGcgcg                                                                   <  1:2245733/100‑1 (MQ=255)
                               tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCga                                                                  >  2:2463447/1‑100 (MQ=255)
                                  tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCgatga                                                               >  1:799690/1‑100 (MQ=255)
                                  tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACGTGCAGGCTGGCAGCGCgatga                                                               <  2:2017941/100‑1 (MQ=255)
                                    ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                             <  1:2813750/100‑1 (MQ=255)
                                     tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTATAGAAAGGGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATct                                                            <  2:229152/100‑1 (MQ=255)
                                     tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATct                                                            <  2:1130039/100‑1 (MQ=255)
                                       tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                          >  2:2492301/1‑100 (MQ=255)
                                          cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                        <  2:1684947/100‑1 (MQ=255)
                                           aGCCCCACTAGCATGGTGGATGCTGCGCACCTGTGGTGTAGAGAAAGTGGCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                        <  1:2106948/98‑1 (MQ=255)
                                           aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  1:2594504/100‑1 (MQ=255)
                                              cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   >  2:1664208/1‑100 (MQ=255)
                                              cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   >  1:1443721/1‑100 (MQ=255)
                                              cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   <  1:2631246/99‑1 (MQ=255)
                                                cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                 >  1:628876/1‑100 (MQ=255)
                                                cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                 >  2:728199/1‑100 (MQ=255)
                                                cACGAGCATGGTGGATGCTGCGCACCTCTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGCGGACTTGCAGGCCGGCAGCGCGATGCCCCGCGGTTag                                                    >  1:2311803/1‑97 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:1920097/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  1:956855/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:1608179/100‑1 (MQ=255)
                                                     gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                            >  2:877256/1‑100 (MQ=255)
                                                     gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                            >  1:112112/1‑100 (MQ=255)
                                                     gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                            >  2:877090/1‑100 (MQ=255)
                                                       aTGGTGGATGCTGCGCACCTTTGGTGTATAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagatcggaagag                                          >  1:2843753/1‑91 (MQ=255)
                                                        tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                         >  1:2130087/1‑100 (MQ=255)
                                                        tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                        >  2:1159887/1‑100 (MQ=255)
                                                        tggtggAGGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                         <  1:2099470/100‑1 (MQ=255)
                                                            ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                     >  2:352646/1‑100 (MQ=255)
                                                              aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                   <  2:1169950/100‑1 (MQ=255)
                                                                 ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                >  2:1406124/1‑100 (MQ=255)
                                                                  tgcGCACCTGTGGTGTAGAGAAAGTGTCGATGCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                >  1:1124123/1‑99 (MQ=255)
                                                                   gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                              <  2:865055/100‑1 (MQ=255)
                                                                    cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                             <  2:2785096/100‑1 (MQ=255)
                                                                      cACCTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           <  2:1649435/100‑1 (MQ=255)
                                                                         cTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGACGGAAg                         >  1:1665605/1‑100 (MQ=255)
                                                                         cTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCTGGAAgg                        >  1:2391470/1‑100 (MQ=255)
                                                                         cTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        <  2:2786/100‑1 (MQ=255)
                                                                         cTTTGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        <  2:1494273/100‑1 (MQ=255)
                                                                           ttGGTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                      <  1:2301669/100‑1 (MQ=255)
                                                                             ggTGTAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCACTAGAGCTGCCGGAAGGagag                    <  1:1783737/100‑1 (MQ=255)
                                                                              gtgtAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                   <  2:2211069/100‑1 (MQ=255)
                                                                              gtgtAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                   >  2:2284045/1‑100 (MQ=255)
                                                                               tgtAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                  >  2:2109224/1‑100 (MQ=255)
                                                                               tgtAGAGAAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                  <  2:1105314/100‑1 (MQ=255)
                                                                               tggagagAAAGTGTCGATTCT‑‑GGGTG‑GGTGCACTTGCAGGCTGGCAGCGCGATGATCTGCGGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                   <  1:2900002/96‑1 (MQ=255)
                                                                                 tAGAGAAAGTGTCGATTCT‑‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                 >  1:1565559/1‑98 (MQ=255)
                                                                                   gagaAAGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              >  2:2819205/1‑100 (MQ=255)
                                                                                   gagaAAGTGTCGATTCT‑‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg             >  1:89307/1‑100 (MQ=255)
                                                                                   gagaAAGTGTCGATTCT‑‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg             <  2:1803887/100‑1 (MQ=255)
                                                                                       aaGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccg          >  2:829553/1‑100 (MQ=255)
                                                                                       aaGTGTCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccg          >  2:1469790/1‑100 (MQ=255)
                                                                                         gtgtCGATTCT‑‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGt       >  2:1559460/1‑100 (MQ=255)
                                                                                          tgtCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGt       >  2:822335/1‑100 (MQ=255)
                                                                                          tgtCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTCAGAAGAAGGGGCAGTGGAGCTGCCAGAAGGAGAGTTTAACGCCGCGt       >  1:1963003/1‑100 (MQ=255)
                                                                                           gtCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      <  2:2691238/100‑1 (MQ=255)
                                                                                           gtCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGAAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  2:2754318/1‑100 (MQ=255)
                                                                                            tCGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt     <  2:664440/100‑1 (MQ=255)
                                                                                             cGATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    >  2:1483260/1‑100 (MQ=255)
                                                                                              gATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa   >  2:2670441/1‑100 (MQ=255)
                                                                                              gATTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTACCGGAAGGAGAGTTTAACGCCGCGTTTa    <  2:2542404/99‑1 (MQ=255)
                                                                                               aTTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGGTGTTAGAGGAGGGTGCAGTAGAGCTGCCGGAAGGAGTGTTTAACGCCGCGTTTaa   >  1:1734358/1‑99 (MQ=255)
                                                                                               aTTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt  <  1:2492290/100‑1 (MQ=255)
                                                                                               aTTCTG‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt  <  1:2591129/100‑1 (MQ=255)
                                                                                                    |                                                                                                 
AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑‑GGG‑G‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652750‑2652943

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: ‑