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.6% 312.0 / 6.8 90coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base G (43/43);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATC  >  NC_000913/2652752‑2652944
                                                                                                  |                                                                                                
gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg                                                                                                <  2:253696/100‑1 (MQ=255)
gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg                                                                                                <  2:2881456/100‑1 (MQ=255)
 aTAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑ggg                                                                                               <  2:1123012/100‑1 (MQ=255)
    aGCACCTGATTGTCTGTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggg                                                                                            <  1:2990099/100‑1 (MQ=255)
    aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggg                                                                                            >  2:627644/1‑100 (MQ=255)
     gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑ggggggg                                                                                           >  1:2939415/1‑99 (MQ=255)
     gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑ggggggg                                                                                           >  2:1849946/1‑99 (MQ=255)
        ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggggtgg                                                                                        >  2:1917190/1‑96 (MQ=255)
         cTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGGGG‑TGGAc                                                                                      <  1:1786001/100‑1 (MQ=255)
           gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACt                                                                                     <  1:1220376/100‑1 (MQ=255)
            aTTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGATCACGTGTCGATTCT‑gggggggtggcgtt                                                                                    >  2:511298/1‑92 (MQ=255)
             ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTg                                                                                   <  1:1083928/100‑1 (MQ=255)
             ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATCGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTg                                                                                   <  2:2129387/100‑1 (MQ=255)
              tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGc                                                                                  <  1:1849962/100‑1 (MQ=255)
              tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGc                                                                                  <  2:2028129/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCa                                                                                 <  1:2311859/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCa                                                                                 >  1:2079780/1‑100 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCa                                                                                 <  1:1493605/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCa                                                                                 <  1:1377391/100‑1 (MQ=255)
                  tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGc                                                                              >  1:1803998/1‑100 (MQ=255)
                    tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTg                                                                            <  2:2057916/100‑1 (MQ=255)
                       cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCa                                                                         <  2:702695/100‑1 (MQ=255)
                            gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGcgcg                                                                    <  2:1492243/100‑1 (MQ=255)
                            gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGcgcg                                                                    >  2:1794458/1‑100 (MQ=255)
                            gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑GGGACTTGCAGGCTGGCAGcgcg                                                                    >  1:608378/1‑100 (MQ=255)
                              aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgat                                                                  <  1:1887535/100‑1 (MQ=255)
                              aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgat                                                                  <  2:2273670/100‑1 (MQ=255)
                               aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatg                                                                 >  2:2382296/1‑100 (MQ=255)
                               aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑GGGACTTGCAGGCTGGCAGCGCgatg                                                                 >  1:908269/1‑100 (MQ=255)
                                tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                                >  2:259915/1‑100 (MQ=255)
                                tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                                >  2:722558/1‑100 (MQ=255)
                                 cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatgat                                                               >  2:572225/1‑100 (MQ=255)
                                  ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATc                                                              >  2:1524726/1‑100 (MQ=255)
                                     tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                           <  1:2620154/100‑1 (MQ=255)
                                      ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                          >  2:1887103/1‑100 (MQ=255)
                                      ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                          >  1:486133/1‑100 (MQ=255)
                                       tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                         <  1:2070986/100‑1 (MQ=255)
                                         aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                       <  2:196872/100‑1 (MQ=255)
                                         aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                       >  2:2934323/1‑100 (MQ=255)
                                          gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                      <  1:1153161/100‑1 (MQ=255)
                                           ccccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACGTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                     <  2:1619810/100‑1 (MQ=255)
                                            cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                    <  2:1021809/100‑1 (MQ=255)
                                             ccACGAGCATGGTGGATGTGGCGCACCTTTGGTGTAGAGAAGGTGTCAATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                     <  2:983731/98‑1 (MQ=255)
                                              cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                  >  2:2839673/1‑100 (MQ=255)
                                              cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                  <  1:921709/100‑1 (MQ=255)
                                                  aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGTAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                              <  2:2177716/100‑1 (MQ=255)
                                                  aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                              <  2:1248841/100‑1 (MQ=255)
                                                   gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                                <  1:1225369/97‑1 (MQ=255)
                                                   gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                             >  1:1426896/1‑100 (MQ=255)
                                                     ctggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                           <  2:1732577/99‑1 (MQ=255)
                                                     aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGGTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTg                                            >  1:606913/1‑99 (MQ=255)
                                                     aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                           <  2:2745888/100‑1 (MQ=255)
                                                      tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                          >  2:91326/1‑100 (MQ=255)
                                                      tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                          >  2:2251433/1‑100 (MQ=255)
                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                         <  1:1129108/100‑1 (MQ=255)
                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                         >  2:1796866/1‑100 (MQ=255)
                                                       ggtggATGCTGCGCACCGTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAACGCGATGGTCTGCTGTTAGAAGAAGGTGCAg                                         >  1:1732592/1‑100 (MQ=255)
                                                        gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                        <  1:2382318/100‑1 (MQ=255)
                                                        gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                        >  1:2451463/1‑100 (MQ=255)
                                                          ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                      <  2:879558/100‑1 (MQ=255)
                                                          ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                      <  2:1797180/100‑1 (MQ=255)
                                                           gATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                     >  1:792830/1‑100 (MQ=255)
                                                           gATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                     >  1:2610673/1‑100 (MQ=255)
                                                            aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                    <  2:2873198/100‑1 (MQ=255)
                                                            aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                    <  1:2820925/100‑1 (MQ=255)
                                                             tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                   >  2:877078/1‑100 (MQ=255)
                                                               ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                 <  2:2079761/100‑1 (MQ=255)
                                                               ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                 >  2:1737320/1‑100 (MQ=255)
                                                                tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCACGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                                >  2:2825410/1‑100 (MQ=255)
                                                                tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGGGGATGGACTTGCAGGCTGGTAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                                >  2:2125434/1‑100 (MQ=255)
                                                                actgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                                <  1:1394242/97‑1 (MQ=255)
                                                                 gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                               <  2:1718153/100‑1 (MQ=255)
                                                                  cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                              >  1:2402510/1‑100 (MQ=255)
                                                                    cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTGGAAGAAGGTGCAGTAGAGCTGCCGGa                            >  2:1525172/1‑100 (MQ=255)
                                                                     aCCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGaa                           <  1:774541/100‑1 (MQ=255)
                                                                         ttGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCACGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                       <  1:2825435/100‑1 (MQ=255)
                                                                           ggTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTGGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                     <  2:523495/100‑1 (MQ=255)
                                                                            gtgtAGAGAAAGTGTCGATTCT‑GGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                    >  2:2895844/1‑99 (MQ=255)
                                                                             tgtAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                   >  2:1399255/1‑100 (MQ=255)
                                                                                 gagaAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc               >  2:2753338/1‑100 (MQ=255)
                                                                                 gagaAAGTGTCGAGTCTGGGGGGG‑TGGACTTGCAGGCTAGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGCGCTGCCGGAAGGagag                     >  2:1873480/1‑94 (MQ=255)
                                                                                  agaAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg              <  1:722565/100‑1 (MQ=255)
                                                                                  agaAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg              >  1:747004/1‑100 (MQ=255)
                                                                                   gaAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc             >  1:1613296/1‑100 (MQ=255)
                                                                                     aaGTGTCGATCCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCGGGTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgcc            <  2:606906/99‑1 (MQ=255)
                                                                                       gtgtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg         >  1:1156884/1‑100 (MQ=255)
                                                                                       gtgtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg         <  1:2757652/100‑1 (MQ=255)
                                                                                         gtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCTGGAAGGAGAGTTTAACGCCGCGtt       >  2:1644478/1‑100 (MQ=255)
                                                                                         gtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt       >  1:1438885/1‑100 (MQ=255)
                                                                                         gtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt       >  1:2097981/1‑100 (MQ=255)
                                                                                         gtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt       >  1:1334537/1‑100 (MQ=255)
                                                                                          tCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt      >  1:2689884/1‑100 (MQ=255)
                                                                                          tCGATTCTGGGGGGG‑TGGACTTGCAGGCGGGCAGGGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt      <  1:487552/100‑1 (MQ=255)
                                                                                           cGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa     >  1:1445359/1‑100 (MQ=255)
                                                                                           cGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa     >  1:2534441/1‑100 (MQ=255)
                                                                                           cGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa     >  1:1656593/1‑100 (MQ=255)
                                                                                            gATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa    <  2:908261/100‑1 (MQ=255)
                                                                                            gATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa    >  2:2280147/1‑100 (MQ=255)
                                                                                             aTTCT‑GGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATc  >  1:2681129/1‑100 (MQ=255)
                                                                                                  |                                                                                                
GATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATC  >  NC_000913/2652752‑2652944

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