Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,468,207 G→A 5.3% S138S (AGC→AGT)  proA_2 ← Gamma‑glutamyl phosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,468,2070GA 5.3% 260.1 / 3.1 76S138S (AGC→AGT) proA_2Gamma‑glutamyl phosphate reductase
Reads supporting (aligned to +/- strand):  ref base G (37/35);  new base A (2/2);  total (39/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 = 9.49e-01

ACCACCTGCACCACCGCCGCCGGCAGGCCGGCTTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGACGGTAGCTCATGTCGCGGATGGCACCGACTGGGTC  >  KT2440/5468073‑5468338
                                                                                                                                      |                                                                                                                                   
accaccTGCACCACCGCCGCTGGCAGGCCGGCTTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGc                                                                                                                                    >  7:330939/1‑136 (MQ=255)
  caccTGCACCACCGCCGCCGGCAGGCCGGCTTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCt                                                                                                                                   <  2:410647/135‑1 (MQ=255)
  caccTGCACCACCGCCGCCGGCAGGCCGGCTTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCt                                                                                                                                   <  2:293963/135‑1 (MQ=255)
  caccGGCACCACCGCCGCCGGCAGGCCGGCTTCGGCCAGGCGACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCt                                                                                                                                   <  1:450946/135‑1 (MQ=255)
       gCACCACCGCCGCCGGCAGGCCGGCTTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACGTGAGGCACAGGCTggcgg                                                                                                                              >  1:365585/1‑135 (MQ=255)
        caccacCGCCGCCGGCAGGCCGGCTTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTggcggc                                                                                                                             >  3:117098/1‑135 (MQ=255)
               gccgccgGCAGGCCGGCTTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATg                                                                                                                      <  2:79599/135‑1 (MQ=255)
                ccgccgGCAGGCCGGCTTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCCGGCGGCATCAATGGt                                                                                                                    >  8:213081/1‑136 (MQ=255)
                   ccgGCAGGCCGGCTTCGGCCAGGCAACGCGGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGACCTCCTTGAGGCACAGGCTGGCGGCTTCAATGGTCa                                                                                                                  >  2:29118/1‑135 (MQ=255)
                   ccgGCAGGCCGGCGTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCAc                                                                                                                 >  7:369566/1‑136 (MQ=255)
                          gCCGGCTTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCgg                                                                                                           <  4:293681/135‑1 (MQ=255)
                          gCCGGCTTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCgg                                                                                                           <  4:96227/135‑1 (MQ=255)
                          gCCGGCTTCGGCCAGGCCACGCTGAATGCAGGGGGCAATGGCGCGGTTGGAGTGGATGGCTTCGTAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCgg                                                                                                           <  1:25171/135‑1 (MQ=255)
                            cGGCTTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGAc                                                                                                         <  6:270302/135‑1 (MQ=255)
                            cGGCTTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGGCACGTTCGGAc                                                                                                         <  2:108698/135‑1 (MQ=255)
                                       aggccaCGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTc                                                                                                  >  6:193755/1‑131 (MQ=255)
                                       aggccaCGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTc                                                                                                  >  2:66027/1‑131 (MQ=255)
                                        ggccaCGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAgat                                                                                            <  8:171477/136‑1 (MQ=255)
                                        ggccaCGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCCGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAgat                                                                                            <  7:213081/136‑1 (MQ=255)
                                         gccaCGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAgat                                                                                            <  6:101804/135‑1 (MQ=255)
                                            aCGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCa                                                                                                                            <  4:230074/100‑1 (MQ=255)
                                            aCGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATc                                                                                        >  7:247234/1‑136 (MQ=255)
                                            aCGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACAGACACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCa                                                                                                                            >  3:230074/1‑100 (MQ=255)
                                               cTGAATGCAGGTGGCAATGGCGCGGGTGGAGGGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATccc                                                                                      <  1:231438/135‑1 (MQ=255)
                                                tGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCGCGACTTGAGGCACATGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCg                                                                                     <  1:29118/135‑1 (MQ=255)
                                                tGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCg                                                                                     <  1:52338/135‑1 (MQ=255)
                                                      cAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCgg                                                                                                           >  2:93088/1‑107 (MQ=255)
                                                      cAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCgg                                                                                                           <  1:93088/107‑1 (MQ=255)
                                                      cAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCAccc                                                                              >  7:59550/1‑136 (MQ=255)
                                                      cAGGGGGCAATGGCGCGGTTGGAGTGGATGGCTTCGTAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCAccc                                                                              <  8:228454/136‑1 (MQ=255)
                                                       aGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCAcccc                                                                             >  7:205713/1‑136 (MQ=255)
                                                         gTGGCAATGGCGCGGTTGGAGTGGCTGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCa                                                                           <  7:319381/136‑1 (MQ=255)
                                                          tGGCTATGGCGCGGTTGGAGTGGATGGCTGCGGAACCGCCACGCAGGATGGGTGCATTGCCAGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCa                                                                           <  4:245383/135‑1 (MQ=255)
                                                           ggCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCCCGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGc                                                                         >  8:442725/1‑136 (MQ=255)
                                                           ggCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGCCGCGATTCATTGATGATCCCGATCACCCCCAg                                                                          >  1:227067/1‑135 (MQ=255)
                                                               aTGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCggg                                                                      <  4:270662/135‑1 (MQ=255)
                                                               aTGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCggg                                                                      >  5:299711/1‑135 (MQ=255)
                                                               aTGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCggg                                                                      >  2:340011/1‑135 (MQ=255)
                                                               aTGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCCACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCggg                                                                      >  2:144134/1‑135 (MQ=255)
                                                               aTGGCGCGGGGGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGTCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCggg                                                                      <  1:330172/135‑1 (MQ=255)
                                                                tGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCgggg                                                                     >  7:36621/1‑135 (MQ=255)
                                                                  gcgcGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGt                                                                    <  1:441267/134‑1 (MQ=255)
                                                                   cgcgGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTAc                                                                  <  4:421917/135‑1 (MQ=255)
                                                                     ctggtGGAGTGGATGGCATCGCAACCGCCACGCAGGATGGGTGCATTCCCCGACGTGAGGCACAGGCTGGAGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGc                                                                <  3:22898/131‑1 (MQ=255)
                                                                     cgGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGc                                                                <  2:316560/135‑1 (MQ=255)
                                                                      ggTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCa                                                               >  3:300330/1‑135 (MQ=255)
                                                                       gTTGGAGTGGATGGCTGCGGAACCGCCACGCAGGATGGGTGCATTGCCCGACTTGAGGCCCAGGCGGGCGGCTTCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCGCAGCGGGGTACGCAt                                                              <  6:294539/135‑1 (MQ=255)
                                                                          ggAGTGGAGGGCTTCGGAACCGCCACGCAGGATGGGTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCtt                                                           <  4:381306/135‑1 (MQ=255)
                                                                                     cTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGaa                                                >  4:220002/1‑135 (MQ=255)
                                                                                     cTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGaa                                                >  6:369916/1‑135 (MQ=255)
                                                                                     cTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGaa                                                >  3:96731/1‑135 (MQ=255)
                                                                                      ttCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAAt                                               <  6:299711/135‑1 (MQ=255)
                                                                                        cGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGc                                             >  4:197086/1‑135 (MQ=255)
                                                                                        cGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCACCGGGGAACGCATCTTCCCAACCAGAATGc                                             >  2:227963/1‑135 (MQ=255)
                                                                                           aaCCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCg                                          >  3:144086/1‑135 (MQ=255)
                                                                                                 cACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGAc                                    <  4:144086/135‑1 (MQ=255)
                                                                                                  aCGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGACg                                   >  6:357254/1‑135 (MQ=255)
                                                                                                   cGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGACgg                                  >  7:52754/1‑135 (MQ=255)
                                                                                                   cGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGACg                                   >  4:338208/1‑134 (MQ=255)
                                                                                                   cGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGACTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCAcccc                                                                             >  6:281944/1‑92 (MQ=255)
                                                                                                   cGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGACTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCAcccc                                                                             <  5:276639/92‑1 (MQ=255)
                                                                                                   cGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGACTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCAcccc                                                                             <  5:281944/92‑1 (MQ=255)
                                                                                                   cGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGACTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCAcccc                                                                             >  6:276639/1‑92 (MQ=255)
                                                                                                                cATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGc                                                                <  5:325124/92‑1 (MQ=255)
                                                                                                                cATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGc                                                                >  6:325124/1‑92 (MQ=255)
                                                                                                                cATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGACGGTAGCTCATGTcgc                     >  4:227155/1‑135 (MQ=255)
                                                                                                                 aTTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGACGGTAGCTCATGTcgc                     >  1:421441/1‑134 (MQ=255)
                                                                                                                     cccGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGACGGTAGCTCATGTCGCGGATg                <  7:135802/135‑1 (MQ=255)
                                                                                                                     cccGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGACGGTAGCTCATGTCGCGGATg                >  2:14290/1‑135 (MQ=255)
                                                                                                                        gACTTGAGGCACAGGCTGGCGGCATCATTGCTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCACAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGACGGTAGCTCATGTCGCGGATGGCAc            <  8:94303/136‑1 (MQ=255)
                                                                                                                        gACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCggg                                                                      >  6:5625/1‑78 (MQ=255)
                                                                                                                        gACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCggg                                                                      <  5:5625/78‑1 (MQ=255)
                                                                                                                        gACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGACGGTAGCTCATGTCGCGGATGGCAc            >  7:281718/1‑136 (MQ=255)
                                                                                                                           tagAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGACGGTAGCTCATGTCGCGGATGGCACCg          <  1:109304/133‑1 (MQ=255)
                                                                                                                                  acaGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGACGGTAGCTCATGTCGCGGATGGCACCGACTGGGt   >  3:33549/1‑135 (MQ=255)
                                                                                                                                   caGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATTGACGGTAGCTCATGTCTCGGATTGCACCGACTGGGTc  >  4:68471/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                   
ACCACCTGCACCACCGCCGCCGGCAGGCCGGCTTCGGCCAGGCCACGCTGAATGCAGGTGGCAATGGCGCGGTTGGAGTGGATGGCTTCGGAACCGCCACGCAGGATGGTTGCATTGCCCGACTTGAGGCACAGGCTGGCGGCATCAATGGTCACGTTCGGACGCGATTCATAGATGATCCCGATCACCCCCAGCGGGGTACGCATCTTGCCTACCTGAATGCCCGATGGACGGTAGCTCATGTCGCGGATGGCACCGACTGGGTC  >  KT2440/5468073‑5468338

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