Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 2,921,458 C→T 5.9% C556C (TGC→TGT rep_1 → ATP‑dependent DNA helicase Rep

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,921,4580CT 5.9% 213.7 / 3.7 68C556C (TGC→TGTrep_1ATP‑dependent DNA helicase Rep
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base T (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

TGTTTGGCGAGAACAACCCAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGCCCAAGTATTTCATG  >  KT2440/2921331‑2921590
                                                                                                                               |                                                                                                                                    
tGTTTGGCGAGAACAACCCAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGTACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCt                                                                                                                                >  1:374542/1‑134 (MQ=255)
tGTTTGGCGAGAACAACCCAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTg                                                                                                                               >  7:341868/1‑135 (MQ=255)
 gTTTGGCGAGAACAACCCAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGa                                                                                                                              >  7:34878/1‑135 (MQ=255)
     ggCGAGAACAACCCAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGCTCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATa                                                                                                                          >  1:427371/1‑135 (MQ=255)
     ggCGAGAACAACCCAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATa                                                                                                                          >  1:21655/1‑135 (MQ=255)
     ggCGAGAACAACCCAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATa                                                                                                                          >  1:437462/1‑135 (MQ=255)
         agaACAACCCAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGgc                                                                                                                      >  4:201265/1‑135 (MQ=255)
          gaACAACCCAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGgcg                                                                                                                     <  5:434599/135‑1 (MQ=255)
            acaacCCAGAAACCCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGcgccg                                                                                                                  <  8:331350/136‑1 (MQ=255)
              accccAGAAACCCTGCACTTACTGCCAGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGCTCGCCGAGGTCAAGGACGTTTCGACGCCGATCGTCAGCGGGTTTCGTGCCTGCCTCGCTGAAATACGGcgccgc                                                                                                                 >  6:400223/3‑135 (MQ=255)
              aacCCAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCcgcg                                                                                                                <  7:229388/136‑1 (MQ=255)
               acCCAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCcgcg                                                                                                                <  3:215144/135‑1 (MQ=255)
                cccAGAAAACCTGCACTTACTGCCCGAGCCGTTGCTGCGGCCGCGGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGcggc                                                                                                              <  8:403105/136‑1 (MQ=255)
                 ccAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGcggc                                                                                                              <  3:6405/135‑1 (MQ=255)
                 ccAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGGGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGcggc                                                                                                              <  2:153020/135‑1 (MQ=255)
                   aGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGCCGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGcggccc                                                                                                            >  1:431184/1‑135 (MQ=255)
                   aGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGTCTCGCTGAAATACGGCGCCGcggccc                                                                                                            >  1:79309/1‑135 (MQ=255)
                   aGAAAACCTGCACTTACTCCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGTCTCGCTGAAATACGGCGCCGcggccc                                                                                                            >  1:79306/1‑135 (MQ=255)
                       aaCCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGcggcccggcc                                                                                                        <  2:10777/135‑1 (MQ=255)
                        aCCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCa                                                                                                       >  6:121897/1‑135 (MQ=255)
                                     gCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGgcgc                                                                                         >  8:71338/1‑136 (MQ=255)
                                       ccGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCCCACCGGGCGcc                                                                                        >  4:292323/1‑135 (MQ=255)
                                         gATCCGTTGCTGCGGCCGCTGGCGCACCTGTTGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCAt                                                                                      <  2:242848/135‑1 (MQ=255)
                                             cGTTGCTGCGGCCGCTGGCGCACCTGTTGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCa                                                                                 <  7:230389/136‑1 (MQ=255)
                                             cGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCCCACCGGGCGCCATGccc                                                                                  >  2:376200/1‑135 (MQ=255)
                                             cGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGc                                                                                    >  5:138555/1‑133 (MQ=255)
                                               ttGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGc                                                                               <  7:5238/136‑1 (MQ=255)
                                                  ctgcGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCa                                                                              >  2:381916/1‑134 (MQ=255)
                                                   tgcGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGcggccc                                                                                                            >  4:394404/1‑103 (MQ=255)
                                                   tgcGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGcggccc                                                                                                            <  3:394404/103‑1 (MQ=255)
                                                   tgcGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCAtt                                                                            >  3:209591/1‑135 (MQ=255)
                                                     cGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATtcg                                                                          <  7:287324/135‑1 (MQ=255)
                                                        ccGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATtcgtcg                                                                       >  5:101478/1‑135 (MQ=255)
                                                        ccGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATtcgtcg                                                                       >  4:188977/1‑135 (MQ=255)
                                                                            cgacgaGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTgacga                                                   >  4:295292/1‑135 (MQ=255)
                                                                               cgaGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTgg                                               >  7:317643/1‑136 (MQ=255)
                                                                                gaGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGc                                              <  7:6974/136‑1 (MQ=255)
                                                                                gaGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCATGCCTCGCTTAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGc                                              <  7:66493/136‑1 (MQ=255)
                                                                                 aGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAATTACGGCGCCGCGGCCCTTCCATGCACACCGGGCGACATGCCCTGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGc                                              <  6:342946/135‑1 (MQ=255)
                                                                                 aGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCTCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGc                                              <  3:435904/135‑1 (MQ=255)
                                                                                 aGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGc                                              <  6:173736/135‑1 (MQ=255)
                                                                                 aGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGc                                              <  1:381916/135‑1 (MQ=255)
                                                                                 aGTTCCAGGACGTGTCGCCGCAGATCGTCAGATGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGc                                              <  6:138555/135‑1 (MQ=255)
                                                                                   ttCCCGGACGGGTCGCCGCAGATCGTCAGCTGGTGGCGGTACTGCCTCGCTTAAATACGGCGACGCGGCCCGGCCATGAACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTTCGTGGGTGACGACTGGc                                              <  5:400223/133‑1 (MQ=255)
                                                                                   ttCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCaa                                            <  5:121897/135‑1 (MQ=255)
                                                                                    tCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATtcg                                                                          >  8:96894/1‑104 (MQ=255)
                                                                                    tCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATtcg                                                                          <  7:96894/104‑1 (MQ=255)
                                                                                                cgccgcAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATtcgtcg                                                                       <  4:411375/95‑1 (MQ=255)
                                                                                                cgccgcAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATtcgtcg                                                                       >  3:411375/1‑95 (MQ=255)
                                                                                                cgccgcAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCt                               >  2:246969/1‑135 (MQ=255)
                                                                                                cgccgcAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCt                               >  3:107163/1‑135 (MQ=255)
                                                                                                 gccgcAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTg                              >  3:129098/1‑135 (MQ=255)
                                                                                                  ccgcAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTgg                             <  2:21655/135‑1 (MQ=255)
                                                                                                       gATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTgg                        <  1:246969/135‑1 (MQ=255)
                                                                                                         tCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCa                      <  5:406027/135‑1 (MQ=255)
                                                                                                         tCGGCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCa                      <  2:374542/135‑1 (MQ=255)
                                                                                                           gTCAGCGGGTTGCGGGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGGGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGt                    <  1:232758/135‑1 (MQ=255)
                                                                                                            tCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGt                    <  2:92280/134‑1 (MQ=255)
                                                                                                            tCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACAACGGGCGCCAGGCCCAGCATTCGTCGTTGCTGTGCGTGGGTTACGGCTGGAAATCGCTCTACGGCTGGCGTGGCAGTTc                  >  8:433464/1‑136 (MQ=255)
                                                                                                              aGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGc                <  7:64814/136‑1 (MQ=255)
                                                                                                               gCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCg                 <  6:390724/134‑1 (MQ=255)
                                                                                                               gCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGcc               <  8:34878/136‑1 (MQ=255)
                                                                                                                cTGGTTGCGTGCCTGTCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGcc               <  2:79309/135‑1 (MQ=255)
                                                                                                                cTGGTTGCGTGCCTGTCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGcc               <  2:79306/135‑1 (MQ=255)
                                                                                                                   gTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGCCCaa            <  8:317643/135‑1 (MQ=255)
                                                                                                                     tgggTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGCCCAAGt          >  3:178684/4‑135 (MQ=255)
                                                                                                                     tGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGCCCAAGt          >  6:143538/1‑135 (MQ=255)
                                                                                                                        gTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGCCCAAGTAtt       <  3:283683/135‑1 (MQ=255)
                                                                                                                          gcctgcctCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGCCCAAGTATTTc     <  2:437462/135‑1 (MQ=255)
                                                                                                                          gcctgcctCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGCCCAAGTATTTc     <  6:429107/135‑1 (MQ=255)
                                                                                                                          gcctgcctCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGCCCAAGTATTTc     >  1:39592/1‑135 (MQ=255)
                                                                                                                          gcctgcctCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGCCCAAGTATTTc     <  1:376200/135‑1 (MQ=255)
                                                                                                                             tgcctCTCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGACATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATATACGGCTGGCGTGGCAGTTCGCCCAAGTATTTCATg  <  3:86417/135‑1 (MQ=255)
                                                                                                                             tgcctCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGCCCAAGTATTTCATg  >  3:72367/1‑135 (MQ=255)
                                                                                                                             tgcctCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGCCCAAGTATTTCATg  <  4:129098/135‑1 (MQ=255)
                                                                                                                             tgcctCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGGGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGCCCAAGTATTTCATg  <  6:101478/135‑1 (MQ=255)
                                                                                                                               |                                                                                                                                    
TGTTTGGCGAGAACAACCCAGAAAACCTGCACTTACTGCCCGATCCGTTGCTGCGGCCGCTGGCGCACCTGATGATCGACGAGTTCCAGGACGTGTCGCCGCAGATCGTCAGCTGGTTGCGTGCCTGCCTCGCTGAAATACGGCGCCGCGGCCCGGCCATGCACACCGGGCGCCATGCCCAGCATTCGTCGTTGCTGTGCGTGGGTGACGACTGGCAATCGATCTACGGCTGGCGTGGCAGTTCGCCCAAGTATTTCATG  >  KT2440/2921331‑2921590

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