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.G95.7% 483.0 / 12.5 140coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base G (62/72);  total (64/76)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 152 total aligned reads displayed.

GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑GGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652751‑2652943
                                                                                                    |                                                                                               
ggATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑g                                                                                                <  2:1228649/100‑1 (MQ=255)
   tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑gggg                                                                                             <  1:780643/100‑1 (MQ=255)
    aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑ggggg                                                                                            <  2:246172/100‑1 (MQ=255)
     aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑gggggg                                                                                           <  1:1070826/100‑1 (MQ=255)
     aGCACCTCATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCGCCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑gggggg                                                                                           >  1:1988655/1‑100 (MQ=255)
      gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCCCGAGCATGGGGGAGGCTGCGCACCGTTGGTGTCGAGAAAGGGGCGAAT‑CG‑gggtggg                                                                                          >  1:1784456/1‑99 (MQ=255)
      gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑ggggggg                                                                                          >  2:2724115/1‑99 (MQ=255)
       cACCTGGTTGTGTAGGACGTAGGTACTCTTTTCTCAGCGCCACGAGCATGGTGGATGCTGTGCACCTTTTGTGTAGTGGAAGTGTGGATT‑CT‑ggg                                                                                              <  1:2656789/95‑1 (MQ=255)
         ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑gggggggtg                                                                                        <  2:3299821/99‑4 (MQ=255)
            gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGGATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGAc                                                                                     <  1:2497068/99‑1 (MQ=255)
            gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACt                                                                                    >  1:21118/1‑100 (MQ=255)
               tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGGTGCGCACCTTTGTGGTAGAGAAAGTGTCGATT‑CT‑gggggggtgg                                                                                       <  1:311860/94‑5 (MQ=255)
                gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCATCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCa                                                                                <  1:2047788/100‑1 (MQ=255)
                gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCa                                                                                <  2:2972428/100‑1 (MQ=255)
                 tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAg                                                                               <  2:1805272/100‑1 (MQ=255)
                 tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAg                                                                               <  2:3711563/100‑1 (MQ=255)
                  cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTCGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAgg                                                                              >  2:2672014/1‑100 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGc                                                                             <  2:2702597/100‑1 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTg                                                                           <  1:2579270/100‑1 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTg                                                                           <  2:1063352/100‑1 (MQ=255)
                       aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGc                                                                         <  1:3379017/100‑1 (MQ=255)
                       aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGc                                                                         <  1:2778799/100‑1 (MQ=255)
                        cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCa                                                                        >  2:3727273/1‑100 (MQ=255)
                           aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAgcg                                                                     >  2:1613085/1‑100 (MQ=255)
                           aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAgcg                                                                     >  2:650834/1‑100 (MQ=255)
                            ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAgcgc                                                                    <  2:572433/100‑1 (MQ=255)
                             gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGcgcg                                                                   <  2:3293137/100‑1 (MQ=255)
                              tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGTTT‑CTGGGGGGG‑TGGACTTACAGGCTGGCAGcgcg                                                                   >  2:2811324/1‑99 (MQ=255)
                                 tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                               >  2:3584229/1‑100 (MQ=255)
                                 tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                               <  1:40110/100‑1 (MQ=255)
                                 tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                               >  2:1873091/1‑100 (MQ=255)
                                  cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCTCCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGGGCgatg                                                                >  1:3597507/1‑98 (MQ=255)
                                   ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATc                                                             <  1:371737/100‑1 (MQ=255)
                                   ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATc                                                             <  2:3769768/100‑1 (MQ=255)
                                      tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgn                                                          <  1:1820511/100‑2 (MQ=255)
                                       ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                         >  1:3226789/1‑100 (MQ=255)
                                        tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                        <  2:298440/100‑1 (MQ=255)
                                         cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGt                                                       <  1:1904927/100‑1 (MQ=255)
                                         cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGt                                                       <  1:1922717/100‑1 (MQ=255)
                                         cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGt                                                       <  2:1443934/100‑1 (MQ=255)
                                           gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAACGGGGCGAGT‑CG‑GGGGGGTTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     >  1:289346/1‑100 (MQ=255)
                                           gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      >  1:2526054/1‑99 (MQ=255)
                                           gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     >  1:1237303/1‑100 (MQ=255)
                                           gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     <  2:3539518/100‑1 (MQ=255)
                                            cctcACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGGGTCGATT‑CTGGGGGGG‑CGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                    <  1:2019559/97‑1 (MQ=255)
                                            ccccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                    <  2:2750616/100‑1 (MQ=255)
                                             cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   >  1:3794718/1‑100 (MQ=255)
                                             cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   <  2:3519557/100‑1 (MQ=255)
                                             cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   <  2:787299/100‑1 (MQ=255)
                                              ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                  <  1:923268/100‑1 (MQ=255)
                                                   aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  1:2968245/100‑1 (MQ=255)
                                                   aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:2822688/100‑1 (MQ=255)
                                                    gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGTAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTCTTATAAGAAgg                                             <  1:1992218/99‑1 (MQ=255)
                                                    gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                            >  2:3358373/1‑100 (MQ=255)
                                                     cATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTg                                           <  1:2595149/100‑1 (MQ=255)
                                                      aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                          <  1:1166611/100‑1 (MQ=255)
                                                        ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                         >  2:1694223/1‑99 (MQ=255)
                                                        ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                        >  2:894855/1‑100 (MQ=255)
                                                         gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       >  2:1796655/1‑100 (MQ=255)
                                                         gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       >  2:517782/1‑100 (MQ=255)
                                                         gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       <  1:82704/100‑1 (MQ=255)
                                                         gtggATGCGGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       <  1:3228606/100‑1 (MQ=255)
                                                           ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                     >  1:3454712/1‑100 (MQ=255)
                                                            gATGCTGCGCGCCTGTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTCGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                    <  2:2877951/100‑1 (MQ=255)
                                                              tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑GGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGTAGAAGGTGCAGGagag                                   >  2:2066932/1‑99 (MQ=255)
                                                               gctgcGCACCTTTGGTGTAGTGAAAGTGTCGATT‑CT‑GGGGGG‑TGGACTTGCAGGCTGGCAGAGTTATGATCTGCTGTTAGAAGAAGGTGCAGTACAGCTg                                <  1:937367/100‑1 (MQ=255)
                                                               gccgcgcACATTTGGTGTGGAGCAAGTGTCGATTGAG‑GGGGGG‑TGGAGTTGCAGGATGGCAGACCGATGCTCTGCTGTTAAAAGAAGGTGCAGTagtg                                   <  1:1745018/95‑3 (MQ=255)
                                                                   cgcACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                             <  2:41048/100‑1 (MQ=255)
                                                                    gcAGCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCgg                            <  1:2964190/100‑1 (MQ=255)
                                                                    gcACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑GGGGGG‑TGGACTTGCAGGCTGACAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           <  2:615293/100‑1 (MQ=255)
                                                                     cACCTTTGGTGGAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           <  1:2411737/100‑1 (MQ=255)
                                                                       ccTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                         <  2:797340/100‑1 (MQ=255)
                                                                        cTTTGGTGTAGAGCAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTCGCAGAAGGTGCAGTAGAGCTGcc                              >  1:3571609/1‑94 (MQ=255)
                                                                        cTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        >  1:1329324/1‑100 (MQ=255)
                                                                         tttGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGga                       <  2:2646418/100‑1 (MQ=255)
                                                                         tttGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGga                       >  1:2876093/1‑100 (MQ=255)
                                                                          ttGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                      >  2:33062/1‑100 (MQ=255)
                                                                          ttGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                      <  1:3727244/100‑1 (MQ=255)
                                                                          ttGGTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                      >  2:1451318/1‑100 (MQ=255)
                                                                            ggTGTAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    >  2:2123606/1‑100 (MQ=255)
                                                                             gtgtAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                   <  2:2599061/100‑1 (MQ=255)
                                                                               gtAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                 >  1:3414810/1‑100 (MQ=255)
                                                                               gtAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                 <  2:259490/100‑1 (MQ=255)
                                                                                tAGAGAAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGGTCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGCAGGAGTGTTTa                >  1:3463354/1‑100 (MQ=255)
                                                                                    gaAAGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc            >  1:237081/1‑100 (MQ=255)
                                                                                     aaaGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgcc           <  2:3224009/100‑1 (MQ=255)
                                                                                       aGTGTCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGTCGGAAGGAGAGTTTAAcgccg          >  1:2730479/1‑99 (MQ=255)
                                                                                        gtgtCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg        <  1:2206264/100‑1 (MQ=255)
                                                                                        gtgtCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg        >  2:1601049/1‑100 (MQ=255)
                                                                                         tgtCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGt       >  2:2791374/1‑100 (MQ=255)
                                                                                          gtCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGACGGAGAGTTTAACGCCGCGtt      >  1:1759128/1‑100 (MQ=255)
                                                                                          gtCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      <  2:51205/100‑1 (MQ=255)
                                                                                          gtCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCAGGAAGGAGAGTTTAACGCCGCGtt      >  1:3140173/1‑100 (MQ=255)
                                                                                           tCGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt     <  2:1055711/100‑1 (MQ=255)
                                                                                            cGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGGAGGAGAGTTTAAAGCCGCGTTTa    >  2:3484621/1‑100 (MQ=255)
                                                                                            cGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    >  2:548815/1‑100 (MQ=255)
                                                                                            cGATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    >  2:1880011/1‑100 (MQ=255)
                                                                                             gATT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa   >  2:1624870/1‑100 (MQ=255)
                                                                                              aTT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCCGCGCGATGATCGGCTGTTAGAAGAAGGTTCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt  <  2:3382352/100‑1 (MQ=255)
                                                                                              aTT‑CTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt  >  2:2262396/1‑100 (MQ=255)
                                                                                                    |                                                                                               
GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑GGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652751‑2652943

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