Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,652,855 (G)6→7 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.G91.4% 315.4 / 24.3 93coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (6/2);  new base G (43/42);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑G‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTG  >  NC_000913/2652752‑2652947
                                                                                                  |                                                                                                     
gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g                                                                                                      <  1:1284743/99‑1 (MQ=255)
   aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g‑gg‑gg                                                                                                <  1:288882/100‑1 (MQ=255)
    aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g‑gg‑gg‑g                                                                                              >  2:1815532/1‑100 (MQ=255)
      cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g‑gg‑gg‑ggt                                                                                            <  1:616340/100‑3 (MQ=255)
       aCCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g‑gg‑gg‑ggt                                                                                            <  1:1371454/99‑3 (MQ=255)
        ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g‑gg‑gg‑gg                                                                                             >  1:897493/1‑96 (MQ=255)
           gATTGTCTATGACGAAGGTAATCTTTTCTCAGTCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACt                                                                                        <  1:2597149/100‑1 (MQ=255)
            aTTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACtt                                                                                       >  2:2131779/1‑100 (MQ=255)
             ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTg                                                                                      <  2:565949/100‑1 (MQ=255)
              tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGc                                                                                     <  2:2710279/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCa                                                                                    <  2:855657/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCa                                                                                    <  2:1160544/100‑1 (MQ=255)
                tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAg                                                                                   <  2:2836485/100‑1 (MQ=255)
                tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAg                                                                                   >  2:654042/1‑100 (MQ=255)
                 ttccgatctaGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAgg                                                                                  <  2:2766960/91‑1 (MQ=255)
                 cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAgg                                                                                  >  2:2790139/1‑100 (MQ=255)
                  tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGc                                                                                 >  2:1115618/1‑100 (MQ=255)
                  tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGc                                                                                 <  1:638648/100‑1 (MQ=255)
                  tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGc                                                                                 <  1:1151081/100‑1 (MQ=255)
                  tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGc                                                                                 <  1:1592731/100‑1 (MQ=255)
                     gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTgg                                                                              >  2:1055028/1‑100 (MQ=255)
                      aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTCTGGTGTAGAGAAAGGGTCGATTCT‑G‑GG‑GGTGTGGACTTGCAGGCTgtc                                                                             >  1:1018395/1‑98 (MQ=255)
                          aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAggagatcggaa                                                                         >  1:2766944/1‑91 (MQ=255)
                             tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGTTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCga                                                                      <  2:1018399/100‑1 (MQ=255)
                             tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGcgcg                                                                       >  2:881339/1‑99 (MQ=255)
                             tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCga                                                                      >  2:517469/1‑100 (MQ=255)
                             tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑G‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCgat                                                                     >  2:1919102/1‑100 (MQ=255)
                              aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCgat                                                                     >  2:2547202/1‑100 (MQ=255)
                               aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCgatg                                                                    <  1:1946785/100‑1 (MQ=255)
                               aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCgatg                                                                    >  2:1586597/1‑100 (MQ=255)
                                  ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATc                                                                 <  2:2731909/100‑1 (MQ=255)
                                  ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATc                                                                 <  2:2906352/100‑1 (MQ=255)
                                  ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATc                                                                 >  2:1218303/1‑100 (MQ=255)
                                      ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                             >  2:373107/1‑100 (MQ=255)
                                       tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                            <  2:1797115/100‑1 (MQ=255)
                                       tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                            <  2:1796229/100‑1 (MQ=255)
                                       tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                            <  2:1279347/100‑1 (MQ=255)
                                        cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑G‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                          >  2:1471471/1‑100 (MQ=255)
                                         aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                          <  1:624648/100‑1 (MQ=255)
                                         aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                          <  1:2216734/100‑1 (MQ=255)
                                         aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                          <  1:70203/100‑1 (MQ=255)
                                         aGCCCCACGAGCATGGTGGATGCTGCGCACCGTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGTCTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                          <  1:1431048/100‑1 (MQ=255)
                                           ccccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                        <  1:2761487/100‑1 (MQ=255)
                                            cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                       <  2:2919938/100‑1 (MQ=255)
                                            cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACCTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                       <  1:1843849/100‑1 (MQ=255)
                                             ccACGAGCCTGGTGGATGCTGCGCACCTTTGTTGTAGAGAAAGTGTCGATTCG‑G‑GGTGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCGGCTGTGagaa                                                      >  1:2556572/1‑100 (MQ=255)
                                             ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                      >  1:2482957/1‑100 (MQ=255)
                                             ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                      >  1:1946502/1‑100 (MQ=255)
                                                cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                                   >  1:790084/1‑100 (MQ=255)
                                                cGAGCATCGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACGTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                                   <  1:1647069/100‑1 (MQ=255)
                                                  aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                                 <  2:1693114/100‑1 (MQ=255)
                                                    ttccgatcttgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑G‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTgn                                              <  1:626735/91‑2 (MQ=255)
                                                     aTGGTGGATGCTGCGCGCCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                              <  1:1838860/100‑1 (MQ=255)
                                                     aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                              <  1:698198/100‑1 (MQ=255)
                                                      tggtggATGCTGCGCAGCTTTGGTGTAGAGAAAGTTTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                             <  1:1796627/100‑1 (MQ=255)
                                                      tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                             >  2:534656/1‑100 (MQ=255)
                                                      tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                             >  2:778094/1‑100 (MQ=255)
                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGGAGAAGGTGCAg                                            >  1:1903784/1‑100 (MQ=255)
                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                            <  1:584723/100‑1 (MQ=255)
                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                            >  2:338206/1‑100 (MQ=255)
                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                            >  1:912165/1‑100 (MQ=255)
                                                        gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                           <  1:191870/100‑1 (MQ=255)
                                                          ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                         <  1:2661060/100‑1 (MQ=255)
                                                           gATGCTGCGCACCTTTGGTGTAGAGCAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                         >  2:2921575/1‑99 (MQ=255)
                                                           gATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                        >  1:2497818/1‑100 (MQ=255)
                                                            aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                       <  1:6627/100‑1 (MQ=255)
                                                            aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                       <  2:2497832/100‑1 (MQ=255)
                                                             tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                       <  1:2321288/99‑1 (MQ=255)
                                                             ngctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑G‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAgatcggaa                                     >  2:626736/2‑93 (MQ=255)
                                                              gctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGCTGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCt                                     <  1:1847835/100‑1 (MQ=255)
                                                                tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                                   >  1:959184/1‑100 (MQ=255)
                                                                   gcACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCgg                                >  2:354127/1‑100 (MQ=255)
                                                                    cACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                               >  2:966270/1‑100 (MQ=255)
                                                                    cACCTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                               >  1:907152/1‑100 (MQ=255)
                                                                       cTTTGGTGTAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                            >  1:2243480/1‑100 (MQ=255)
                                                                            gtgtAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                       >  2:2409471/1‑100 (MQ=255)
                                                                            gtgtAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                       <  2:790086/100‑1 (MQ=255)
                                                                             tgtAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                      <  2:795252/100‑1 (MQ=255)
                                                                              gtAGAGAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                     >  2:1957117/1‑100 (MQ=255)
                                                                               taaagaAAGTCTTCATTCT‑GAGA‑GG‑CTGTAATTTCAGGCTGGCTGCGAGATGATCTGCTGTTAGAAGAAGGTGCAGTAAAGCGGCTGGAAGGAGAGTGTa                    <  1:3012941/97‑1 (MQ=255)
                                                                                agagAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa                   >  2:2792050/1‑100 (MQ=255)
                                                                                agagAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa                   >  2:948672/1‑100 (MQ=255)
                                                                                agagAAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa                   >  2:2143608/1‑100 (MQ=255)
                                                                                 gtgaAAGTGTCGACTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc                  >  2:2074476/3‑100 (MQ=255)
                                                                                 gagaAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc                  <  2:2013283/100‑1 (MQ=255)
                                                                                 gagaAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc                  >  2:2790470/1‑100 (MQ=255)
                                                                                 gagaAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc                  >  2:2593102/1‑100 (MQ=255)
                                                                                  agaAAGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg                 >  1:942163/1‑100 (MQ=255)
                                                                                   gaAAGTGTCGAATCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc                >  1:995350/1‑100 (MQ=255)
                                                                                      aGTGTCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccgc             >  2:522232/1‑100 (MQ=255)
                                                                                       gtgtCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg            <  2:1703102/100‑1 (MQ=255)
                                                                                         gtCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCTATGATCTGCTGTTAGGAGAAGGTGCCGTAGAGCTGCCGGAAGGAGAGATTAACGCCGCGtt          >  1:86387/1‑100 (MQ=255)
                                                                                         gtCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt          >  1:2840319/1‑100 (MQ=255)
                                                                                         gtCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt          >  1:29587/1‑100 (MQ=255)
                                                                                         gtCGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt          <  1:1055024/100‑1 (MQ=255)
                                                                                           cGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa        <  1:354126/100‑1 (MQ=255)
                                                                                           cGATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa        >  2:1430373/1‑100 (MQ=255)
                                                                                            gATTCTGG‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa       >  2:2389720/1‑100 (MQ=255)
                                                                                                cT‑G‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTg  >  1:259051/1‑100 (MQ=255)
                                                                                                  |                                                                                                     
GATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑G‑GG‑GG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTG  >  NC_000913/2652752‑2652947

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