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.G97.8% 325.0 / 0.8 89coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (46/41);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 101 total aligned reads displayed.

GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTT  >  NC_000913/2652751‑2652940
                                                                                                   |                                                                                             
ggATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g                                                                                              >  2:2114470/1‑100 (MQ=255)
 gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg                                                                                             <  1:2364923/100‑1 (MQ=255)
   tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg                                                                                           <  2:1915815/100‑1 (MQ=255)
   tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg                                                                                           <  2:463521/100‑1 (MQ=255)
   tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg                                                                                           >  2:2344728/1‑100 (MQ=255)
   tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg                                                                                           >  1:926565/1‑100 (MQ=255)
   tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg                                                                                           <  2:1300601/100‑1 (MQ=255)
    aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg‑g                                                                                         >  1:2588027/1‑100 (MQ=255)
    aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg‑g                                                                                         <  2:2603687/100‑1 (MQ=255)
       cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg‑gggt                                                                                      >  1:102223/1‑98 (MQ=255)
       cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg‑gggt                                                                                      >  1:1659672/1‑98 (MQ=255)
        aCCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg‑gggtg                                                                                     <  2:1686753/100‑4 (MQ=255)
           tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACTGTTGGTGTAGAGAAAGTGGCGATGCTGGGGG‑GG‑TGGAc                                                                                  >  1:930508/1‑100 (MQ=255)
            gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACt                                                                                 >  2:2644843/1‑100 (MQ=255)
              ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTg                                                                               >  1:79207/1‑100 (MQ=255)
              ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTg                                                                               >  2:2246804/1‑100 (MQ=255)
               tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGc                                                                              <  1:2469481/100‑1 (MQ=255)
                gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCa                                                                             <  1:1158335/100‑1 (MQ=255)
                gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCa                                                                             <  1:550707/100‑1 (MQ=255)
                gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCa                                                                             <  2:2160931/100‑1 (MQ=255)
                 tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAg                                                                            >  2:1782605/1‑100 (MQ=255)
                  cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAgg                                                                           >  1:2695258/1‑100 (MQ=255)
                  cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAgg                                                                           >  1:1451362/1‑100 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGc                                                                          >  1:1599768/1‑100 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGGTGG‑TGGACTTGCAGGc                                                                          >  2:2837760/1‑100 (MQ=255)
                    aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCt                                                                         >  1:1374166/1‑100 (MQ=255)
                    aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCt                                                                         >  2:1076957/1‑100 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTg                                                                        >  2:1520492/1‑100 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTg                                                                        >  1:2701786/1‑100 (MQ=255)
                      gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTgg                                                                       >  2:2858358/1‑100 (MQ=255)
                       aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGc                                                                      <  2:2187552/100‑1 (MQ=255)
                       aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGc                                                                      <  1:1339373/100‑1 (MQ=255)
                       aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGc                                                                      >  2:2360098/1‑100 (MQ=255)
                        cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCa                                                                     >  1:476826/1‑100 (MQ=255)
                         gAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAg                                                                    <  2:23427/100‑1 (MQ=255)
                          aaGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGCCGATTCTGGGGG‑GG‑GGGACTTGCAGGGTGGCAgc                                                                   >  2:1415236/1‑100 (MQ=255)
                           aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAgcg                                                                  >  1:1507245/1‑100 (MQ=255)
                            ggTAATCTTTTCTCATCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAgcgc                                                                 <  2:2840164/100‑1 (MQ=255)
                            ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAgcgc                                                                 <  2:1398404/100‑1 (MQ=255)
                            ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAgcgc                                                                 <  1:437075/100‑1 (MQ=255)
                             gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGcgcg                                                                >  1:1138800/1‑100 (MQ=255)
                             gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGcgcg                                                                <  1:1959134/100‑1 (MQ=255)
                              tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCga                                                               <  2:1374174/100‑1 (MQ=255)
                                aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatg                                                             >  2:2239707/1‑100 (MQ=255)
                                  cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTGTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatgat                                                           >  2:2027145/1‑100 (MQ=255)
                                  cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatgat                                                           >  2:1160272/1‑100 (MQ=255)
                                  cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTGTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatgat                                                           >  1:1478538/1‑100 (MQ=255)
                                       ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                      >  2:2809543/1‑100 (MQ=255)
                                       ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgcc                                                      <  2:102223/100‑2 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                   >  2:786546/1‑100 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                   <  1:1090403/100‑1 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                   >  1:2329886/1‑100 (MQ=255)
                                           gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                  <  1:1659247/100‑1 (MQ=255)
                                           gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                  <  1:1660108/100‑1 (MQ=255)
                                           gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                  >  1:1264325/1‑100 (MQ=255)
                                           gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                  <  2:780961/100‑1 (MQ=255)
                                           gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                  <  2:610332/100‑1 (MQ=255)
                                            ccccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                 <  1:1520480/100‑1 (MQ=255)
                                                aCGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaaga                                             <  1:2560602/100‑1 (MQ=255)
                                                 cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                            >  1:821135/1‑100 (MQ=255)
                                                 cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                            >  2:570979/1‑100 (MQ=255)
                                                   aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                          >  2:2725308/1‑100 (MQ=255)
                                                   aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                          <  1:2837734/100‑1 (MQ=255)
                                                   aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                          <  2:1274358/100‑1 (MQ=255)
                                                    gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                         >  2:1330668/1‑100 (MQ=255)
                                                      aTGGTGGATGCTGCGCCCCTTTGGTGTAGAGAAAGTGGCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                       <  1:1342791/100‑1 (MQ=255)
                                                      aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                       >  2:1840649/1‑100 (MQ=255)
                                                      aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                       <  2:538678/100‑1 (MQ=255)
                                                      aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                       <  1:2163070/100‑1 (MQ=255)
                                                        ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                     >  1:254886/1‑100 (MQ=255)
                                                         gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                    <  2:1451373/100‑1 (MQ=255)
                                                           ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                  <  2:327198/100‑1 (MQ=255)
                                                             aTGCTGCGCACCTGTGGTGTAGGGAAAGAGTCGATTCT‑GGGG‑GGTTGGAGTTGGATTCTGGCGGCGCGATGATCTGGTGTTAGAAGAAGGTGCCGAagag                                >  1:958946/1‑100 (MQ=255)
                                                                ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCGGACTGGCAGCGCGCTGATCTGCTGTTAGCAGAAGGTGCCGTAGAGc                               >  2:695389/1‑98 (MQ=255)
                                                                  gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                           >  2:381405/1‑100 (MQ=255)
                                                                   cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                          <  1:2344710/100‑1 (MQ=255)
                                                                   cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                          <  1:2881052/100‑1 (MQ=255)
                                                                     cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCgg                         >  2:2490039/1‑99 (MQ=255)
                                                                     cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                        >  2:2134456/1‑100 (MQ=255)
                                                                       ccTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                      <  1:786543/100‑1 (MQ=255)
                                                                       ccTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                      <  2:1278805/100‑1 (MQ=255)
                                                                       ccTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                      <  2:2715040/100‑1 (MQ=255)
                                                                        cTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                     <  1:2239689/100‑1 (MQ=255)
                                                                        cTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                     >  1:2203115/1‑100 (MQ=255)
                                                                           tGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgaga                  <  1:1817762/100‑1 (MQ=255)
                                                                             gtgtAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                <  2:1659684/100‑1 (MQ=255)
                                                                               gtAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt              <  2:2701810/100‑1 (MQ=255)
                                                                               gtAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt              <  1:1160267/100‑1 (MQ=255)
                                                                                 tcttccgatctCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa            <  2:484244/90‑1 (MQ=255)
                                                                                 agagAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa            >  1:2387542/1‑100 (MQ=255)
                                                                                  gagaAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGCTCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc           >  1:152904/1‑100 (MQ=255)
                                                                                   agaAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg          <  1:2114452/100‑1 (MQ=255)
                                                                                   agaAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCAGAAGGAGAGTATaatg          >  1:1689009/1‑98 (MQ=255)
                                                                                      aaGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccg       >  1:1317944/1‑100 (MQ=255)
                                                                                       aGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccgc      >  2:788605/1‑100 (MQ=255)
                                                                                        gtgtCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg     >  2:306704/1‑100 (MQ=255)
                                                                                         tgtCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGt    >  1:547943/1‑100 (MQ=255)
                                                                                          gtCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt   >  1:2626965/1‑100 (MQ=255)
                                                                                           tCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt  >  1:1419749/1‑100 (MQ=255)
                                                                                           tCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt  <  2:1478550/100‑1 (MQ=255)
                                                                                                   |                                                                                             
GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTT  >  NC_000913/2652751‑2652940

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