Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,652,855:1 (G)6→7 100% 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.G95.0% 334.2 / 9.7 99coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base G (52/41);  total (56/43)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 104 total aligned reads displayed.

GATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCT‑GGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652752‑2652943
                                                                                                   |                                                                                               
gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCT‑gg                                                                                               <  1:1161024/100‑1 (MQ=255)
 aTAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCT‑ggg                                                                                              <  2:1724536/100‑1 (MQ=255)
  tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCT‑gggg                                                                                             <  2:633489/100‑1 (MQ=255)
    aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCT‑gggg‑gg                                                                                          >  2:1324064/1‑100 (MQ=255)
     gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCT‑gggg‑ggg                                                                                         >  2:1034139/1‑99 (MQ=255)
        ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACATTTGGTGTAGAGAAAGAGTCGA‑TTCT‑GGGG‑GGTGGa                                                                                      >  2:2775047/1‑100 (MQ=255)
          tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGAc                                                                                     >  2:1022937/1‑100 (MQ=255)
          tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGAc                                                                                     <  1:98678/100‑1 (MQ=255)
             ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACATTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTg                                                                                  >  2:1698974/1‑100 (MQ=255)
              tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGc                                                                                 <  2:1122877/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCCTGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGTAAGTGTAGA‑TCTT‑GGGG‑GGCgggactt                                                                                  >  1:1921168/1‑92 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCa                                                                                <  2:612995/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCa                                                                                <  2:1560075/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCa                                                                                >  1:911767/1‑100 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCa                                                                                <  1:2400687/100‑1 (MQ=255)
                 cTATGACGACGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAgg                                                                              <  1:1154603/100‑1 (MQ=255)
                 cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAgg                                                                              >  2:2153394/1‑100 (MQ=255)
                 cTATCACTAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGATGCGCACCTTTAATGCAGAGAACGTGTCGA‑TTCT‑GGGGTGGTTGAATTGCAgg                                                                              >  1:2617579/1‑100 (MQ=255)
                  tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGc                                                                             <  2:1829984/100‑1 (MQ=255)
                   aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCt                                                                            <  1:2170066/100‑1 (MQ=255)
                   aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTGGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCt                                                                            <  1:1925109/100‑1 (MQ=255)
                     gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTgg                                                                          >  2:2857130/1‑100 (MQ=255)
                      aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTTGTGTAGAGCAAGTGTCGC‑TTCTGGGGG‑GGTGGACTTGCAGGATGGc                                                                         >  1:2354799/1‑100 (MQ=255)
                      aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGc                                                                         <  1:2692375/100‑1 (MQ=255)
                      aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGAGGACTTGCAGGCTGGc                                                                         >  1:1702118/1‑100 (MQ=255)
                          aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAgcg                                                                     >  2:1132614/1‑100 (MQ=255)
                           ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAgcgc                                                                    <  1:2446971/100‑1 (MQ=255)
                            gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGcgcg                                                                   <  1:2739492/100‑1 (MQ=255)
                            gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCCGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGcgcg                                                                   >  1:65657/1‑100 (MQ=255)
                             tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TGCTCGGGG‑GGTGGACTTGCAGGCTGGCAGCGCga                                                                  <  1:1921796/100‑1 (MQ=255)
                             tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGAGTGGAGTTTG‑GGGG‑GGTGGACTTGCAGGCTGGCAGCGCga                                                                  >  1:385520/1‑100 (MQ=255)
                              aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCgat                                                                 >  2:2210141/1‑100 (MQ=255)
                              aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCgat                                                                 <  2:796575/100‑1 (MQ=255)
                              aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCgat                                                                 <  2:797105/100‑1 (MQ=255)
                               aTCTTTTCTCAGCCCCACGTGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTTGACTTGCAGGCTGGCAGCGCgatg                                                                <  1:313647/100‑1 (MQ=255)
                               aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCgatg                                                                <  2:353946/100‑1 (MQ=255)
                               aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCgatg                                                                >  2:2208713/1‑100 (MQ=255)
                               aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCgatg                                                                >  1:2803193/1‑100 (MQ=255)
                                tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCgatga                                                               >  2:1697932/1‑100 (MQ=255)
                                 cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCgatgat                                                              >  1:897799/1‑100 (MQ=255)
                                   tttCTCAGCCCCATGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATct                                                            <  1:1385166/100‑1 (MQ=255)
                                   tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATct                                                            <  1:2857111/100‑1 (MQ=255)
                                   tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                             <  1:610543/99‑1 (MQ=255)
                                    ttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATctg                                                           >  2:1541067/1‑100 (MQ=255)
                                     tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                          >  2:1726029/1‑100 (MQ=255)
                                      ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                         <  2:1438231/100‑1 (MQ=255)
                                      ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                         >  2:2180783/1‑100 (MQ=255)
                                      ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                          >  2:1676123/1‑99 (MQ=255)
                                      ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGACctgct                                                         >  1:1930954/1‑100 (MQ=255)
                                      ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTAGCAGCGCGATGATctgct                                                         >  1:2237502/1‑100 (MQ=255)
                                      ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑GTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATcatt                                                          >  1:1274333/1‑96 (MQ=255)
                                        cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGt                                                       >  2:239794/1‑100 (MQ=255)
                                         aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  1:1812173/100‑1 (MQ=255)
                                         aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  2:2171529/100‑1 (MQ=255)
                                          gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     >  1:823697/1‑100 (MQ=255)
                                          gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     >  1:2299988/1‑100 (MQ=255)
                                            cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   <  2:1228918/100‑1 (MQ=255)
                                            cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   >  2:1333639/1‑100 (MQ=255)
                                            cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGAGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaaa                                                   >  2:673121/1‑98 (MQ=255)
                                             ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                  <  2:396902/100‑1 (MQ=255)
                                                   gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATGTGCTGTTAGAAGAAGGt                                            >  1:1355184/1‑100 (MQ=255)
                                                     aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                          <  1:2062823/100‑1 (MQ=255)
                                                      tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                         >  1:1230111/1‑100 (MQ=255)
                                                        gtggATGCTGCGCACCTTTGGTGTCGAGACAGCGTCGC‑TTATGGGGG‑GGTGGTATTGCAGGGTGGCAGAGCGATAATCTGCTGTTAAAAGAAGGTGCAgg                                       <  1:241784/100‑2 (MQ=255)
                                                        gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       >  1:1837062/1‑100 (MQ=255)
                                                          ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                     <  2:1460298/100‑1 (MQ=255)
                                                            aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                   >  1:519039/1‑100 (MQ=255)
                                                            aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                   >  1:2168973/1‑100 (MQ=255)
                                                             tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                  >  2:2177803/1‑100 (MQ=255)
                                                              gctgcGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCt                                 <  2:1294025/100‑1 (MQ=255)
                                                               ctgcGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                >  2:264366/1‑100 (MQ=255)
                                                               ctgcGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                >  2:1280535/1‑100 (MQ=255)
                                                                tgcGCACCTTTGGTGTAGAGAATGTGTCGC‑TTCT‑GGGG‑GGTGGACTTGCGGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                              <  1:2775029/100‑1 (MQ=255)
                                                                tgcGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                               >  1:2537614/1‑100 (MQ=255)
                                                                tgcGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                               >  1:2539578/1‑100 (MQ=255)
                                                                  cgcACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                             <  2:1806511/100‑1 (MQ=255)
                                                                  cgcACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                             >  2:1427854/1‑100 (MQ=255)
                                                                  cgcACCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                             >  2:1134426/1‑100 (MQ=255)
                                                                   gcACCTTTGATGTAGAGAAAGTGCCGA‑TTCTGGGGG‑GGTGGACTAGCAGGCGGTCAGCGCGATGATCCGCTGTTAGAAGAAGGTGCAGTAGAGCGGCCgg                            <  1:1953133/100‑1 (MQ=255)
                                                                     aCCTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTGGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           <  1:2247367/99‑1 (MQ=255)
                                                                      ccTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                         <  1:2380944/100‑1 (MQ=255)
                                                                       cTTTGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        >  1:1544450/1‑100 (MQ=255)
                                                                        tttGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGga                       <  2:1290989/100‑1 (MQ=255)
                                                                         ttGGTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                      <  1:1017072/100‑1 (MQ=255)
                                                                           ggTGTAGAGAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    >  1:2749357/1‑100 (MQ=255)
                                                                                agagAAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  2:2232006/1‑100 (MQ=255)
                                                                                 gagaAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGTTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              >  2:8056/1‑100 (MQ=255)
                                                                                 gagaAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              >  2:746283/1‑100 (MQ=255)
                                                                                 gagaAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              <  1:2337981/100‑1 (MQ=255)
                                                                                  agaAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGTCTGGCAGGGGGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg             <  1:657195/100‑1 (MQ=255)
                                                                                   gaAAGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGCCTGACAGCGCGATGATCTGCTGTTAGAAGGAGGTGCAGTAGAGCTGCCGGAAAGAGAGTTTAAcgc            >  1:1337361/1‑100 (MQ=255)
                                                                                      aGTGTCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccgc         >  2:2575564/1‑100 (MQ=255)
                                                                                       gtgtCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg        <  1:2153381/100‑1 (MQ=255)
                                                                                         gtCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  1:1805471/1‑100 (MQ=255)
                                                                                         gtCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  2:4036/1‑100 (MQ=255)
                                                                                         gtCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      <  2:560874/100‑1 (MQ=255)
                                                                                         gtCGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      <  2:1800849/100‑1 (MQ=255)
                                                                                           cGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    >  2:375610/1‑100 (MQ=255)
                                                                                           cGA‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    >  1:386194/1‑100 (MQ=255)
                                                                                             a‑TTCTGGGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt  <  2:65658/100‑1 (MQ=255)
                                                                                                   |                                                                                               
GATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGA‑TTCT‑GGGG‑GGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652752‑2652943

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: ‑