Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,654,144 T→A 15.5% intergenic (+805/+13) sseA → / ← sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,654,1440TA15.5% 508.9 / 73.7 198intergenic (+805/+13)sseA/sseB3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base T (92/75);  new base A (2/29);  total (94/104)
Fisher's exact test for biased strand distribution p-value = 2.15e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 205 total aligned reads displayed.

GCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTG  >  NC_000913/2654045‑2654243
                                                                                                   |                                                                                                   
acTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCAGGAGACAATTTACTGGTAGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCttt                                                                                                     >  2:214077/2‑100 (MQ=255)
 tttattaccgaacatattgcgccattctatgaacgtcgctggggtggttttttgcgtgactttaaacagaatcggataatctaaaatATGCCGCCCTTTg                                                                                                    >  1:1239084‑M1/88‑100 (MQ=255)
  ttattaccgaacatattgcgccattctatgaacgtcgctggggtggttttttgcgtgactttaaacagaatcggataatctaaaatATGCCGCCCTTTgg                                                                                                   <  1:267897‑M1/14‑1 (MQ=255)
  tGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTgg                                                                                                   >  1:3703513/1‑100 (MQ=255)
   gTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTggg                                                                                                  >  2:653/1‑100 (MQ=255)
     tAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcg                                                                                                >  1:1311568/1‑100 (MQ=255)
     tAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcg                                                                                                >  2:2322777/1‑100 (MQ=255)
     tAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcg                                                                                                >  2:546559/1‑100 (MQ=255)
     tAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTAggg                                                                                                  >  1:2598403/1‑98 (MQ=255)
      aaCCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcgg                                                                                               >  2:1510227/1‑100 (MQ=255)
       aCCGAGGTAAATAATAACCGGAGTCTCTCCGGCGATAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcggc                                                                                              <  1:2115166/100‑1 (MQ=255)
       aCCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcggc                                                                                              >  2:3470555/1‑100 (MQ=255)
       aCCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGAAAGCCCGAATGCTGCCCAAAGggcggc                                                                                              <  2:1878639/100‑1 (MQ=255)
        ccgaacatattgcgccattctatgaacgtcgctggggtggttttttgcgtgactttaaacagaatcggataatctaaaatATGCCGCCCTTTGGGCGGCa                                                                                             >  1:3567198‑M1/81‑100 (MQ=255)
        ccGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCGGCCCTTTGGGCGGCa                                                                                             >  1:3873631/1‑100 (MQ=255)
        ccGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcggc                                                                                              >  2:1658163/1‑99 (MQ=255)
        ccGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCa                                                                                             <  2:277932/100‑1 (MQ=255)
        ccGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCa                                                                                             >  1:3027484/1‑100 (MQ=255)
        ccGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCa                                                                                             <  1:2341062/100‑1 (MQ=255)
         cgaacatattgcgccattctatgaacgtcgctggggtggttttttgcgtgactttaaacagaatcggataatctaaaatATGCCGCCCTTTGGGCGGCa                                                                                             >  1:1791632‑M1/80‑99 (MQ=255)
         cGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCat                                                                                            <  1:804840/100‑1 (MQ=255)
         cGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCat                                                                                            <  2:2693665/100‑1 (MQ=255)
            ggTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattcg                                                                                         <  2:3660415/100‑4 (MQ=255)
            ggTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTAAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCata                                                                                           >  2:3261063/1‑98 (MQ=255)
              tAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcggagatcgga                                                                                       >  1:716964/1‑92 (MQ=255)
              tAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATAtttt                                                                                       >  2:2559143/1‑100 (MQ=255)
              tAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATAtttt                                                                                       <  2:1111505/100‑1 (MQ=255)
              tAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattcggg                                                                                       <  2:1065226/100‑6 (MQ=38)
              tAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattcg                                                                                         <  1:2366078/98‑4 (MQ=255)
                aataatAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAg                                                                                     <  1:3230943/100‑1 (MQ=255)
                aataatAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAg                                                                                     <  2:1026940/100‑1 (MQ=255)
                aataatAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAg                                                                                     >  1:3093314/1‑100 (MQ=255)
                aataatAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAg                                                                                     <  2:3890273/100‑1 (MQ=255)
                     taaCCGGAGTCTCTCCGGCGATAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                 <  1:1601564/99‑1 (MQ=255)
                     taaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGAGCGGCATATTTTAGATTAt                                                                                >  1:1691815/1‑100 (MQ=255)
                      aaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAATCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                 <  2:1286747/98‑1 (MQ=255)
                        ccgtattctttcCGGTGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCg                                                                             <  1:1394402/94‑1 (MQ=255)
                          ggAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAATCTTCTGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGAt                                                                           <  1:1032962/100‑1 (MQ=255)
                            ttccgatctggCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTc                                                                         <  2:3869816/91‑1 (MQ=255)
                                  tcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTa                                                                   <  1:3595503/100‑1 (MQ=255)
                                  tcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTa                                                                   <  2:3612658/100‑1 (MQ=255)
                                    cGGCGACAATTTACTGGTGGATAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTaaa                                                                 >  2:2397874/1‑100 (MQ=255)
                                     ggCGACAATTTACTGGTGGTTAACAATCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTAGCCGATTCTGTTTAAAg                                                                <  1:1305300/100‑1 (MQ=255)
                                       cGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTc                                                              >  1:2807589/1‑100 (MQ=255)
                                          cAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACg                                                           <  1:3552007/100‑1 (MQ=255)
                                           aaTTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGc                                                          >  1:973575/1‑100 (MQ=255)
                                                  tggtggTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAac                                                   <  2:2512550/100‑1 (MQ=255)
                                                   ggtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAacc                                                  <  1:2431180/100‑1 (MQ=255)
                                                   ggtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAac                                                   >  1:3816148/1‑99 (MQ=255)
                                                     tggtTAACAACCTTCAGAGTAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAaccac                                                <  2:1908434/100‑1 (MQ=255)
                                                     tggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTAACGCAAAAAaccac                                                >  1:1386991/1‑100 (MQ=255)
                                                        ttAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAaccacc                                               >  1:2866050/1‑98 (MQ=255)
                                                        ttAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAaccacc                                               >  2:1483329/1‑98 (MQ=255)
                                                        ttAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCAcccc                                             >  1:203218/1‑100 (MQ=255)
                                                         tAACAACCTTCAGAGAAGCAAGTAAGCCCGAATGCCGCCCTTTGGACGGCATATTTTAGATTATACAATTCTGTTTAAAGTCACGCAAAAAACCACCCCa                                            >  2:3319667/1‑100 (MQ=255)
                                                          aacaacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAACACACCGCa                                            >  2:1460738/1‑99 (MQ=255)
                                                          aacaacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAg                                           >  1:3807690/1‑100 (MQ=255)
                                                           acaacCTTCATAGCAGCAAGTAAGCCCGAATTCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGc                                          <  1:1159268/100‑1 (MQ=255)
                                                           acaacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGc                                          >  1:2460713/1‑100 (MQ=255)
                                                            caacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCg                                         >  2:462150/1‑100 (MQ=255)
                                                              acCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGAc                                       >  2:1353320/1‑100 (MQ=255)
                                                               ccTTCAGAGCAGCAAGTAAGCCCGAATGCCGTCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGAc                                       <  2:2856812/99‑1 (MQ=255)
                                                                  tCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGAGCGGCATATTTTAGATCATCCTATTATGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGtt                                    <  2:257820/99‑1 (MQ=255)
                                                                    agagCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCAt                                 <  1:2484377/100‑1 (MQ=255)
                                                                     gagCAGCAAGTAATCCCGAATGCCGCCTTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAATTCACGCAAAATTCCACCCAAGCGACGTTCATa                                <  2:446723/100‑1 (MQ=255)
                                                                     gagCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATa                                >  1:343350/1‑100 (MQ=255)
                                                                        cagcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGaa                             <  1:509189/100‑1 (MQ=255)
                                                                        cagcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGaa                             >  1:3857189/1‑100 (MQ=255)
                                                                        cagcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGaa                             >  2:2908196/1‑100 (MQ=255)
                                                                         agcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAAt                            <  2:826326/100‑1 (MQ=255)
                                                                         agcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAAt                            >  1:2881427/1‑100 (MQ=255)
                                                                         agcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCAAGCAAAAAACCACCCAAGCGACGTTCATAGAAt                            >  1:2240531/1‑100 (MQ=255)
                                                                          gcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCTACGTTCATAGAATg                           <  2:1561997/100‑1 (MQ=255)
                                                                          gcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATg                           >  1:3138610/1‑100 (MQ=255)
                                                                            aaGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGgc                         >  2:2665933/1‑100 (MQ=255)
                                                                             aGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGgcg                        <  2:2522551/100‑1 (MQ=255)
                                                                              gTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACATTCATAGaa                             >  2:1191178/1‑94 (MQ=255)
                                                                               atctaaaatATGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCa                      <  1:1568732‑M1/91‑1 (MQ=255)
                                                                                 ctaaaatATGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCGGCGACGTTCATAGAATGGCGCAat                    <  2:2936766‑M1/93‑1 (MQ=255)
                                                                                 aGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAat                    >  1:1921211/1‑100 (MQ=255)
                                                                                  taaaatATGCCGCCCAAAGGGCGGCATATTTTAGATTGTCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAata                   <  2:2079906‑M1/94‑1 (MQ=255)
                                                                                  gCCCGAATGTCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAata                   <  1:1894730/100‑1 (MQ=255)
                                                                                  gCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAata                   <  1:2193292/100‑1 (MQ=255)
                                                                                  gCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAata                   >  2:3167673/1‑100 (MQ=255)
                                                                                   cccGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAatat                  >  1:2915937/1‑100 (MQ=255)
                                                                                    ccGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCATGCAAAAAACCACCCCAGCGGCGTTCATAGAATGGCGCAATATg                 <  2:720557/100‑1 (MQ=255)
                                                                                     cGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTAACGCAAAAAACAACCCAAACGACGTTCATAGAATGGCGCAATATGt                >  2:796422/1‑100 (MQ=255)
                                                                                       aaTGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTc              <  2:2814895/100‑1 (MQ=255)
                                                                                        ATGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCg             <  1:3389256‑M1/100‑1 (MQ=255)
                                                                                         tgctgctttttGGGCGGTATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCg             <  1:845806/91‑1 (MQ=255)
                                                                                          GCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGt           >  2:2460690‑M1/1‑100 (MQ=255)
                                                                                          GCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGt           >  2:978133‑M1/1‑100 (MQ=255)
                                                                                           CCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGta          <  1:2734592‑M1/100‑1 (MQ=255)
                                                                                           CCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGta          <  2:1741279‑M1/100‑1 (MQ=255)
                                                                                            CGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACTCCAGCGACGTTCATAGAATTGCTCAATATGTTCGGta          <  2:869456‑M1/99‑1 (MQ=255)
                                                                                              CCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGtaata       <  1:2968151‑M1/100‑1 (MQ=255)
                                                                                               CCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCAGCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGtaataa      <  2:2754279‑M1/100‑1 (MQ=255)
                                                                                                  ttGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGt   <  2:2348941/100‑1 (MQ=255)
                                                                                                   tGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTg  <  1:3161026/100‑1 (MQ=255)
                                                                                                   tGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTg  <  2:1132421/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTG  >  NC_000913/2654045‑2654243

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: