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.9% 339.5 / 0.1 96coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (42/52);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 110 total aligned reads displayed.

AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAA  >  NC_000913/2652750‑2652942
                                                                                                     |                                                                                              
aGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCt                                                                                                 >  1:2504171/1‑100 (MQ=255)
aGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCt                                                                                                 >  2:1271345/1‑100 (MQ=255)
    tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggg                                                                                           <  1:3026689/100‑1 (MQ=255)
     aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTggggg                                                                                          <  1:965144/100‑1 (MQ=255)
      aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggg                                                                                         >  1:2308723/1‑100 (MQ=255)
      aGCACCTGATTATCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggg                                                                                         >  2:2576668/1‑100 (MQ=255)
       gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTGAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTggggggg                                                                                        <  2:614113/100‑2 (MQ=255)
        cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggt                                                                                       >  1:605433/1‑98 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggtgg                                                                                     >  1:648296/1‑96 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggtgg                                                                                     >  1:648744/1‑96 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggtgg                                                                                     >  1:1917038/1‑96 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggtgg                                                                                     <  1:1638244/100‑5 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggtgg                                                                                     >  2:1592195/1‑96 (MQ=255)
             gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACt                                                                                   >  1:2595720/1‑100 (MQ=255)
             gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACt                                                                                   <  2:478688/100‑1 (MQ=255)
               ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTg                                                                                 <  2:259977/100‑1 (MQ=255)
                 gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCa                                                                               <  1:478727/100‑1 (MQ=255)
                    tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGc                                                                            <  2:1906832/100‑1 (MQ=255)
                     aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCt                                                                           >  1:1501738/1‑100 (MQ=255)
                     aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCt                                                                           >  2:1510127/1‑100 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTg                                                                          <  1:1599789/100‑1 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTg                                                                          >  1:1781658/1‑100 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTgg                                                                         >  2:2509237/1‑100 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTgg                                                                         >  1:2044784/1‑100 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTgg                                                                         >  2:141983/1‑100 (MQ=255)
                       gACGAAGGTAATCGTGTCTCATCCCCACTAGCAGGGTGGATGCTGCGCACCTTTGGGGTAGAGAAAGTGGCGATTCTGGGGGGGTGGACTTGCAGGCTg                                                                          <  2:2588891/99‑1 (MQ=255)
                        aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGc                                                                        <  1:1666620/100‑1 (MQ=255)
                         cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCa                                                                       >  2:2875727/1‑100 (MQ=255)
                         cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCa                                                                       <  1:529645/100‑1 (MQ=255)
                          gAAGTTAATCTTTTCTCAGGCGCACGAGCATGGTGGATGCCGCGCACCTTCGGTGTAGAGAGAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTgg                                                                         <  2:214256/97‑1 (MQ=255)
                             ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAgcgc                                                                   >  2:1788548/1‑100 (MQ=255)
                              gTAATCTTTTCTCAGCCCCACGAGCATGGTTGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAgcgc                                                                   <  2:334701/99‑1 (MQ=255)
                               tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCga                                                                 >  1:2606489/1‑100 (MQ=255)
                               tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCga                                                                 >  1:1643355/1‑100 (MQ=255)
                                 aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatg                                                               <  2:2466185/100‑1 (MQ=255)
                                 aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatg                                                               <  2:2504137/100‑1 (MQ=255)
                                   cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatgat                                                             >  1:2410212/1‑100 (MQ=255)
                                   cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatgat                                                             >  1:1773340/1‑100 (MQ=255)
                                    ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                            <  2:1685863/100‑1 (MQ=255)
                                     tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATct                                                           <  1:2644464/100‑1 (MQ=255)
                                       tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                         <  2:311956/100‑1 (MQ=255)
                                       tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                         >  2:1119450/1‑100 (MQ=255)
                                           aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     <  1:1477252/100‑1 (MQ=255)
                                           aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     >  1:1866913/1‑100 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGTAGCGCGATGATCTGCTGTTa                                                    >  1:2798958/1‑100 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    <  1:880004/100‑1 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    >  1:2983878/1‑100 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    <  2:1643336/100‑1 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    <  2:1721256/100‑1 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    <  2:1723507/100‑1 (MQ=255)
                                               ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                 <  1:1732889/100‑1 (MQ=255)
                                               ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                 >  1:1584021/1‑100 (MQ=255)
                                                cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                >  1:856514/1‑100 (MQ=255)
                                                 aCGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaaga                                               <  2:866620/100‑1 (MQ=255)
                                                   gAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGaagaag                                             <  2:1501722/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            >  2:929265/1‑100 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:1193049/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:648736/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:648288/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:197298/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:568063/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            >  2:503723/1‑100 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  1:2457521/100‑1 (MQ=255)
                                                     gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCAGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                           >  2:1700158/1‑100 (MQ=255)
                                                       aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                         <  2:2595683/100‑1 (MQ=255)
                                                       aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                         >  1:2345001/1‑100 (MQ=255)
                                                        tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                        <  1:1495486/100‑1 (MQ=255)
                                                        tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                        >  2:134029/1‑100 (MQ=255)
                                                              aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                  >  2:2767506/1‑100 (MQ=255)
                                                               tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                 <  2:2161878/100‑1 (MQ=255)
                                                                 ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                               <  1:1271357/100‑1 (MQ=255)
                                                                 ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                               <  1:2551739/100‑1 (MQ=255)
                                                                   gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                             <  1:1668738/100‑1 (MQ=255)
                                                                    cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                            <  2:221418/100‑1 (MQ=255)
                                                                      cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                          >  2:1176718/1‑100 (MQ=255)
                                                                      cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                          <  2:2980940/100‑1 (MQ=255)
                                                                        ccTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                        <  2:173874/100‑1 (MQ=255)
                                                                         cTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                       <  1:2509271/100‑1 (MQ=255)
                                                                         cTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGACCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                       >  1:1219219/1‑100 (MQ=255)
                                                                          tttGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGCTGCAGTAGAGCTGCCGGAAGga                      <  2:1773320/100‑1 (MQ=255)
                                                                           ttGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                     <  2:856504/100‑1 (MQ=255)
                                                                             ggTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                   <  1:1196891/100‑1 (MQ=255)
                                                                              gtgtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                  <  1:2161093/100‑1 (MQ=255)
                                                                              gtgtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                  >  1:289668/1‑100 (MQ=255)
                                                                               tgtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                 >  2:1617354/1‑100 (MQ=255)
                                                                                gtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                >  2:1278852/1‑100 (MQ=255)
                                                                                gtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                <  1:1510143/100‑1 (MQ=255)
                                                                                 tAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa               <  2:605425/100‑1 (MQ=255)
                                                                                  agagAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGGAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa              >  2:3005130/1‑100 (MQ=255)
                                                                                   gagaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc             <  1:296378/100‑1 (MQ=255)
                                                                                   gagaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc             >  2:2460564/1‑100 (MQ=255)
                                                                                   gagaAAGTGTCGATTCTGTGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc             <  2:2766217/100‑1 (MQ=255)
                                                                                     gaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc           <  1:760870/100‑1 (MQ=255)
                                                                                        aGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccgc        >  1:2143136/1‑100 (MQ=255)
                                                                                          tgtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCGGCTGTTAGAAGAAGGTGCAGGAGAGCTGCCGGAAGGAGAGTTTAACgcggcg       >  2:1710961/1‑99 (MQ=255)
                                                                                           gtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt     >  2:283376/1‑100 (MQ=255)
                                                                                           gtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt     >  1:424094/1‑100 (MQ=255)
                                                                                             cGATTCTGGGGGGGGGGACTTGCAGGCTGGCAGCGCGCTGATCTGCTGTTAGAAGAAAGTGCAGTACAGCTGCCGGAAGGAGAGTTTAACTCCGCGTTTa   >  1:2512646/1‑100 (MQ=255)
                                                                                              gATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa  >  1:158488/1‑100 (MQ=255)
                                                                                              gATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa  >  2:1359284/1‑100 (MQ=255)
                                                                                                     |                                                                                              
AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAA  >  NC_000913/2652750‑2652942

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: