Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 305,344 A→C 19.5% G27G (GGT→GGG ydeP_1 ← Protein YdeP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440305,3440AC19.5% 114.7 / 8.2 64G27G (GGT→GGGydeP_1Protein YdeP
Reads supporting (aligned to +/- strand):  ref base A (21/29);  new base C (2/12);  total (23/41)
Fisher's exact test for biased strand distribution p-value = 6.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGGGAGTCCCCCCAGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTGGG  >  KT2440/305210‑305477
                                                                                                                                      |                                                                                                                                     
tCCGGGGAGTCCCCCCAGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATATTCATCAGCCCATTATCGCTGCCGACCCAGGCCGTGGCCACGCTACGCAGTGCGCCCCaa                                                                                                                                       >  5:422311/1‑135 (MQ=255)
        ggccccccAGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGGTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCcggcg                                                                                                                              <  7:290548/134‑1 (MQ=255)
          ccccccAGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCcggcgg                                                                                                                             <  4:177484/135‑1 (MQ=255)
           cccccAGGCGCAGCCGGGGCAGGCGAAACCGCCGTTCTGGGTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCggg                                                                                                                            <  6:48981/135‑1 (MQ=255)
            ccccAGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCgggg                                                                                                                           >  4:247135/1‑135 (MQ=255)
            ccccAGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGc                                                                                                                          >  7:291854/1‑136 (MQ=255)
            ccccAGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTAGGCGCATTTTCGCTGCCGCCCCAGGCCGTGGCCACGCTACGCAGTGCGCCCCACCCGCCGCCgggg                                                                                                                           >  3:255045/1‑135 (MQ=25)
               cAGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGGGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGt                                                                                                                       <  8:30176/136‑1 (MQ=255)
               cAGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGt                                                                                                                       <  7:330147/136‑1 (MQ=255)
               cAGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGGTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGt                                                                                                                       <  7:390031/136‑1 (MQ=255)
                aGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGt                                                                                                                       <  2:272021/135‑1 (MQ=255)
                aGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGt                                                                                                                       <  4:406054/135‑1 (MQ=255)
                aGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTGCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGt                                                                                                                       <  1:47740/135‑1 (MQ=255)
                aGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGGTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGt                                                                                                                       <  2:71910/135‑1 (MQ=255)
                  gcgcAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTgg                                                                                                                     <  4:391072/135‑1 (MQ=255)
                    gcAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCcggcgg                                                                                                                             >  7:404739/1‑125 (MQ=255)
                    gcAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCcggcgg                                                                                                                             <  8:404739/125‑1 (MQ=255)
                      aGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCACCCGCCGGCAGGGCCGTGGTAggg                                                                                                                >  8:219954/1‑136 (MQ=255)
                       gCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAggg                                                                                                                <  6:81565/135‑1 (MQ=255)
                        ccGGGGCAGTCGAAACCGCCGTTCGGGTTGGGCTTGAGCAGGGCGCGAATGTTCGTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCCACCGCCGGCGGGGCCGTgggcgg                                                                                                                 >  2:136928/1‑129 (MQ=38)
                         cGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGcg                                                                                                              >  5:150220/1‑134 (MQ=255)
                            ggCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTa                                                                                                          <  7:251279/136‑1 (MQ=255)
                              cAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAg                                                                                                         <  1:254579/135‑1 (MQ=255)
                                   gAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGt                                                                                                                       >  2:93908/1‑116 (MQ=255)
                                   gAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGt                                                                                                                       <  1:93908/116‑1 (MQ=255)
                                   gAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTAAGCGCATTAACCCTTCCGACCCCAGCCGTGGCCCCGCTACGCAGTGCGCCCCCACCGCCGGCGGGGCCCTGGTAGGGGTGGTTGGGAggg                                                                                                   >  7:182017/1‑136 (MQ=25)
                                   gAAACCGCCGTTCTGGTTGGGCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAgg                                                                                                    <  6:65169/135‑1 (MQ=255)
                                     acccgccgTTCTGGGTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAggggg                                                                                                 <  7:394530/134‑1 (MQ=255)
                                     aaCCGCCGTTCTGGTTGGTCTTGAGCTGGGCGCGAATGTTCTTCAGCGCATTATCGCGGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCGGTGGTAGGGCTTGTAGCGAgggg                                                                                                  <  4:407071/135‑1 (MQ=255)
                                         gccgTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGGAGGGCGTGTAGCGAgggggggg                                                                                              >  5:33995/1‑134 (MQ=255)
                                              tCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCATCGCATTATCGCTGCCGACCCAGGCCTTGGCCACCCTACGCAGTGCGCCCCAACCCCCGGCCGGGCCGTGGTAGGGGATGTAGCGAGGGGGGGAGGCTg                                                                                         >  6:36542/1‑135 (MQ=255)
                                                   ttggtCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTggg                                                                                       <  4:89116/132‑1 (MQ=255)
                                                         ttGAGCCGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTCCGCAGTGCGCCCCAACCGCGGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTcc                                                                              <  6:359732/135‑1 (MQ=255)
                                                                    gcgAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGttgttg                                                                   >  5:172639/1‑135 (MQ=255)
                                                                    gcgAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGttgttg                                                                   >  6:137683/1‑135 (MQ=255)
                                                                    gccatggTTCTTCAGCTCTTTTTGGCTGCGCACCCTGGCCTGGGCCTCTCTACGCCTGGCGCCCCACGCGCCGGCGGGGGGGTGGTAGTGATTGTCGCTTGGGGGGGAGGATGGGGTGTTGTCCGGAAGTTGTTgg                                                                  <  8:182017/130‑1 (MQ=255)
                                                                           ttcttcAGCGCATTCTCGCGGCCGCCCCTGGCGTGGGCCACGCTCCGCAGTGCGCCCCACCCGCCGGGGGGGCGGGGGTTGGGCGGGTTGCGAGGGGGGGAGGCTGGGGTGGTGTCCGGAAGTTGTTGGTACGa                                                             <  6:371904/134‑1 (MQ=255)
                                                                              ttcAGCGCATTAACGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCCCCCCAACCCCCCGCGGGGCCGGGGGAGGGCGGGTTGAGCGGGGGGGAGGCGGGGGGGTTGTGGGGTAGTTTTTTGTACGAGGTc                                                         >  2:191106/1‑135 (MQ=255)
                                                                                cAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCAc                                                       >  3:118154/1‑135 (MQ=255)
                                                                                      aTTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTg                                                  >  1:88502/1‑134 (MQ=255)
                                                                                      aTTATCCCTGCGCACCACGGCGTTCGCCACGCACCGCCGGGCCCCCCTCCCGGCGGGGTGGCTGTGGAAGGCGTTGTAGCGAGGGGGGGAGGCTTGGGTGTTGCCCGGATGTTGTTGGAACGAGGCCACGACTGt                                                 <  4:14039/135‑1 (MQ=255)
                                                                                        tATCGCTGCCGACCCAGGCCTGGGCCACGCTACGCAGGGCGCCCCACCCGCCGGCGGGGCCGTGGTAGGGGTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTa                                               <  3:435825/135‑1 (MQ=255)
                                                                                            gcggCCGACCCCGGCCTGGGCCCCGCTCCGCGGTGCCCCCCCCCCGGGGGGGGGGGGGGGTTTGGGGTGGTAGCGCGGGGGGGAGGCGGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGa                                            <  1:212222/131‑1 (MQ=255)
                                                                                             cTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGGTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACt                                          <  4:377780/135‑1 (MQ=255)
                                                                                                ccGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTgg                                                                  <  4:378590/108‑1 (MQ=255)
                                                                                                ccGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTgg                                                                  >  3:378590/1‑108 (MQ=255)
                                                                                                ccGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTccgc                                      >  7:326159/1‑136 (MQ=255)
                                                                                                  cccccGGGCGTGGGCCCCGCTACGCCGGGCGCCCCGCCGGCCGGGGGGGCGGTGGTTGGGCTTGTAGGGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTccgccg                                    <  7:392391/134‑1 (MQ=255)
                                                                                                   cccccGGCCTGGGCCCCGCCCCGCGGGGCGCCCCCACCGCGGGGGGGGCCGGGGGAGGGCGGGCAGCGCGGGGGGGCGGGGGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGGCACGACTGTTACGACTccgccg                                    <  6:172639/134‑1 (MQ=255)
                                                                                                   aCCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCGGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTccgccg                                    <  4:184483/135‑1 (MQ=255)
                                                                                                   aCCCAGGCCTGGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGGTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTccgccg                                    <  2:295709/135‑1 (MQ=255)
                                                                                                           cTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTGCGGAAGTTGGTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTg                            >  5:189546/1‑135 (MQ=255)
                                                                                                             tGGCCACGCACCGCTGGGCGCCCCCCCCGCCGGGGGGGCGGGGGGAGGGCTGGGCGCGCGGGGGGGAGGGGGGGGGGGTGTCCGGAAGGTGTTGTTACGAGGTCACGACTGGTACGACTCCGCCGCTGGGCTGTa                          <  5:248085/135‑1 (MQ=255)
                                                                                                                 cACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGGGTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTTGTACGACGTCACCCCTGTTACGACACCGCCGCGGGGCGGTAGccc                      >  5:301121/1‑135 (MQ=255)
                                                                                                                 cACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGAccc                     >  8:120416/1‑136 (MQ=255)
                                                                                                                 cACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGAcc                      >  5:171833/1‑135 (MQ=255)
                                                                                                                      gaCGCCGTGCGCCCCAACCCGCGGGGGGGGCGTGGTTGGGCGTGTAGCGCGGGGGGGCGGCGGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGGCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGt                <  8:103597/135‑1 (MQ=255)
                                                                                                                        cgctgggcgcCCCCCCCGCCGGGGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTg               <  1:381967/129‑1 (MQ=255)
                                                                                                                        cGCAGGGCGCCCCAACCGCCGGCGGGGCCGTGGGTGGGCTTGTCGCGTGGGGGGGTGGCGGGGGTGTTGTCCGGACGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGt                <  8:86514/134‑1 (MQ=255)
                                                                                                                         gccgggcgcCCCACCCGCCGGGGGGGCGGGGGTCGGGCTGGCCGCGAGGGGGGGAGGGGGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGt                <  3:247135/128‑1 (MQ=255)
                                                                                                                         gccgTGCGCCCCCCCCGCCGGGGGGGCGGTGGTAGGGGTGGTCGCGAGGGGGGGAGGGTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGt                <  5:256325/130‑1 (MQ=255)
                                                                                                                           aGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTgcgc            >  3:3330/1‑135 (MQ=255)
                                                                                                                             ggcgcCCCCCCCCCGGGCGGGGCGGGGGGGGGGGTGGTCGGGCGGGGGGGCGGGGGGGGGGGTGGCCGGAAGTTGTTGGTCCGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCtg          <  6:422311/134‑1 (MQ=255)
                                                                                                                              ccgcCCCACCCGCGGGGGGGGCGGGGGTTGGGGTTTTAGGGTGGGGGGGAGGCTGGGGTGTTGTCCGGATGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCtgt         <  5:85349/134‑1 (MQ=255)
                                                                                                                                    cAACCGCCGGGGGGGCGGTGGTAGGGCTTGTCGCGAGGGGGGGTGGCTGGGGTGTTGGCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTggg  <  8:275298/136‑1 (MQ=255)
                                                                                                                                     aaCCGCCGGGGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTggg  <  3:147473/135‑1 (MQ=255)
                                                                                                                                     aaCCGCCGGCGGGGCCGTGGTAGGGGTTGTAGCGTGGGGGGGAGGCTGGGGTGGTGTCCGGGAGTTGGGGGTAGGAGGTCCCGGCTGTTACGCCTCCGCCGCTGGGCTGTTGCCCCGCCGGGCGCTGTGCGTggg  >  2:43511/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                     
TCCGGGGAGTCCCCCCAGGCGCAGCCGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTGGG  >  KT2440/305210‑305477

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