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.G96.9% 230.3 / 1.3 64coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (27/35);  total (28/36)
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.96e-01

CGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAAT  >  NC_000913/2652704‑2652994
                                                                                                                                                   |                                                                                                                                                  
cgcCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggg                                                                                                                                                <  1:1118597/150‑1 (MQ=255)
cgcCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggg                                                                                                                                                <  2:516559/150‑1 (MQ=255)
cgcCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggg                                                                                                                                                <  2:1034309/150‑1 (MQ=255)
    cGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggtg                                                                                                                                           <  1:761261/151‑4 (MQ=255)
        aaCGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACt                                                                                                                                       <  2:674025/151‑1 (MQ=255)
                ccGTGGCGATGCGTGAATTAGGCGTTTACCAGGTTAAGCGCCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGTGCATGGTTGATGCTGCGCACCTTTGGTGTAAAGAAAGTGTCTATTCGGGGGGGGTGGACTTGCAGGc                                                                                                                                <  2:791802/150‑1 (MQ=255)
                ccGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGc                                                                                                                                <  1:1060244/150‑1 (MQ=255)
                         tGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGcgcg                                                                                                                      >  2:400345/1‑151 (MQ=255)
                         tGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGACGTTTCTGGGGGGAGGGACTTGCAGGCTGGCAGcgcg                                                                                                                      >  1:560510/1‑151 (MQ=255)
                             tGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGGCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTTGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatga                                                                                                                  <  2:499818/151‑1 (MQ=255)
                             tGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatga                                                                                                                  <  2:297533/151‑1 (MQ=255)
                               aaTTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatgat                                                                                                                 >  1:1079866/1‑150 (MQ=255)
                               aaTTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                                                                                <  2:605458/151‑1 (MQ=255)
                               aaTTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGAGGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                                                                                <  2:973168/151‑1 (MQ=255)
                                    ggCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctg                                                                                                              >  1:18929/1‑148 (MQ=255)
                                     gCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                                                                           <  1:1063452/150‑1 (MQ=255)
                                     gCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                                                                           <  2:970426/150‑1 (MQ=255)
                                         tAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                                                                      <  2:822322/151‑1 (MQ=255)
                                         tAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                                                                      <  1:501774/151‑1 (MQ=255)
                                          acccgggATAAGCACATGATTGTTTATGACGAAGGTAATAATTCCTCAGCCCCACGAGCATGGTGGATGCCGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGAACTGCTGTTagaa                                                                                                     <  2:62965/146‑1 (MQ=255)
                                          aaCCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                                                                     >  1:455133/1‑151 (MQ=255)
                                          aaCCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                                                                     >  1:512040/1‑151 (MQ=255)
                                          aaCCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                                                                     >  2:260661/1‑151 (MQ=255)
                                            ccAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaaga                                                                                                   <  2:786894/151‑1 (MQ=255)
                                            ccAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaaga                                                                                                   <  2:283000/151‑1 (MQ=255)
                                            ccAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaaga                                                                                                   >  1:759222/1‑151 (MQ=255)
                                                gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCTCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                                                                                  <  2:925733/148‑1 (MQ=255)
                                                gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAggg                                                                                               >  1:1241051/1‑150 (MQ=255)
                                                     gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTTCAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                                                                          >  2:1221371/1‑151 (MQ=255)
                                                      cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                                                                          <  2:177268/150‑1 (MQ=255)
                                                       aCCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                                                                        <  2:453967/151‑1 (MQ=255)
                                                          tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                                                                     <  2:316076/151‑1 (MQ=255)
                                                               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                                                                                >  1:340700/1‑151 (MQ=255)
                                                               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                                                                                >  1:1260028/1‑151 (MQ=255)
                                                               gTCTATGACGAAGGAAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                                                                                >  2:1051012/1‑151 (MQ=255)
                                                                      aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                                                                         <  2:1200406/151‑1 (MQ=255)
                                                                                tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc                                                               <  2:697198/151‑1 (MQ=255)
                                                                                tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc                                                               <  2:1251686/151‑1 (MQ=255)
                                                                                   tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccg                                                             >  1:744845/1‑150 (MQ=255)
                                                                                            cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATcc                                                   <  1:393727/151‑1 (MQ=255)
                                                                                            cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATcc                                                   <  2:454987/151‑1 (MQ=255)
                                                                                                    cATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGc                                             <  1:705941/149‑1 (MQ=255)
                                                                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCgtgg                                        >  1:232391/1‑151 (MQ=255)
                                                                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCgtgg                                        >  1:381501/1‑151 (MQ=255)
                                                                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCgtgg                                        >  1:265671/1‑151 (MQ=255)
                                                                                                            aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGaaa                                   >  1:678344/1‑151 (MQ=255)
                                                                                                               ctgcGCACCTTTGGTGTAGAGACAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGACGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGTGTTTAATCCTGCAGCAGTGGTGAAAGTa                                >  2:970530/1‑151 (MQ=255)
                                                                                                               ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTa                                >  2:423526/1‑151 (MQ=255)
                                                                                                               ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTa                                >  1:538469/1‑151 (MQ=255)
                                                                                                                  cgcACCTTTGGTGTAGAGAAATTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTTTTGTGAGAATAAGGTGTAGTAGAGCTGCAGGAAGGAGAGTTTAACTCCGCTTTTAATCCTGAAGCCGTTGTGaaattgac                              <  2:560281/150‑6 (MQ=255)
                                                                                                                  cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAAc                              <  2:124231/150‑1 (MQ=255)
                                                                                                                  cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAAc                              <  2:756606/150‑1 (MQ=255)
                                                                                                                      ccTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGGGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGGGAAAGTAACCGATg                         >  1:188053/1‑151 (MQ=255)
                                                                                                                           ggTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGtatta                    >  1:616130/1‑151 (MQ=255)
                                                                                                                            gtgtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGtattat                   >  1:339130/1‑151 (MQ=255)
                                                                                                                            gtgtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGtattat                   >  1:709218/1‑151 (MQ=255)
                                                                                                                             tgtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGtattat                   <  1:1073395/150‑1 (MQ=255)
                                                                                                                             tgtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGtattat                   <  1:573554/150‑1 (MQ=255)
                                                                                                                              gtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGACAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTg                 >  1:1114866/1‑151 (MQ=255)
                                                                                                                              gtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTGCCAGGTTGGCAGCGCGATGATTGGCTGTTAGAAGAAGGGGCAGTAAAGCTGCCGGAAGGTGAGTTTAAGGCCGGGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATCATTg                 >  2:742516/1‑151 (MQ=255)
                                                                                                                               tAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTg                 <  1:568761/150‑1 (MQ=255)
                                                                                                                                  agaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGGGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGc               <  1:290279/149‑1 (MQ=255)
                                                                                                                                  agaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGc               <  2:1000261/149‑1 (MQ=255)
                                                                                                                                  agaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCa              <  1:260721/150‑1 (MQ=255)
                                                                                                                                  agaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGACGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGc               <  2:648011/149‑1 (MQ=255)
                                                                                                                                           cGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGaaa    <  1:602475/151‑1 (MQ=255)
                                                                                                                                             aTTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATTTATTATTGGCAAGCCATGaaaa   <  2:682155/150‑1 (MQ=255)
                                                                                                                                             aTTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAAt  >  1:85503/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCCCGGAAACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAAT  >  NC_000913/2652704‑2652994

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: