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

GATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATC  >  NC_000913/2652752‑2652944
                                                                                                  |                                                                                                
gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgg                                                                                                <  2:253696/100‑1 (MQ=255)
gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgg                                                                                                <  2:2881456/100‑1 (MQ=255)
 aTAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTggg                                                                                               <  2:1123012/100‑1 (MQ=255)
    aGCACCTGATTGTCTGTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggg                                                                                            <  1:2990099/100‑1 (MQ=255)
    aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggg                                                                                            >  2:627644/1‑100 (MQ=255)
     gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTggggggg                                                                                           >  1:2939415/1‑99 (MQ=255)
     gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTggggggg                                                                                           >  2:1849946/1‑99 (MQ=255)
        ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggtgg                                                                                        >  2:1917190/1‑96 (MQ=255)
         cTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGAc                                                                                      <  1:1786001/100‑1 (MQ=255)
           gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACt                                                                                     <  1:1220376/100‑1 (MQ=255)
            aTTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGATCACGTGTCGATTCTgggggggtggcgtt                                                                                    >  2:511298/1‑92 (MQ=255)
             ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTg                                                                                   <  1:1083928/100‑1 (MQ=255)
             ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATCGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTg                                                                                   <  2:2129387/100‑1 (MQ=255)
              tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGc                                                                                  <  1:1849962/100‑1 (MQ=255)
              tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGc                                                                                  <  2:2028129/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCa                                                                                 <  1:2311859/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCa                                                                                 >  1:2079780/1‑100 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCa                                                                                 <  1:1493605/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCa                                                                                 <  1:1377391/100‑1 (MQ=255)
                  tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGc                                                                              >  1:1803998/1‑100 (MQ=255)
                    tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTg                                                                            <  2:2057916/100‑1 (MQ=255)
                       cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCa                                                                         <  2:702695/100‑1 (MQ=255)
                            gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGcgcg                                                                    <  2:1492243/100‑1 (MQ=255)
                            gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGcgcg                                                                    >  2:1794458/1‑100 (MQ=255)
                            gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGGGGACTTGCAGGCTGGCAGcgcg                                                                    >  1:608378/1‑100 (MQ=255)
                              aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgat                                                                  <  1:1887535/100‑1 (MQ=255)
                              aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgat                                                                  <  2:2273670/100‑1 (MQ=255)
                               aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatg                                                                 >  2:2382296/1‑100 (MQ=255)
                               aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGGGGACTTGCAGGCTGGCAGCGCgatg                                                                 >  1:908269/1‑100 (MQ=255)
                                tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatga                                                                >  2:259915/1‑100 (MQ=255)
                                tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatga                                                                >  2:722558/1‑100 (MQ=255)
                                 cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatgat                                                               >  2:572225/1‑100 (MQ=255)
                                  ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                              >  2:1524726/1‑100 (MQ=255)
                                     tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                           <  1:2620154/100‑1 (MQ=255)
                                      ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                          >  2:1887103/1‑100 (MQ=255)
                                      ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                          >  1:486133/1‑100 (MQ=255)
                                       tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                         <  1:2070986/100‑1 (MQ=255)
                                         aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                       <  2:196872/100‑1 (MQ=255)
                                         aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                       >  2:2934323/1‑100 (MQ=255)
                                          gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                      <  1:1153161/100‑1 (MQ=255)
                                           ccccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACGTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                     <  2:1619810/100‑1 (MQ=255)
                                            cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                    <  2:1021809/100‑1 (MQ=255)
                                             ccACGAGCATGGTGGATGTGGCGCACCTTTGGTGTAGAGAAGGTGTCAATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                     <  2:983731/98‑1 (MQ=255)
                                              cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                  >  2:2839673/1‑100 (MQ=255)
                                              cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                  <  1:921709/100‑1 (MQ=255)
                                                  aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGTAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                              <  2:2177716/100‑1 (MQ=255)
                                                  aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                              <  2:1248841/100‑1 (MQ=255)
                                                   gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                                <  1:1225369/97‑1 (MQ=255)
                                                   gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                             >  1:1426896/1‑100 (MQ=255)
                                                     ctggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                           <  2:1732577/99‑1 (MQ=255)
                                                     aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGGTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTg                                            >  1:606913/1‑99 (MQ=255)
                                                     aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                           <  2:2745888/100‑1 (MQ=255)
                                                      tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                          >  2:91326/1‑100 (MQ=255)
                                                      tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                          >  2:2251433/1‑100 (MQ=255)
                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                         <  1:1129108/100‑1 (MQ=255)
                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                         >  2:1796866/1‑100 (MQ=255)
                                                       ggtggATGCTGCGCACCGTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAACGCGATGGTCTGCTGTTAGAAGAAGGTGCAg                                         >  1:1732592/1‑100 (MQ=255)
                                                        gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                        <  1:2382318/100‑1 (MQ=255)
                                                        gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                        >  1:2451463/1‑100 (MQ=255)
                                                          ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                      <  2:879558/100‑1 (MQ=255)
                                                          ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                      <  2:1797180/100‑1 (MQ=255)
                                                           gATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                     >  1:792830/1‑100 (MQ=255)
                                                           gATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                     >  1:2610673/1‑100 (MQ=255)
                                                            aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                    <  2:2873198/100‑1 (MQ=255)
                                                            aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                    <  1:2820925/100‑1 (MQ=255)
                                                             tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                   >  2:877078/1‑100 (MQ=255)
                                                               ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                 <  2:2079761/100‑1 (MQ=255)
                                                               ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                 >  2:1737320/1‑100 (MQ=255)
                                                                tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCACGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                                >  2:2825410/1‑100 (MQ=255)
                                                                tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGATGGACTTGCAGGCTGGTAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                                >  2:2125434/1‑100 (MQ=255)
                                                                actgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                                <  1:1394242/97‑1 (MQ=255)
                                                                 gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                               <  2:1718153/100‑1 (MQ=255)
                                                                  cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                              >  1:2402510/1‑100 (MQ=255)
                                                                    cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTGGAAGAAGGTGCAGTAGAGCTGCCGGa                            >  2:1525172/1‑100 (MQ=255)
                                                                     aCCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGaa                           <  1:774541/100‑1 (MQ=255)
                                                                         ttGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCACGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                       <  1:2825435/100‑1 (MQ=255)
                                                                           ggTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTGGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                     <  2:523495/100‑1 (MQ=255)
                                                                            gtgtAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                    >  2:2895844/1‑99 (MQ=255)
                                                                             tgtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                   >  2:1399255/1‑100 (MQ=255)
                                                                                 gagaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc               >  2:2753338/1‑100 (MQ=255)
                                                                                 gagaAAGTGTCGAGTCTGGGGGGGTGGACTTGCAGGCTAGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGCGCTGCCGGAAGGagag                     >  2:1873480/1‑94 (MQ=255)
                                                                                  agaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg              <  1:722565/100‑1 (MQ=255)
                                                                                  agaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg              >  1:747004/1‑100 (MQ=255)
                                                                                   gaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc             >  1:1613296/1‑100 (MQ=255)
                                                                                     aaGTGTCGATCCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCGGGTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgcc            <  2:606906/99‑1 (MQ=255)
                                                                                       gtgtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg         >  1:1156884/1‑100 (MQ=255)
                                                                                       gtgtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg         <  1:2757652/100‑1 (MQ=255)
                                                                                         gtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCTGGAAGGAGAGTTTAACGCCGCGtt       >  2:1644478/1‑100 (MQ=255)
                                                                                         gtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt       >  1:1438885/1‑100 (MQ=255)
                                                                                         gtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt       >  1:2097981/1‑100 (MQ=255)
                                                                                         gtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt       >  1:1334537/1‑100 (MQ=255)
                                                                                          tCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt      >  1:2689884/1‑100 (MQ=255)
                                                                                          tCGATTCTGGGGGGGTGGACTTGCAGGCGGGCAGGGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt      <  1:487552/100‑1 (MQ=255)
                                                                                           cGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa     >  1:1445359/1‑100 (MQ=255)
                                                                                           cGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa     >  1:2534441/1‑100 (MQ=255)
                                                                                           cGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa     >  1:1656593/1‑100 (MQ=255)
                                                                                            gATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa    <  2:908261/100‑1 (MQ=255)
                                                                                            gATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa    >  2:2280147/1‑100 (MQ=255)
                                                                                             aTTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATc  >  1:2681129/1‑100 (MQ=255)
                                                                                                  |                                                                                                
GATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATC  >  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: