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.G94.5% 385.2 / 13.9 110coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base G (55/49);  total (59/51)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑GGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAA  >  NC_000913/2652750‑2652942
                                                                                                     |                                                                                              
aGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑t                                                                                                 >  2:3205889/1‑100 (MQ=255)
 ggATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g                                                                                               <  2:2108744/100‑1 (MQ=255)
 ggATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g                                                                                               >  1:187016/1‑100 (MQ=255)
  gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gg                                                                                              <  2:774481/100‑1 (MQ=255)
  gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gg                                                                                              <  1:55373/100‑1 (MQ=255)
   aTAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑ggg                                                                                             <  2:618144/100‑1 (MQ=255)
   aTAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑ggg                                                                                             >  1:557697/1‑100 (MQ=255)
    tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gggg                                                                                            >  2:2351597/1‑100 (MQ=255)
     aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gggg‑g                                                                                          <  1:3051392/100‑1 (MQ=255)
       gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gggg‑ggg                                                                                        >  2:3164267/1‑99 (MQ=255)
       gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gggg‑ggg                                                                                        <  2:2068090/100‑2 (MQ=255)
       gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gggg‑ggg                                                                                        >  1:1613616/1‑99 (MQ=255)
        ccccTTATTGTCTATGACGAAGGTAATCTTTTCTAAGAAACAAGAACATGGTGGATGTCGCGCAACTTTGGTGAAGAGAAATAGTCGATTC‑T‑gg                                                                                              >  2:1480648/3‑94 (MQ=255)
         aCCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gggg‑gggtg                                                                                      >  2:1145986/1‑97 (MQ=255)
         aCCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gggg‑gggtg                                                                                      >  1:339321/1‑97 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gggg‑gggtgg                                                                                     <  2:1928002/100‑5 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gggg‑gggtgg                                                                                     <  1:2528826/100‑5 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gggg‑gggtgg                                                                                     <  1:2100877/100‑5 (MQ=255)
            tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gggg‑gggtgga                                                                                    >  1:48025/1‑94 (MQ=255)
            tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gggg‑gggtgga                                                                                    >  2:2022385/1‑94 (MQ=255)
             gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gggg‑gtgggg                                                                                     >  1:1878607/1‑92 (MQ=255)
               ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTg                                                                                 >  1:29015/1‑100 (MQ=255)
               atctCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTg                                                                                 <  1:1254227/97‑1 (MQ=255)
                tGTCGATGAGGAAGGTAATCTTTTCTCAGTGCCACGAGCATGGTGGATGCTGCGCACCTGTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTg                                                                                 <  2:711479/99‑1 (MQ=255)
                 gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCa                                                                               <  2:557695/100‑1 (MQ=255)
                 gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCATTTGTGTAGGGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCa                                                                               >  2:2280674/1‑100 (MQ=255)
                   cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAgg                                                                             >  2:419158/1‑100 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTg                                                                          <  2:335895/100‑1 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTg                                                                          >  1:1561424/1‑100 (MQ=255)
                      ngACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTg                                                                          >  2:382684/2‑100 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTgg                                                                         >  1:736540/1‑100 (MQ=255)
                        aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGc                                                                        >  2:778555/1‑100 (MQ=255)
                           aaGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAgc                                                                     >  2:3113592/1‑100 (MQ=255)
                             ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAgcgc                                                                   >  2:655662/1‑100 (MQ=255)
                             ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑GGGGCGGTGGACTTGCAGGCTGGCAgcg                                                                    >  1:863634/1‑99 (MQ=255)
                              gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCAGCGCACCTTTGGAGTAGAGAAACTGTGAAATG‑T‑GGGGAGATGGACTTGAAGGCTAGCCgc                                                                     >  2:2950891/1‑97 (MQ=255)
                                aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCgat                                                                >  1:11237/1‑100 (MQ=255)
                                 aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCgatg                                                               >  1:2013598/1‑100 (MQ=255)
                                 aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCgatg                                                               >  2:1287489/1‑100 (MQ=255)
                                 aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCgatg                                                               <  1:914806/100‑1 (MQ=255)
                                     tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATct                                                           >  2:2101418/1‑100 (MQ=255)
                                      ttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATctg                                                          <  1:1267813/100‑1 (MQ=255)
                                       tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                         >  2:583267/1‑100 (MQ=255)
                                       tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                         <  1:1993373/100‑1 (MQ=255)
                                        ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                        <  1:1529652/100‑1 (MQ=255)
                                           aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     <  1:917808/100‑1 (MQ=255)
                                           aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑GGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    <  1:3111610/100‑1 (MQ=255)
                                              cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                  >  2:304530/1‑100 (MQ=255)
                                              cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                  >  1:1354381/1‑100 (MQ=255)
                                              cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                  >  1:2470476/1‑100 (MQ=255)
                                                cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                >  2:2643305/1‑100 (MQ=255)
                                                  cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                              <  1:2073728/100‑1 (MQ=255)
                                                   gAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGaagaag                                             <  1:1552218/100‑1 (MQ=255)
                                                     gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                           >  1:1075107/1‑100 (MQ=255)
                                                      cATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTg                                          <  2:3248572/100‑1 (MQ=255)
                                                       aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                           >  2:1292821/1‑98 (MQ=255)
                                                       aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                         <  2:1684299/100‑1 (MQ=255)
                                                       aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                         <  1:2171251/100‑1 (MQ=255)
                                                       aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                         <  1:2643326/100‑1 (MQ=255)
                                                       aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGATGTTAGAAGAAGGTg                                          >  2:2225616/1‑99 (MQ=255)
                                                        tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                        >  1:290251/1‑100 (MQ=255)
                                                         ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                       <  1:1574150/100‑1 (MQ=255)
                                                         ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                       <  1:1145991/100‑1 (MQ=255)
                                                          gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                      <  1:3164295/100‑1 (MQ=255)
                                                           tggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGATGCAGTa                                     >  1:1064054/1‑100 (MQ=255)
                                                             gATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                   <  2:474891/100‑1 (MQ=255)
                                                              aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                  >  1:313130/1‑100 (MQ=255)
                                                              aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                  >  2:683764/1‑100 (MQ=255)
                                                              aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTA‑TGGGGG‑GGGGGAATTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                   >  2:757944/1‑99 (MQ=255)
                                                                gctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCt                                <  1:3166036/100‑1 (MQ=255)
                                                                 ctgcGCACCTTTGGTGTAGGGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                               >  2:2868222/1‑100 (MQ=255)
                                                                  tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                              >  1:1823378/1‑100 (MQ=255)
                                                                  tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                              <  2:2767709/100‑1 (MQ=255)
                                                                   gcgcACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGTTGCAGTAGAGCTGcc                             >  1:621749/1‑100 (MQ=255)
                                                                   gcgcACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                             <  1:1044581/100‑1 (MQ=255)
                                                                    cgcACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                            <  1:2365011/100‑1 (MQ=255)
                                                                    cgcACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                            <  2:1683373/100‑1 (MQ=255)
                                                                        ccTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                        <  1:1287497/100‑1 (MQ=255)
                                                                        ccTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                        >  2:2901718/1‑100 (MQ=255)
                                                                         cTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                       >  1:2504850/1‑100 (MQ=255)
                                                                          tttGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                       <  2:398132/99‑1 (MQ=255)
                                                                          tttGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGga                      >  1:1992177/1‑100 (MQ=255)
                                                                            tGGTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgaga                    <  1:1449142/100‑1 (MQ=255)
                                                                             ggTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgaga                    >  1:2610456/1‑99 (MQ=255)
                                                                             ggTGTAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                   <  1:3051276/100‑1 (MQ=255)
                                                                                gtAGAGAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                <  1:2101433/100‑1 (MQ=255)
                                                                                  agagAAAGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa              <  2:1866934/100‑1 (MQ=255)
                                                                                     gaAAGTGTCGATTC‑TGGGGG‑GGTAGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc           >  2:1249211/1‑100 (MQ=255)
                                                                                      aaaGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgcc          <  1:655662/100‑1 (MQ=255)
                                                                                        aGTGTCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccgc        <  2:339319/100‑1 (MQ=255)
                                                                                        aGTGTAGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAAAGTTTAAcgcc          >  2:901101/1‑98 (MQ=255)
                                                                                          tgtCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGt      >  2:569157/1‑100 (MQ=255)
                                                                                           gtCGATTCTT‑GGGG‑GGTGGACTTACAGACTGGCAGCGCGATGATCTGCCGTTAGAAGAAGGTCCAGTGGGGATGCCGCAAGGAGAGGTTaaagatgtgt      >  2:2576761/1‑91 (MQ=255)
                                                                                           gtCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt     >  1:1487682/1‑100 (MQ=255)
                                                                                            tCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt    >  1:2749874/1‑100 (MQ=255)
                                                                                            tCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt    <  1:1492045/100‑1 (MQ=255)
                                                                                            tCGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt    <  1:1180460/100‑1 (MQ=255)
                                                                                             cGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa   >  1:1657573/1‑100 (MQ=255)
                                                                                             cGATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa   >  2:883552/1‑100 (MQ=255)
                                                                                              gATTC‑TGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGCAAGGAGAGTTTAACGCCGCGTTTaa  >  1:422431/1‑100 (MQ=255)
                                                                                                     |                                                                                              
AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑GGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAA  >  NC_000913/2652750‑2652942

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