Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,729,011 T→G 21.1% Q190P (CAG→CCG)  csdA ← Cysteine desulfurase CsdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,729,0110TG21.1% 120.8 / 3.0 78Q190P (CAG→CCG) csdACysteine desulfurase CsdA
Reads supporting (aligned to +/- strand):  ref base T (11/45);  new base G (13/3);  total (30/48)
Fisher's exact test for biased strand distribution p-value = 1.10e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAGAGCACGCCGACGCCGTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCCTGGTCCTGGTCGATACGAC  >  ACB122/1728878‑1729140
                                                                                                                                     |                                                                                                                                 
aTAGAGCACGCCGACGCCGTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTgg                                                                                                                                 <  7:366798/136‑1 (MQ=255)
aTAGAGCACGCCGACGCCGTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGGCAGCGCGCCCTgg                                                                                                                                 <  7:11979/136‑1 (MQ=255)
      cACGCCGACGCCGTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCgg                                                                                                                           <  3:185864/136‑1 (MQ=255)
      cACGCCGACGCCGTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCGGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCgg                                                                                                                           <  3:249277/136‑1 (MQ=255)
       acgccgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGGCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCgg                                                                                                                           <  1:226910/135‑1 (MQ=255)
       acgccgacgccgTCCGGACCGTACAGCGTGTGGCTGGAAAACACATAGACGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCCCCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCgg                                                                                                                           <  5:191459/135‑1 (MQ=255)
       acgccgacgccgGCCGGACCGTACAGCTGGTGGCGGGAAAACACAAACAAGGCGCCGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCgg                                                                                                                           <  5:254214/135‑1 (MQ=255)
         gccgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCCCCATCGCCCCCCGTCAGCGCCCCCGGGGCGGGGGc                                                                                                                         >  6:157745/1‑135 (MQ=255)
           cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCCCCATCGACCACCGTCAGCGCGCCCTGGGCGCggggagg                                                                                                                      >  7:238638/1‑132 (MQ=255)
           cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCCTCGACCCCCGTCAGCCCCCCCGCGGGGCGGGCCTg                                                                                                                      >  3:335269/1‑136 (MQ=255)
           cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCCTCGACCACCCTCAGCCCCCCCGGGGCGCGGGCAtt                                                                                                                      >  7:353753/1‑135 (MQ=255)
           cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCCGCGCGCCCGGGGCGCCGGcct                                                                                                                       >  1:14202/1‑133 (MQ=255)
           cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGGCCGGGGCGCGGGcct                                                                                                                       >  5:285703/1‑133 (MQ=255)
           cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGcct                                                                                                                       >  5:155833/1‑133 (MQ=255)
           cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGCAt                                                                                                                       >  1:35189/1‑135 (MQ=255)
           cgacgccgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGATGCATGTCCTGGCGACCATGAACCACGCCCGGGGCACCATCGACCACCCACCCGCCGCCCCGGGCGCGGGcca                                                                                                                       >  1:349504/1‑133 (MQ=255)
                 cgTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCAGGGCGGGGGCATGGGccc                                                                                                                 >  2:132462/1‑134 (MQ=255)
                      ggACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCTCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACTCCCTGGGCACCATCGACCACCGACAGCCCCCCCAGGCCGCGGGAATGGGCCcccagc                                                                                                            >  1:79418/1‑135 (MQ=255)
                      ggACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAgcg                                                                                                           >  8:51677/1‑136 (MQ=255)
                      ggACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCCGGGCGCGGGCCTGGGCCcgcagc                                                                                                            >  1:325341/1‑135 (MQ=255)
                              cAGCTTGTGGCTGGAAAACACCTCGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAg                                                                                                    <  8:340486/135‑1 (MQ=255)
                              cAGCTTGTGGCGGGAAAACACATAGAAGTCGCAGCCCAGTAGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCg                                                                                                  <  3:433/136‑1 (MQ=255)
                               aGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAg                                                                                                    <  2:237794/134‑1 (MQ=255)
                               aGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCg                                                                                                  <  7:352614/136‑1 (MQ=255)
                               aGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGGTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGAGCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCg                                                                                                  <  7:234152/136‑1 (MQ=255)
                                     tgGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGGGGGCATGGGCCAGCAGCGCTGGCAGCGGtt                                                                                               <  2:229895/133‑1 (MQ=255)
                                       gCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCCGCACCCCTGGCCGCGGGTGCCAg                                                                                          >  8:74270/1‑136 (MQ=255)
                                           gAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGt                                                                                        <  7:25833/134‑1 (MQ=255)
                                           gAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTAc                                                                                      <  3:330123/136‑1 (MQ=255)
                                           gAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTAc                                                                                      <  4:42697/136‑1 (MQ=255)
                                            aaaaCACATTGAAGGCGCCGCCCTGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTAc                                                                                      <  1:176897/135‑1 (MQ=255)
                                            aaaaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTAcc                                                                                     >  8:229608/1‑136 (MQ=255)
                                            aaaaCACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTAc                                                                                      <  6:381300/135‑1 (MQ=255)
                                                             cAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGGCTGGGCCCGCCGCGCGGGCCGCGGGTGCCCGGTTCCCCGCACGTTGGGCAgct                                                                    >  7:200582/1‑136 (MQ=255)
                                                               gCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGGGCGGGCATGGGCCCGCAGCGCGGGCCGCGGTTGCCAGGTAACCCGCACCTTGGACCgctgg                                                                  >  4:142165/1‑136 (MQ=255)
                                                                cccAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCGGGCCGCGGTTGCCAGGTACCCAGCACGTTGGACAgctggc                                                                 >  7:205296/1‑136 (MQ=255)
                                                                    gTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGAt                                                             <  3:294786/136‑1 (MQ=255)
                                                                     ttGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCGCGGGGCGCGGGCTTGGGCCAGCAGCGCTGGGAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGAt                                                             <  1:132462/135‑1 (MQ=255)
                                                                     ttGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGAt                                                             <  2:114121/135‑1 (MQ=255)
                                                                     gtgctgcATGGCCAGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGAt                                                             <  1:183990/134‑1 (MQ=255)
                                                                             tGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCCCCCGGGGCGCGGGCATGGGCCAGAAGCGCGGGCAGCGGTTGGCAGGTACCCAGCACGTTGGACAGCTGGCGGATGGCGAGCAc                                                    >  3:154602/1‑136 (MQ=255)
                                                                                    gCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGt                                              <  8:353753/135‑1 (MQ=255)
                                                                                      gACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTgc                                            <  2:213914/135‑1 (MQ=255)
                                                                                        ccATGCACCCCGCCCTGGGCACCAGCGACCACCGGCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGcgcg                                         <  8:366167/136‑1 (MQ=255)
                                                                                         cATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCAGGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGcgcg                                         <  5:157745/135‑1 (MQ=255)
                                                                                          aTGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCggg                                       <  3:292093/136‑1 (MQ=255)
                                                                                          aTGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCggg                                       <  3:78947/136‑1 (MQ=255)
                                                                                              accacGCCCTGGGCACCCTCGCCCCCCGTCAGCCCGCCCCGGGGCCGGGCCCGGGCCCCCCGCCCAGGCCCCGGTTGCCAGGTACCCCGCCCCGTGGACACCAGGCTGATGGCGAGCCATCGGGTGCGCGGGccc                                    >  4:66593/1‑134 (MQ=255)
                                                                                              accacGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGcc                                     <  5:349605/134‑1 (MQ=255)
                                                                                               ccacGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGAt                                  <  4:50708/136‑1 (MQ=255)
                                                                                                   gCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGt                              <  4:332743/136‑1 (MQ=255)
                                                                                                    cccTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGtt                             <  8:43282/136‑1 (MQ=255)
                                                                                                    cccTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGtt                             <  3:63426/136‑1 (MQ=255)
                                                                                                     ccTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCTATCACTCTTTTGCGCGGGCCGATCAGtt                             <  5:273669/135‑1 (MQ=255)
                                                                                                     ccTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGtt                             <  1:143465/135‑1 (MQ=255)
                                                                                                     ccTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGGGCGGGCATGGGCCAGCAGCGCTGGCAGCGGGTGCCACGTACCCAGCACGTTGGACAGCTGGCTGATTGCGAGCACTCGGGTGGGCGGGCCGAGCAGTTg                            >  4:288060/1‑136 (MQ=255)
                                                                                                     ccTGGGCACCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGCCTGGGCCCCCCGCGCTGGCCGCGGGTGGCAGGTACCCCGCCCGTTGGCCAGCGGGCTGGTGGCGAGCCCTCGGGGGCGCGGGCCGCTGAGTTg                            >  4:7232/1‑136 (MQ=255)
                                                                                                     ccTCGGCACCCTCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGtt                             <  2:349504/135‑1 (MQ=255)
                                                                                                         ggCACCATCGCCCACCGTCAGCGCGCCCGGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGt                        <  3:252824/136‑1 (MQ=255)
                                                                                                         ggCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAg                         <  7:376814/135‑1 (MQ=255)
                                                                                                          gCACCATCGACGACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGTCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGCTGCAGt                        <  2:79418/135‑1 (MQ=255)
                                                                                                          gCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGt                        <  2:35189/135‑1 (MQ=255)
                                                                                                          gCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGt                        <  2:367310/135‑1 (MQ=255)
                                                                                                           cACCATCGACCACCGTCAGCGCGCCCGGGGCGCGGGCATGGGCCAGCCGCGCGGGCAGCGGGTGCCAGGTACCCAGCACCTTGGACAGCTGGGGGATGGCGCGCACTCGGGGGCGCGGGCCGAGCAGGTGCAGTGc                      >  8:188107/1‑136 (MQ=255)
                                                                                                            aCCATCGACCACCGTCAGCCCGCCCGGGGCGCGGGGATGGGCCCGCCGCGCTGGCAGCCGTTGCCAGGTACCCCGCACGTTGGACAGCTGGCTGATGGCGCGCCCTCGGGTGGGCGGGCCGATCAGTTTCAGTGcc                     >  4:263594/1‑136 (MQ=255)
                                                                                                                 cgaactccGTCAGCGGGCCCTGGGGGCGGGCATGGGGCTGCATCGCTGGCATCGGTTGCCAGGTACCCAGCACGTTGGAGAGCTGGCTTATGGTGAGCACTCGTGTGCGCGGGCCGATCAGTTGCAGTGcctggtc                <  8:291071/130‑1 (MQ=255)
                                                                                                                 cGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCCGCACGTTGGACAGCTGGCTGCTGGCGAGCCCTCGGGGGCGCGGGCCGCTCAGTTGCAGTGcctggt                 >  5:360847/1‑135 (MQ=255)
                                                                                                                 cGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGCTCAGTTGCAGTGccgggt                 >  6:299470/1‑135 (MQ=255)
                                                                                                                 cGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGcctggtc                >  7:392237/1‑136 (MQ=255)
                                                                                                                   accaccGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGcctggtcct              <  3:288060/136‑1 (MQ=255)
                                                                                                                   accaccGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGcctggtcct              <  8:392237/136‑1 (MQ=255)
                                                                                                                      accGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGCGCACTCGGGTGCGCGGGCCGATCAGTTTCAGTGcctggtcctggt           >  3:86279/1‑136 (MQ=255)
                                                                                                                       ccGTCAGCGCGCCCGGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCctggtcctggtc          <  3:290393/136‑1 (MQ=255)
                                                                                                                       ccGGCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCctggtcctggtc          <  7:309955/136‑1 (MQ=255)
                                                                                                                         gTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCCTGGTCCTGGTCg         <  6:285703/135‑1 (MQ=255)
                                                                                                                             gcgcgcCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGt                                              <  5:88716/94‑1 (MQ=255)
                                                                                                                             gcgcgcCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGt                                              >  6:88716/1‑94 (MQ=255)
                                                                                                                                 gcCCAGGGCGCGGGCATGGGCCAGCAGCCCTGGCAGCGGTTGCCAGGTACCCAGCCCCGTGGCCAGCTGGCGGATGGCGCGCACTCGGGTGCGCGGGGCGCTCAGTTGGAGTGCCTGGTCCTGGTCGATACGAc  >  6:302487/1‑134 (MQ=255)
                                                                                                                                     |                                                                                                                                 
ATAGAGCACGCCGACGCCGTCCGGACCGTACAGCTTGTGGCTGGAAAACACATAGAAGTCGCAGCCCAGTTGCTGCATGTCCTGGCGACCATGCACCACGCCCTGGGCACCATCGACCACCGTCAGCGCGCCCTGGGCGCGGGCATGGGCCAGCAGCGCTGGCAGCGGTTGCCAGGTACCCAGCACGTTGGACAGCTGGCTGATGGCGAGCACTCGGGTGCGCGGGCCGATCAGTTGCAGTGCCTGGTCCTGGTCGATACGAC  >  ACB122/1728878‑1729140

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: