Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,367 G→A 17.0% intergenic (+28/+790) 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,653,3670GA17.0% 313.8 / 48.1 129intergenic (+28/+790)sseA/sseB3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base G (51/56);  new base A (16/6);  total (67/62)
Fisher's exact test for biased strand distribution p-value = 3.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 142 total aligned reads displayed.

GTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAC  >  NC_000913/2653268‑2653465
                                                                                                   |                                                                                                  
gTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAAccgccg                                                                                                    <  2:102496/100‑1 (MQ=255)
gTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAAccgccg                                                                                                    <  2:2585403/100‑1 (MQ=255)
 tGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGa                                                                                                   >  1:2050485/1‑100 (MQ=255)
   aaaCTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTCTACAGGCAATAAAAAACCGCCAAAt                                                                                                 >  1:169445/1‑100 (MQ=255)
    aaCTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAAtt                                                                                                >  1:2173619/1‑100 (MQ=255)
       tGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTgg                                                                                             <  2:126688/100‑1 (MQ=255)
       tGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTgg                                                                                             <  2:1230604/100‑1 (MQ=255)
        gTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGc                                                                                            <  1:714471/100‑1 (MQ=255)
         tACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCg                                                                                           >  2:2617020/1‑100 (MQ=255)
          acgacgGCGCCTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGGGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                          >  1:1240536/1‑100 (MQ=255)
          acgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                          >  2:1460841/1‑100 (MQ=255)
          acgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                          <  1:854535/100‑1 (MQ=255)
          acgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                          <  1:853872/100‑1 (MQ=255)
           cgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                          <  2:1020170/99‑1 (MQ=255)
           cgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                          >  1:2029468/1‑99 (MQ=255)
               ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGttttt                                                                                     <  1:2608827/100‑1 (MQ=255)
               ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGttttt                                                                                     <  2:2605960/100‑1 (MQ=255)
               ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGttttt                                                                                     <  2:2016961/100‑1 (MQ=255)
                  gcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTAt                                                                                  <  2:1716847/100‑1 (MQ=255)
                   cATGGGGTGAATGGGGCGCGCGGGCAGTTTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTa                                                                                   <  1:93610/98‑1 (MQ=38)
                   cATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAtt                                                                                 <  2:2641313/100‑1 (MQ=255)
                   cATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAtt                                                                                 >  2:2025533/1‑100 (MQ=255)
                   cATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAtt                                                                                 <  1:276493/100‑1 (MQ=255)
                    aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                >  2:1809732/1‑100 (MQ=255)
                    aTGGAGTGAATGGGGAGCGCAGTCGGAGTTCCGGGTTGAACCAGTGAAATAGGTATTTTACTGGTCATAATAAACCGCCGAATTTGGCga                                                                                          <  1:1736023/90‑2 (MQ=255)
                     tGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGc                                                                               >  2:1902922/1‑100 (MQ=255)
                      ggAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCt                                                                              >  2:2566570/1‑100 (MQ=255)
                          tGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTc                                                                          >  2:2551942/1‑100 (MQ=255)
                           gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCt                                                                         >  2:454173/1‑100 (MQ=255)
                           gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCt                                                                         >  1:2258119/1‑100 (MQ=255)
                           gAATGGGGCGCGAGGGCAGATTTACCGGATGAGCCAGTGAAATAAGTATTTTAAAGGCAATAAAAAACCGCCGAATTTGGCGGATTTTTATAGCTAGTc                                                                          >  1:260057/1‑99 (MQ=38)
                             aTGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTgg                                                                       <  1:93491/100‑1 (MQ=255)
                                  gcgcgcGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTcgc                                                                  <  2:2183608/100‑1 (MQ=255)
                                    gcgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCgg                                                                <  2:2451671/100‑1 (MQ=255)
                                     ctcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGTAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCgg                                                                <  1:2025336/97‑1 (MQ=255)
                                     cgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTCTTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGCTCGCGGc                                                               >  1:2629776/1‑100 (MQ=255)
                                      gcgTGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGcc                                                              <  2:2459615/100‑1 (MQ=255)
                                       cgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCt                                                             >  2:1621259/1‑100 (MQ=255)
                                       cgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCt                                                             >  2:1339390/1‑100 (MQ=255)
                                        gggCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCtt                                                            >  1:353778/1‑100 (MQ=255)
                                          gCAGATTTACCGGTTGAGCCAGTGAAATCAGTATTTTACAGGCAATAAAAACCCGCCGAATTTCGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTc                                                          <  1:578769/100‑1 (MQ=255)
                                           cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTcc                                                         <  1:1902913/100‑1 (MQ=255)
                                           cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTcc                                                         <  2:1644686/100‑1 (MQ=255)
                                           cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTcc                                                         <  1:2566555/100‑1 (MQ=255)
                                           cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTc                                                          >  1:230926/1‑99 (MQ=255)
                                             gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                       >  1:2718354/1‑100 (MQ=255)
                                             gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                       <  1:2436795/100‑1 (MQ=255)
                                              aTTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagc                                                      <  2:1247987/100‑1 (MQ=255)
                                                ttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagcag                                                    >  2:368085/1‑100 (MQ=255)
                                                 tACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAgg                                                   <  2:46798/100‑1 (MQ=255)
                                                     ggTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGa                                               >  2:824642/1‑100 (MQ=255)
                                                     ggTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGa                                               >  2:842990/1‑100 (MQ=255)
                                                     ggTTGAGCCAGTGAAATAAGTATGTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTg                                                <  1:1809724/99‑1 (MQ=255)
                                                         gAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtg                                           >  2:2063004/1‑100 (MQ=255)
                                                         gAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtg                                           >  2:2657764/1‑100 (MQ=255)
                                                           gCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtg                                         >  2:2612009/1‑100 (MQ=255)
                                                           gCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtg                                         >  1:2059651/1‑100 (MQ=255)
                                                           gCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCCGCCTTTCCAGCAGGTTGACtt                                            >  1:618865/1‑97 (MQ=255)
                                                            ccAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTgtgt                                        <  2:1271953/100‑1 (MQ=255)
                                                            ccAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTgtgt                                        <  2:45895/100‑1 (MQ=255)
                                                            ccAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTgtgt                                        <  1:905621/100‑1 (MQ=255)
                                                             cAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTgtgt                                        <  1:961767/99‑1 (MQ=255)
                                                             cAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGtt                                       >  1:693529/1‑100 (MQ=255)
                                                             cAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGtt                                       >  2:276781/1‑100 (MQ=255)
                                                              aGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTa                                      >  1:2804696/1‑100 (MQ=255)
                                                                tGAAAGAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACa                                    <  2:900258/100‑1 (MQ=255)
                                                                  aaaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGTAGGTTGACTTGTGTTACATg                                  <  1:1781101/100‑1 (MQ=255)
                                                                  aaaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATg                                  <  1:107527/100‑1 (MQ=255)
                                                                     tAAGTATTTTACAGGCAATAAAAACCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGc                               <  1:1015952/100‑1 (MQ=255)
                                                                          accagactagCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGc                          >  2:629687/10‑100 (MQ=21)
                                                                           ttttACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCa                         <  2:1412394/100‑1 (MQ=255)
                                                                            tttACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAg                        <  1:362893/100‑1 (MQ=255)
                                                                              tACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAgatcggaagag                      >  2:2251761/1‑90 (MQ=255)
                                                                                 aGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCt                   <  2:187375/100‑1 (MQ=255)
                                                                                  ggCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCtt                  >  1:1367428/1‑100 (MQ=255)
                                                                                   gCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTc                 <  2:1703197/100‑1 (MQ=255)
                                                                                    cAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTca                >  2:2126318/1‑100 (MQ=255)
                                                                                    cAATAAAAAACCGCCAAATTCGGCGGGTTTTTATTGCTAGTCTGTTTCACGGCCTTTCCTTCAGGTTGACTGGTGTGACATGAGCAACGCAGGTGGTTca                >  2:860166/1‑100 (MQ=21)
                                                                                      aTAAAAAGCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTGcaca              <  1:786161/100‑1 (MQ=255)
                                                                                      aTAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTGCAGCAGTTTGAATTGTGTTTCATGAGCAACCCAGGTGCTTc                 >  1:2461819/1‑97 (MQ=25)
                                                                                      aTAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGCGCAACGCAGGTGCTTcaca              >  1:735964/1‑100 (MQ=255)
                                                                                      aTAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCCAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcaca              >  2:2204264/1‑100 (MQ=37)
                                                                                       tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg             >  1:2719245/1‑100 (MQ=255)
                                                                                       tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg             <  2:1291544/100‑1 (MQ=255)
                                                                                       tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg             <  2:751515/100‑1 (MQ=255)
                                                                                       tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGAAGGTTGACTTGTGTTACATGAGCCACGCAGGTGCTTCACAg             >  2:537594/1‑100 (MQ=255)
                                                                                       tAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg             >  2:2377230/1‑100 (MQ=255)
                                                                                       tAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg             >  1:1327472/1‑100 (MQ=255)
                                                                                       tAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg             >  2:2377091/1‑100 (MQ=255)
                                                                                        aaaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGc            >  2:2891880/1‑100 (MQ=255)
                                                                                        aaaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACCTGTGTTACATGAGCAACGCAGGTGCTTCACAGc            >  1:2329494/1‑100 (MQ=255)
                                                                                          aaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaa          >  2:2498288/1‑100 (MQ=255)
                                                                                           aaaCCGCCGAATTTGGCGGTTTTTTATTGCGAGTCTGGTTCCGGGCCTTTCCAACAGGTTGACTTTTGTTACATGAGCATCGCAGGTGCTTCACAACaaa         >  1:840444/1‑100 (MQ=255)
                                                                                             aCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGATACATGAGCAACGCAGGTGCTTCACAGCAAAAc       <  1:1177661/100‑1 (MQ=255)
                                                                                              ccgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc       <  1:397436/99‑1 (MQ=255)
                                                                                               cgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGAATTGTGTTACATGAGCAACGCAGGTGCTTCACAGCCAAACaa     >  1:1562715/1‑100 (MQ=255)
                                                                                                gcggAATTTGGCGGTTTTTTATTGCTAGTCAGGTTCGCGCCCTTCCCAGCAGGTTGCCTTGTGTTACATGAACAACGCAGGTGCTTCACAGCAAAAc       <  1:431563/94‑1 (MQ=255)
                                                                                                gccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAAt    <  2:2657499/100‑1 (MQ=255)
                                                                                                gccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAAt    >  1:2516990/1‑100 (MQ=255)
                                                                                                  cgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAc  <  2:693532/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAC  >  NC_000913/2653268‑2653465

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