Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,201,776 G→A 6.4% R301W (CGG→TGG)  ghrA ← Glyoxylate/hydroxypyruvate reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,201,7760GA 6.4% 226.8 / 2.4 72R301W (CGG→TGG) ghrAGlyoxylate/hydroxypyruvate reductase A
Reads supporting (aligned to +/- strand):  ref base G (38/28);  new base A (2/2);  total (42/30)
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.85e-01

GGCCTCTTCGCGGGCAAGCCCGCTCCTACGGGGGGTATCACCACCCTGTGCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCGGTCCAGAACGGGTGTTGCTTGGGCAA  >  ACB122/5201650‑5201911
                                                                                                                              |                                                                                                                                       
ggCCTCTTCGCGGGCAAGCCCGCTCCTACGGGGGGTATCACCACCCTGTGCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTg                                                                                                                                <  3:355625/136‑1 (MQ=255)
  ccgctgcgcgGCCAAGCCCGCTCATACGGGGGGTATCACCACCCTGTGCTAAAGCGGTTTTACCGCTACGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGcc                                                                                                                              <  3:237845/130‑1 (MQ=255)
  ccTCTTCGCGGGCAAGCCCGCTCCTAGGGGGGGTAGCACCACCCTGTGCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGcc                                                                                                                              <  3:242623/136‑1 (MQ=255)
   ctctTCGCGGGCAAGCCCGCTCCTACGGGGGGTATCACCACCCTGTGCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGcc                                                                                                                              <  2:216156/135‑1 (MQ=255)
     ggtCGCGGGCAAGCCCGCTCCAACGGGGGGTATCCCCCCCCGGTGGAAAACCGGTTTTAACCCCAAGAGGCCGGGAACGGCCCCAGCCGCCCGTAGCCCCGGGCGCAATCGACATCACCCCCCACCCCCAGcccc                                                                                                                            >  6:307407/3‑135 (MQ=255)
         gcgGGCCAGCCCGCTCCTAAGGGGGGGATCACCACCCTGGGATAAAGCCGGTTTACAGCGCAGACGCCGGTACCGGCCGCAGCGAACAGTAGGCCCCGGCGCAATCGACGTCACCCCGCAACCCCTGCCCCGCGTc                                                                                                                       >  8:75274/1‑136 (MQ=255)
                        ccTACGGGGGGTATCACCACCCTGTGCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACAtt                                                                                                        >  4:184573/1‑136 (MQ=255)
                          tACGGGGGGTATCACCACCCTGTGCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATtct                                                                                                      >  4:290274/1‑136 (MQ=255)
                          tACGGGGGGTATCACCACCCTGTGCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATtct                                                                                                      >  7:112909/1‑136 (MQ=255)
                               ggggTATCACCACCCTGTGCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCGGCCCCGCTTCATACTCCCGCACATTCTCGa                                                                                                   >  1:122754/1‑134 (MQ=255)
                                    aTCACCACCCTGTGCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCATCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGacaaa                                                                                               >  5:271236/1‑133 (MQ=255)
                                    aTCACCACCCTGTGCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGacaaac                                                                                              >  8:208416/1‑134 (MQ=255)
                                       accaccCTGTGCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAgc                                                                                          >  2:129881/1‑135 (MQ=255)
                                              tgtgCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCACCCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCAt                                                                                   >  5:211035/1‑135 (MQ=255)
                                              tgtgCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGacaaac                                                                                              <  5:391819/124‑1 (MQ=255)
                                              tgtgCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGacaaac                                                                                              >  6:391819/1‑124 (MQ=255)
                                              tgtgCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATc                                                                                  >  7:303499/1‑136 (MQ=255)
                                                tgCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGCCAAACAAGCGCACCATCGc                                                                                >  8:334891/1‑136 (MQ=255)
                                                 gCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCt                                                                               >  8:78331/1‑136 (MQ=255)
                                                  cTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCACCCCCTGCCCCGCTTCATAAGCCCGCCCATTCTCGCCAAACAAGCGCACCATCGCt                                                                               >  5:312814/1‑135 (MQ=255)
                                                   tAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCACCCCCTGCCCCGCTTCATAAGCCCGCACAGTCTCGACAAACAAGCGCACCATCGCTg                                                                              >  1:199132/1‑135 (MQ=255)
                                                   tAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCt                                                                                                                         <  3:357753/92‑1 (MQ=255)
                                                   tAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCt                                                                                                                         >  4:357753/1‑92 (MQ=255)
                                                   tAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGt                                                                             >  3:197950/1‑136 (MQ=255)
                                                   tAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGt                                                                             >  3:76329/1‑136 (MQ=255)
                                                   tAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGt                                                                             >  4:120872/1‑136 (MQ=255)
                                                   tAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGt                                                                             >  4:181056/1‑136 (MQ=255)
                                                   tAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCAACGCTg                                                                              >  5:294465/1‑135 (MQ=255)
                                                   tAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAACCGCACCATCGCTg                                                                              >  6:278763/1‑135 (MQ=255)
                                                         ggTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCgg                                                                          >  1:39337/1‑133 (MQ=255)
                                                         ggTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGa                                                                       <  7:208416/136‑1 (MQ=255)
                                                         ggTTTTACCGCGAAGAGGACGGTACAGGAAGCATCGATCAGTAGCCAAGTTCGAAATAGACTTCAACACCCAACCCATGACCCGCTTCATAAGACCGCACATTCTAGAAAAACAAGAGCACCATAGCTGtatga                                                                         >  8:188953/1‑130 (MQ=255)
                                                          gTTTAACCGCTAAGCGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGa                                                                       <  6:312814/135‑1 (MQ=255)
                                                           ttttACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACTTTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAgg                                                                     <  4:69694/136‑1 (MQ=255)
                                                           ttgtcccGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAgg                                                                     <  4:8584/131‑1 (MQ=255)
                                                             ttACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCCACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTg                                                                   >  3:317291/1‑136 (MQ=255)
                                                             ttACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCACAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGt                                                                    >  5:123375/1‑135 (MQ=255)
                                                             ttACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCACAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTg                                                                   >  7:282435/1‑136 (MQ=255)
                                                                ccGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAACTCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTgg                                                                  <  1:144153/134‑1 (MQ=255)
                                                                ccGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCATCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTggg                                                                 <  6:271236/135‑1 (MQ=255)
                                                                ccGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGgc                                                                <  8:305283/136‑1 (MQ=255)
                                                                    gAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCACAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCg                                                            <  8:282435/136‑1 (MQ=255)
                                                                     aaGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCACAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCg                                                            <  6:123375/135‑1 (MQ=255)
                                                                         ggCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCtg                                                       <  3:120872/136‑1 (MQ=255)
                                                                                       gCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACACGCGCCCCATCGCTGTCGGCGAGGTGGGCGCCGCGCTGTGCCCGGTcagca                                          >  2:213275/1‑135 (MQ=255)
                                                                                       gCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTTTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTcagcag                                         >  4:70857/1‑136 (MQ=255)
                                                                                       gCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTcagcag                                         >  7:393313/1‑136 (MQ=255)
                                                                                       gCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTcagca                                          >  6:32960/1‑135 (MQ=255)
                                                                                       gCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGAGCCCGGTcagca                                          >  6:242330/1‑135 (MQ=255)
                                                                                            gATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAg                                     >  5:152197/1‑135 (MQ=255)
                                                                                            gATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACAATCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGAGCCCGGTCAGCAGCaac                                     >  5:278064/1‑134 (MQ=255)
                                                                                               cAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACACGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGcccc                                 <  7:149055/136‑1 (MQ=255)
                                                                                                aGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGcccc                                 <  8:303499/135‑1 (MQ=255)
                                                                                                aGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCa                                >  8:66041/1‑136 (MQ=255)
                                                                                                 gTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGccccc                                >  2:126266/1‑134 (MQ=255)
                                                                                                 gTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAg                               >  4:176822/1‑136 (MQ=255)
                                                                                                  tAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAgg                              <  8:301131/136‑1 (MQ=255)
                                                                                                  tAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAgg                              <  8:295988/136‑1 (MQ=255)
                                                                                                    gCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCg                            <  3:158522/136‑1 (MQ=255)
                                                                                                     ccccGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCgg                           >  4:126720/1‑136 (MQ=255)
                                                                                                     ccccGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCg                            <  6:294465/135‑1 (MQ=255)
                                                                                                        cGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCGGt                          >  6:368374/1‑134 (MQ=255)
                                                                                                        cGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCGGTc                         >  7:141055/1‑135 (MQ=255)
                                                                                                         gggCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCGGTCCa                       <  4:197950/136‑1 (MQ=255)
                                                                                                         gggCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCGGTCCa                       <  7:334891/136‑1 (MQ=255)
                                                                                                         gggCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCGGTCCa                       <  7:78331/136‑1 (MQ=255)
                                                                                                             gAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGacaaac                                                                                              <  8:154426/61‑1 (MQ=255)
                                                                                                             gAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGacaaac                                                                                              >  7:154426/1‑61 (MQ=255)
                                                                                                                         aCCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCGGTCCAGAACGGGTGTTGc          <  5:32960/133‑1 (MQ=255)
                                                                                                                         aCCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCGGTCCAGAACGGGTGTTGc          <  2:67758/133‑1 (MQ=255)
                                                                                                                            ccGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCCTCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCGGTCCAGAACGGGTGTTGCTTggg     <  2:250509/135‑1 (MQ=255)
                                                                                                                            ccGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCGGTCCAGAACGGGTGTTGCTTggg     <  5:351892/135‑1 (MQ=255)
                                                                                                                              gCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCGGTCCAGAACGGGTGTTGCTTGGGCaa  <  4:76329/136‑1 (MQ=255)
                                                                                                                              |                                                                                                                                       
GGCCTCTTCGCGGGCAAGCCCGCTCCTACGGGGGGTATCACCACCCTGTGCTAAAGCGGTTTTACCGCGAAGAGGCCGGTACAGGCCGCAGCGATCAGTAGCCCCGGGCGAAATCGACTTCACCCCGCAACCCCTGCCCCGCTTCATAAGCCCGCACATTCTCGACAAACAAGCGCACCATCGCTGTCGGCGAGGTGGGCGCCGAGCTGTGCCCGGTCAGCAGCAAGCCCCAGGCGGTCCAGAACGGGTGTTGCTTGGGCAA  >  ACB122/5201650‑5201911

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