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.

AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑G‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAA  >  NC_000913/2652750‑2652942
                                                                                                     |                                                                                              
aGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑t                                                                                                 >  1:2504171/1‑100 (MQ=255)
aGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑t                                                                                                 >  2:1271345/1‑100 (MQ=255)
    tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g‑ggg                                                                                           <  1:3026689/100‑1 (MQ=255)
     aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g‑gggg                                                                                          <  1:965144/100‑1 (MQ=255)
      aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g‑ggggg                                                                                         >  1:2308723/1‑100 (MQ=255)
      aGCACCTGATTATCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g‑ggggg                                                                                         >  2:2576668/1‑100 (MQ=255)
       gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTGAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g‑gggggg                                                                                        <  2:614113/100‑2 (MQ=255)
        cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g‑ggggggt                                                                                       >  1:605433/1‑98 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g‑ggggggtgg                                                                                     >  1:648296/1‑96 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g‑ggggggtgg                                                                                     >  1:648744/1‑96 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g‑ggggggtgg                                                                                     >  1:1917038/1‑96 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g‑ggggggtgg                                                                                     <  1:1638244/100‑5 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g‑ggggggtgg                                                                                     >  2:1592195/1‑96 (MQ=255)
             gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACt                                                                                   >  1:2595720/1‑100 (MQ=255)
             gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACt                                                                                   <  2:478688/100‑1 (MQ=255)
               ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTg                                                                                 <  2:259977/100‑1 (MQ=255)
                 gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCa                                                                               <  1:478727/100‑1 (MQ=255)
                    tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGc                                                                            <  2:1906832/100‑1 (MQ=255)
                     aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCt                                                                           >  1:1501738/1‑100 (MQ=255)
                     aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCt                                                                           >  2:1510127/1‑100 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTg                                                                          <  1:1599789/100‑1 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTg                                                                          >  1:1781658/1‑100 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTgg                                                                         >  2:2509237/1‑100 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTgg                                                                         >  1:2044784/1‑100 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTgg                                                                         >  2:141983/1‑100 (MQ=255)
                       gACGAAGGTAATCGTGTCTCATCCCCACTAGCAGGGTGGATGCTGCGCACCTTTGGGGTAGAGAAAGTGGCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTg                                                                          <  2:2588891/99‑1 (MQ=255)
                        aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGc                                                                        <  1:1666620/100‑1 (MQ=255)
                         cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCa                                                                       >  2:2875727/1‑100 (MQ=255)
                         cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCa                                                                       <  1:529645/100‑1 (MQ=255)
                          gAAGTTAATCTTTTCTCAGGCGCACGAGCATGGTGGATGCCGCGCACCTTCGGTGTAGAGAGAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTgg                                                                         <  2:214256/97‑1 (MQ=255)
                             ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAgcgc                                                                   >  2:1788548/1‑100 (MQ=255)
                              gTAATCTTTTCTCAGCCCCACGAGCATGGTTGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAgcgc                                                                   <  2:334701/99‑1 (MQ=255)
                               tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCga                                                                 >  1:2606489/1‑100 (MQ=255)
                               tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCga                                                                 >  1:1643355/1‑100 (MQ=255)
                                 aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCgatg                                                               <  2:2466185/100‑1 (MQ=255)
                                 aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCgatg                                                               <  2:2504137/100‑1 (MQ=255)
                                   cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCgatgat                                                             >  1:2410212/1‑100 (MQ=255)
                                   cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCgatgat                                                             >  1:1773340/1‑100 (MQ=255)
                                    ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                            <  2:1685863/100‑1 (MQ=255)
                                     tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATct                                                           <  1:2644464/100‑1 (MQ=255)
                                       tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                         <  2:311956/100‑1 (MQ=255)
                                       tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                         >  2:1119450/1‑100 (MQ=255)
                                           aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     <  1:1477252/100‑1 (MQ=255)
                                           aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     >  1:1866913/1‑100 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGTAGCGCGATGATCTGCTGTTa                                                    >  1:2798958/1‑100 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    <  1:880004/100‑1 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    >  1:2983878/1‑100 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    <  2:1643336/100‑1 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    <  2:1721256/100‑1 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    <  2:1723507/100‑1 (MQ=255)
                                               ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                 <  1:1732889/100‑1 (MQ=255)
                                               ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                 >  1:1584021/1‑100 (MQ=255)
                                                cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                >  1:856514/1‑100 (MQ=255)
                                                 aCGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaaga                                               <  2:866620/100‑1 (MQ=255)
                                                   gAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGaagaag                                             <  2:1501722/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            >  2:929265/1‑100 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:1193049/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:648736/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:648288/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:197298/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:568063/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            >  2:503723/1‑100 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  1:2457521/100‑1 (MQ=255)
                                                     gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCAG‑G‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                           >  2:1700158/1‑100 (MQ=255)
                                                       aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                         <  2:2595683/100‑1 (MQ=255)
                                                       aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                         >  1:2345001/1‑100 (MQ=255)
                                                        tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                        <  1:1495486/100‑1 (MQ=255)
                                                        tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                        >  2:134029/1‑100 (MQ=255)
                                                              aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                  >  2:2767506/1‑100 (MQ=255)
                                                               tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                 <  2:2161878/100‑1 (MQ=255)
                                                                 ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                               <  1:1271357/100‑1 (MQ=255)
                                                                 ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                               <  1:2551739/100‑1 (MQ=255)
                                                                   gcgcACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                             <  1:1668738/100‑1 (MQ=255)
                                                                    cgcACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                            <  2:221418/100‑1 (MQ=255)
                                                                      cACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                          >  2:1176718/1‑100 (MQ=255)
                                                                      cACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                          <  2:2980940/100‑1 (MQ=255)
                                                                        ccTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                        <  2:173874/100‑1 (MQ=255)
                                                                         cTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                       <  1:2509271/100‑1 (MQ=255)
                                                                         cTTTGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGACCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                       >  1:1219219/1‑100 (MQ=255)
                                                                          tttGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGCTGCAGTAGAGCTGCCGGAAGga                      <  2:1773320/100‑1 (MQ=255)
                                                                           ttGGTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                     <  2:856504/100‑1 (MQ=255)
                                                                             ggTGTAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                   <  1:1196891/100‑1 (MQ=255)
                                                                              gtgtAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                  <  1:2161093/100‑1 (MQ=255)
                                                                              gtgtAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                  >  1:289668/1‑100 (MQ=255)
                                                                               tgtAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                 >  2:1617354/1‑100 (MQ=255)
                                                                                gtAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                >  2:1278852/1‑100 (MQ=255)
                                                                                gtAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                <  1:1510143/100‑1 (MQ=255)
                                                                                 tAGAGAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa               <  2:605425/100‑1 (MQ=255)
                                                                                  agagAAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGGAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa              >  2:3005130/1‑100 (MQ=255)
                                                                                   gagaAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc             <  1:296378/100‑1 (MQ=255)
                                                                                   gagaAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc             >  2:2460564/1‑100 (MQ=255)
                                                                                   gagaAAGTGTCGATTC‑T‑GTGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc             <  2:2766217/100‑1 (MQ=255)
                                                                                     gaAAGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc           <  1:760870/100‑1 (MQ=255)
                                                                                        aGTGTCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccgc        >  1:2143136/1‑100 (MQ=255)
                                                                                          tgtCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCGGCTGTTAGAAGAAGGTGCAGGAGAGCTGCCGGAAGGAGAGTTTAACgcggcg       >  2:1710961/1‑99 (MQ=255)
                                                                                           gtCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt     >  2:283376/1‑100 (MQ=255)
                                                                                           gtCGATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt     >  1:424094/1‑100 (MQ=255)
                                                                                             cGATTC‑TGG‑GGGGGGGGACTTGCAGGCTGGCAGCGCGCTGATCTGCTGTTAGAAGAAAGTGCAGTACAGCTGCCGGAAGGAGAGTTTAACTCCGCGTTTa   >  1:2512646/1‑100 (MQ=255)
                                                                                              gATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa  >  1:158488/1‑100 (MQ=255)
                                                                                              gATTC‑TGG‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa  >  2:1359284/1‑100 (MQ=255)
                                                                                                     |                                                                                              
AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑G‑GGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAA  >  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: ‑