Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,290,634 G→A 6.3% T185T (ACC→ACT)  dgcN ← Diguanylate cyclase DgcN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,290,6340GA 6.3% 175.1 / 3.0 63T185T (ACC→ACT) dgcNDiguanylate cyclase DgcN
Reads supporting (aligned to +/- strand):  ref base G (20/39);  new base A (2/2);  total (22/41)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

GTTCATCCAGCAGGGCATTGAAGTCGTTACCCAGCTCGTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGCAGGAACAGCAACAGGCTCCTGCCCTGGCCTGCCAGTT  >  KT2440/5290501‑5290761
                                                                                                                                     |                                                                                                                               
gTTCATCCAGCAGGGCATTGAAGTCGTTACCCAGCTCGGTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGcgggcggg                                                                                                                                <  6:28698/135‑1 (MQ=255)
   cATCCAGCAGGGCATTGAAGTCGTTACCCAGCTCGTTCAGCTCGGCAATCGGCGCGTCTGGTACGCGCCGGTCGAAGCAGCGCACGCGGGCGGCGGCATGGGTGCCACGGGCCCCGCCACGCAGccggcggggg                                                                                                                              >  1:383154/1‑132 (MQ=255)
        agcagGGCATTGAAGTCGTTACCCAGCTCGTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGCCACTGGCCAGGCCACGCAGCGGGCGGGTGATATCg                                                                                                                        >  3:70446/1‑135 (MQ=255)
        agcagGGCATTGAAGTCGTTACCCAGCTCGTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACCCGGGCCAGGCCCCGCAGCGGGCGGGGGaaatag                                                                                                                        >  5:236948/1‑133 (MQ=255)
        agcagGGCATTGAAGTCGTTACCCAGCTCGTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCg                                                                                                                        >  8:214046/1‑135 (MQ=255)
             ggCATTGAAGTCGTTACCCAGCTCGTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAg                                                                                                                   >  4:361980/1‑135 (MQ=255)
                aTTGAAGTCGTTACCCAGCTCGTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAg                                                                                                                >  3:163209/1‑135 (MQ=255)
                  tgcaGTCGTTACCCCGCTCGAGCAGCTCGGCAACCGGCGCTGCTGGTCCGCGCCGGTCGACGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCGGAGCAGgc                                                                                                              <  2:218642/132‑1 (MQ=255)
                   gacggcGTTCCCCACCGCGGTCCGCGCGGCAAGCGGCCCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGgc                                                                                                              <  2:179842/129‑1 (MQ=255)
                       tCGTTACCCAGCTCGTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGcgccgc                                                                                                         <  4:332546/135‑1 (MQ=255)
                       gcGTTCCCCAGCTCGTTCCCCGCGCCATGGGGCGCTGCGGGTCCGCGCCGGGCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACGGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGcgccgc                                                                                                         <  1:200775/134‑1 (MQ=255)
                          ttACCCAGCTCGTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCACGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCcgcg                                                                                                        >  1:142732/1‑133 (MQ=255)
                            aCCCAGCTCGTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGCCACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCcgcg                                                                                                        >  6:200043/1‑131 (MQ=255)
                             cccAGCTCGTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAg                                                                                                    <  3:51900/134‑1 (MQ=255)
                                 gCTCGTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACt                                                                                               <  5:90073/135‑1 (MQ=255)
                                  cTCGTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTg                                                                                              <  4:315364/135‑1 (MQ=255)
                                     gTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTgg                                                                                             >  8:25172/1‑133 (MQ=255)
                                     gTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGc                                                                                           >  8:17028/1‑135 (MQ=255)
                                      ttCAGCTCGGCAATCGGCGCTTCTGGTACCCGCCGGTCGAAGCTGCGCACGCGGCGGGCGGCCTGGGCGACACTGGCCAGGCCACCCCGCGGGCGGGGGATATCGCTGAGCCGGCGCCGCCACAGATAATTGGCa                                                                                          >  6:241764/1‑135 (MQ=255)
                                         aGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGcc                                                                                       <  7:25172/135‑1 (MQ=255)
                                         aGCTCGGCAATCGGCGCGTCGGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGcc                                                                                       <  8:142280/135‑1 (MQ=255)
                                         aGCGCGGCACTCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGcc                                                                                       <  1:74267/135‑1 (MQ=255)
                                              ggcacgcGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCGCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGgc                                                                                  <  1:10394/129‑1 (MQ=255)
                                              ggCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGgc                                                                                  >  2:122051/1‑135 (MQ=255)
                                              ggCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGCGCAGCCAGGgc                                                                                  >  2:106967/1‑135 (MQ=255)
                                                     ggCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAgg                                                                           <  4:22716/135‑1 (MQ=255)
                                                        gcTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCTGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGAcc                                                                        <  8:23879/135‑1 (MQ=255)
                                                        gcTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGAcc                                                                        <  1:301989/135‑1 (MQ=255)
                                                        gcTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGAcc                                                                        <  5:129223/135‑1 (MQ=255)
                                                          ttCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGt                                                                      <  1:307236/135‑1 (MQ=255)
                                                             tGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGCCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCa                                                                   >  7:351352/1‑135 (MQ=255)
                                                                tccgcgCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAg                                                                  <  3:231066/131‑1 (MQ=255)
                                                                tACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCa                                                                   <  4:354776/132‑1 (MQ=255)
                                                                tACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCa                                                                   <  3:193062/132‑1 (MQ=255)
                                                                tACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAg                                                                  <  3:284633/133‑1 (MQ=255)
                                                                tACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAg                                                                  <  8:324166/133‑1 (MQ=255)
                                                                    cgcCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAgg                                                            <  7:51835/135‑1 (MQ=255)
                                                                    cgcCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAgg                                                            <  3:249212/135‑1 (MQ=255)
                                                                    cgcCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAgg                                                            <  4:243015/135‑1 (MQ=255)
                                                                       cGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGGGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATa                                                         >  6:60269/1‑135 (MQ=255)
                                                                       cGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGAGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATa                                                         >  8:357059/1‑135 (MQ=255)
                                                                       cGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGAGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATa                                                         >  8:361736/1‑135 (MQ=255)
                                                                                 tGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGc                                               <  3:190592/135‑1 (MQ=255)
                                                                                 ggcgcTCGCGGCGGGCGGCAGGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATAACGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGACAAGc                                               <  6:126518/134‑1 (MQ=255)
                                                                                    gcTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCAc                                            <  6:236948/135‑1 (MQ=255)
                                                                                    ccaccaggcgggccGCATGTGCGACACTGCCCAGGCCAAGCAGCCGGCGGGTGATATCGCTGACCACGCGCCGCGACAGATACTGGGCAGCCAGGCCGCTCAGGACCGTTCAGAACAGGATACCGCCCAAGCCAc                                            >  2:10342/7‑135 (MQ=255)
                                                                                         cggcggGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGAGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGc                                       <  7:361736/135‑1 (MQ=255)
                                                                                         cggcggGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGAGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGc                                       <  7:357059/135‑1 (MQ=255)
                                                                                          ggcgggcgGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGCa                                      <  2:142732/135‑1 (MQ=255)
                                                                                           gcgggcggCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGCAg                                     <  2:72721/135‑1 (MQ=255)
                                                                                                 ggcATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGCAGGAACAg                               <  8:150495/135‑1 (MQ=255)
                                                                                                  gcATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGCAGGAACAGc                              <  3:307438/135‑1 (MQ=255)
                                                                                                   cATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGCAGGAACAGCa                             <  7:17028/135‑1 (MQ=255)
                                                                                                                  ccaggccaCGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGCAGGAACAGCAACAGGCTCCTGcc                <  1:106967/133‑1 (MQ=255)
                                                                                                                    aggccaCGCAGCGGGCGGGTGATATCGCTGAGCAGgc                                                                                                              >  7:269792/1‑37 (MQ=255)
                                                                                                                    aggccaCGCAGCGGGCGGGTGATATCGCTGAGCAGgc                                                                                                              <  8:269792/37‑1 (MQ=255)
                                                                                                                    aggccaCGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAg                                                                                                    >  7:39240/1‑47 (MQ=255)
                                                                                                                    aggccaCGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAg                                                                                                    <  8:39240/47‑1 (MQ=255)
                                                                                                                    aggccaCGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGcc                                                                                       >  1:325845/1‑60 (MQ=255)
                                                                                                                    aggccaCGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGcc                                                                                       <  2:325845/60‑1 (MQ=255)
                                                                                                                      gccaCGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGCAGGAACAGCAACAGGCTCCTGCCCTGGcc          <  4:163209/135‑1 (MQ=255)
                                                                                                                         aCGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGa                                                                             <  1:352638/65‑1 (MQ=255)
                                                                                                                         aCGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGa                                                                             >  2:352638/1‑65 (MQ=255)
                                                                                                                         aCGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGCAGGAACAGCAACAGGCTCCTGCCCTGGCCTg        >  6:173746/1‑134 (MQ=255)
                                                                                                                         aCGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGCAGGAACAGCAACAGGCTCCTGCCCTGGCCTGc       <  2:120579/135‑1 (MQ=255)
                                                                                                                              gcgggcggGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGa                                                                             <  6:353970/60‑1 (MQ=255)
                                                                                                                              gcgggcggGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGa                                                                             >  5:353970/1‑60 (MQ=255)
                                                                                                                              gcgggcggGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGCAGGAACAGCAACAGGCTCCTGCCCTGGCCTGCCAGtt  >  8:380538/1‑135 (MQ=255)
                                                                                                                              gcgggcggGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGCAGGAACAGCAACAGCCTCCTGCCCTCGCCTGCCAggt  >  6:10490/1‑133 (MQ=255)
                                                                                                                                    gggTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGa                                                                             >  3:53126/1‑54 (MQ=255)
                                                                                                                                    gggTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGa                                                                             <  4:53126/54‑1 (MQ=255)
                                                                                                                                     |                                                                                                                               
GTTCATCCAGCAGGGCATTGAAGTCGTTACCCAGCTCGTTCAGCTCGGCAATCGGCGCTTCTGGTACGCGCCGGTCGAAGCTGCGCTCGCGGCGGGCGGCATGGGCGACACTGGCCAGGCCACGCAGCGGGCGGGTGATATCGCTGAGCAGGCGCCGCGACAGATACTGGGCAGCCAGGGCGCTGAGGACCGTGCAGAACAGGATACCGGCCAAGCCACTGAGCAGGAACAGCAACAGGCTCCTGCCCTGGCCTGCCAGTT  >  KT2440/5290501‑5290761

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