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.G95.5% 378.6 / 10.4 110coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base G (48/57);  total (52/58)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑GG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652750‑2652943
                                                                                                     |                                                                                                
aGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑Ct                                                                                                   >  2:1540393/1‑100 (MQ=255)
 ggATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑g                                                                                                 >  1:2282046/1‑100 (MQ=255)
  gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑gg                                                                                                <  2:1382266/100‑1 (MQ=255)
  gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑gg                                                                                                <  2:333975/100‑1 (MQ=255)
         aCCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑Gg‑ggg‑gggg                                                                                       >  2:2797490/1‑100 (MQ=255)
           cTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑gg‑ggg‑ggtgga                                                                                     <  1:1029310/100‑6 (MQ=255)
            tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑gg‑ggg‑gggggac                                                                                    >  2:1086924/1‑94 (MQ=255)
             gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACt                                                                                    >  1:32841/1‑100 (MQ=255)
             gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACt                                                                                    <  2:1686593/100‑1 (MQ=255)
             gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGACCATGGTGGATGCTGCGCAACTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACt                                                                                    >  2:804764/1‑100 (MQ=255)
                 gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCa                                                                                <  1:1670040/100‑1 (MQ=255)
                 gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCa                                                                                >  1:2811859/1‑100 (MQ=255)
                   cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAgg                                                                              <  2:1849614/100‑1 (MQ=255)
                   cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAgg                                                                              <  1:1229176/100‑1 (MQ=255)
                   cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑GG‑GGG‑GTGGACTTGCAGGc                                                                             <  1:998110/100‑1 (MQ=255)
                    tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGc                                                                             <  1:2111541/100‑1 (MQ=255)
                    tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGc                                                                             >  2:821796/1‑100 (MQ=255)
                    tATGACGAAGGTAATCTTTTCTCAGCACCAAGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAgg                                                                              <  2:130671/99‑1 (MQ=255)
                     aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCt                                                                            >  2:1164807/1‑100 (MQ=255)
                     aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCt                                                                            <  1:842561/100‑1 (MQ=255)
                     aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCt                                                                            <  1:552496/100‑1 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTg                                                                           <  1:2118146/100‑1 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTg                                                                           <  1:2118102/100‑1 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTgg                                                                          >  1:1919444/1‑100 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTgg                                                                          >  1:1406987/1‑100 (MQ=255)
                        aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGc                                                                         <  2:2481024/100‑1 (MQ=255)
                         cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCa                                                                        >  1:427507/1‑100 (MQ=255)
                         cGAAGGTAATCTTTTCTCAGCCCCACGAGCAGGGTGGATTCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑GGCGGG‑GTGGACTTGCAGGCTGGCa                                                                        >  1:2045787/1‑100 (MQ=255)
                            aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAgcg                                                                     <  1:1238770/100‑1 (MQ=255)
                             ggTAATATTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAgcg                                                                     >  1:1691256/1‑99 (MQ=255)
                              gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGcgcg                                                                   <  1:1525946/100‑1 (MQ=255)
                               tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGcgcg                                                                   >  1:645553/1‑99 (MQ=255)
                               tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCga                                                                  <  1:1086908/100‑1 (MQ=255)
                                 aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCgatg                                                                <  1:2311379/100‑1 (MQ=255)
                                 aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCgatg                                                                <  2:1735389/100‑1 (MQ=255)
                                  tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCgatga                                                               >  2:1580524/1‑100 (MQ=255)
                                   nttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCgatgat                                                              >  1:450485/2‑100 (MQ=255)
                                   cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCgatgat                                                              >  2:2407622/1‑100 (MQ=255)
                                   cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCgatgat                                                              >  2:1534115/1‑100 (MQ=255)
                                   cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCgatgat                                                              >  2:1534123/1‑100 (MQ=255)
                                   cTTTTCTCAGCCCCACGAGCATGGTGCATGCTGCGCACCTTTGGTGTAGCGACAGTGTCGCTTCCG‑GG‑GGG‑GTGGACTTGCACGCTGGCAGCGCgatgat                                                              >  2:903428/1‑100 (MQ=255)
                                       tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                          <  1:1237811/100‑1 (MQ=255)
                                         tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                        <  1:35713/100‑1 (MQ=255)
                                         tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                        >  2:456086/1‑100 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  1:815750/99‑1 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     <  1:456081/100‑1 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     <  1:2766513/100‑1 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     <  2:865656/100‑1 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     >  2:2299786/1‑100 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     <  1:1540369/100‑1 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     <  2:555137/100‑1 (MQ=255)
                                             ccccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                    <  1:1336440/100‑1 (MQ=255)
                                              ctcACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   <  1:1764893/98‑1 (MQ=255)
                                              cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   <  1:504489/100‑1 (MQ=255)
                                               ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                  <  2:2886532/100‑1 (MQ=255)
                                                cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAAT‑CTGGG‑GGG‑GAGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                 >  2:397397/1‑100 (MQ=255)
                                                  cGAGCATGGTGGATGCTGCGCACCTTAGGTGGAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGGAGCGCGATGATCTGCTGTTggaag                                                 >  1:1825951/1‑98 (MQ=255)
                                                   gAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGaagaag                                              <  2:838268/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  1:274300/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             >  1:683732/1‑100 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:1407783/100‑1 (MQ=255)
                                                        tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                         >  1:2472817/1‑100 (MQ=255)
                                                        tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                         >  2:1921735/1‑100 (MQ=255)
                                                        tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                         >  2:1431231/1‑100 (MQ=255)
                                                          gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       <  1:621641/100‑1 (MQ=255)
                                                          gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       <  2:1526100/100‑1 (MQ=255)
                                                          gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       >  2:2860843/1‑100 (MQ=255)
                                                           tggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTa                                      <  1:1536486/100‑1 (MQ=255)
                                                            ngATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                     >  2:2247364/2‑100 (MQ=255)
                                                            ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                     <  1:2726321/100‑1 (MQ=255)
                                                             gATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                    <  2:1013563/100‑1 (MQ=255)
                                                              aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                   <  1:872575/100‑1 (MQ=255)
                                                               tgctgcGCACCTTTGGTGTAGAGTAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATTATCTGCTGTTAGAAGAAGGTGCAGTag                                     <  1:1002363/97‑1 (MQ=255)
                                                                 ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                <  1:903420/100‑1 (MQ=255)
                                                                  tgcGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                               >  2:1946289/1‑100 (MQ=255)
                                                                  tgcGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                               >  2:2271920/1‑100 (MQ=255)
                                                                  tgcGCACATTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                               >  1:2788068/1‑100 (MQ=255)
                                                                   gcgcACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                              <  1:877332/100‑1 (MQ=255)
                                                                     gcACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCgg                            <  1:1796151/100‑1 (MQ=255)
                                                                      cACCTTTGGTGTAGAGATAGTGGCGCTT‑CA‑GG‑GGGTTTGGACTTGCAGGCTGGCAGTGCGATGTTCAGACGTAAGAAAAAGGTGAAGTATAGCTGc                               >  1:1104840/1‑96 (MQ=255)
                                                                           ttGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                      <  1:111243/100‑1 (MQ=255)
                                                                           ttGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                      <  2:544945/100‑1 (MQ=255)
                                                                           ttGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                      <  1:2075808/100‑1 (MQ=255)
                                                                           ttGGTGTAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                      <  1:2292203/100‑1 (MQ=255)
                                                                              gtgtAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                   <  1:1376065/100‑1 (MQ=255)
                                                                               tgtAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                  >  1:1764024/1‑100 (MQ=255)
                                                                               tgtAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGATGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                  >  2:141325/1‑100 (MQ=255)
                                                                                ttAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                 <  2:2723197/99‑1 (MQ=255)
                                                                                gtAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                  >  2:1094224/1‑99 (MQ=255)
                                                                                 tAGAGAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa                >  1:2723161/1‑100 (MQ=255)
                                                                                  agagAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTTCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  2:2473556/1‑100 (MQ=255)
                                                                                  agagAAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  1:1040499/1‑100 (MQ=255)
                                                                                   gagaAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              >  1:1656542/1‑100 (MQ=255)
                                                                                   gagaAAGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              <  1:195614/100‑1 (MQ=255)
                                                                                      aaaGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgcc           >  2:1696836/1‑100 (MQ=255)
                                                                                       aaGTGTCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccg          >  2:451193/1‑100 (MQ=255)
                                                                                         gtgtCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg        <  1:2860805/100‑1 (MQ=255)
                                                                                           gtCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  2:2829954/1‑100 (MQ=255)
                                                                                            tCGATT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt     >  2:1600508/1‑100 (MQ=255)
                                                                                               aTT‑CTGGG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGCTGAACGCCGCGTTTAAt  >  2:2863400/1‑100 (MQ=255)
                                                                                                     |                                                                                                
AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑GG‑GGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652750‑2652943

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