Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,778,681 C→T 5.5% T343T (ACC→ACT selO → Protein adenylyltransferase SelO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,778,6810CT 5.5% 242.5 / 3.1 73T343T (ACC→ACTselOProtein adenylyltransferase SelO
Reads supporting (aligned to +/- strand):  ref base C (42/26);  new base T (2/2);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

ACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAGGCGCTGAAGGTGGCGCGCGAT  >  ACB122/5778548‑5778813
                                                                                                                                     |                                                                                                                                    
aCCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCAccg                                                                                                                                     >  1:96163/1‑135 (MQ=255)
      gcgcgcTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGTCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGa                                                                                                                                 <  7:243285/133‑1 (MQ=255)
       cgcgcTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGTAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAga                                                                                                                              <  6:184810/135‑1 (MQ=255)
       cgcgcTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAga                                                                                                                              >  3:50284/1‑135 (MQ=255)
       cgagatGGACCAGGCACTGACCACCGTGATCGTTCTTTAGCGTCTGATGGTTGGGTGGGGCCTGTTCGTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGGGGGCGCGTAGGGCTGATCACAGCCGTAga                                                                                                                              <  8:444466/130‑1 (MQ=255)
        gcgcTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCCTTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAga                                                                                                                              >  6:168882/1‑134 (MQ=255)
         cgcTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGCTGCGCCGGCGCCTAGGGCTGACCCCCGCCGAAgatg                                                                                                                            >  2:145687/1‑135 (MQ=255)
           cTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAgatgat                                                                                                                          <  5:27542/135‑1 (MQ=255)
             ggCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGAtgatga                                                                                                                        <  3:429865/135‑1 (MQ=255)
              gCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGAc                                                                                                                       >  6:85672/1‑135 (MQ=255)
               cccAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGCCCTGATGCGCCGGCGCCTAGGGCTGCCCCCCGCCGAAGAGGATGACa                                                                                                                      >  1:420100/1‑135 (MQ=255)
                      aCTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCt                                                                                                               >  7:69627/1‑135 (MQ=255)
                             accGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAg                                                                                                        <  3:184682/135‑1 (MQ=255)
                             accGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAg                                                                                                        >  3:229237/1‑135 (MQ=255)
                             accGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAg                                                                                                        >  3:25889/1‑135 (MQ=255)
                             accGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAg                                                                                                        >  4:46179/1‑135 (MQ=255)
                             accGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAg                                                                                                        >  1:131012/1‑135 (MQ=255)
                              ccGTGATCGAAGTTGAGCCACTGAAGGAAGCGGTGGGCCTGTTCCTGCCGCTGTACCAGGCCCTTTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAgc                                                                                                       <  6:155907/135‑1 (MQ=255)
                                  gATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAga                                                                                                                              <  2:352509/108‑1 (MQ=255)
                                  gATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAga                                                                                                                              >  1:352509/1‑108 (MQ=255)
                                       aaGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCGGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCtgctgc                                                                                              <  6:348436/135‑1 (MQ=255)
                                          tggAGCCACTGACGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGt                                                                                           <  7:321004/133‑1 (MQ=255)
                                            gAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGtgc                                                                                         <  6:202477/135‑1 (MQ=255)
                                                cACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGtgcatg                                                                                      >  6:211948/1‑134 (MQ=255)
                                                cACTGAAGGAAGCGTTGGGCCTGTTCCTGCAGCTGTACCACGACCATTACCTGGACCTGATGCTCCGGCGCCTAGGGCTGACCACCGCCCAAGATGATGACATGGTGCTGGTCGATCGCCTGCTGCACtgcatg                                                                                      >  8:62790/1‑134 (MQ=255)
                                                   tGAAGGAAGCGTTGGGCCGGTTCCTGCCGCTGTACCAGGCCCATTACCTGGCCCTGACGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAAc                                                                                  <  1:240971/135‑1 (MQ=255)
                                                      aggaagCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAAcg                                                                                 <  4:460508/133‑1 (MQ=255)
                                                         aagCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTg                                                                            <  4:25889/135‑1 (MQ=255)
                                                          agCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTgg                                                                           <  2:243979/135‑1 (MQ=255)
                                                             gTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGt                                                                        <  8:147451/135‑1 (MQ=255)
                                                             gTTGGGCCTGTTCCTGCCGCTGTACAAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGt                                                                        >  6:165152/1‑135 (MQ=255)
                                                               tGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAAc                                                                                  <  8:280334/123‑1 (MQ=255)
                                                               tGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAAc                                                                                  >  7:280334/1‑123 (MQ=255)
                                                               tGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTgg                                                                      >  8:44189/1‑135 (MQ=255)
                                                                  gCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGAtgatga                                                                                                                        >  2:53435/1‑82 (MQ=255)
                                                                  gCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGAtgatga                                                                                                                        <  1:53435/82‑1 (MQ=255)
                                                                         ccTGCCGCTGTACCAGGCCCATTACCTGGCCCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCt                                                            >  4:140512/1‑135 (MQ=255)
                                                                         ccTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACTGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCt                                                            >  3:115083/1‑135 (MQ=255)
                                                                         ccTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACTGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCt                                                            >  3:108729/1‑135 (MQ=255)
                                                                         ccTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCt                                                            >  6:383962/1‑135 (MQ=255)
                                                                          cTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCCCCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGCCTACAGCCTg                                                           >  1:306967/1‑135 (MQ=255)
                                                                          cTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACTGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCt                                                            <  4:115083/134‑1 (MQ=255)
                                                                          cTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACTGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCt                                                            <  4:108729/134‑1 (MQ=255)
                                                                          cTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTg                                                           >  2:36841/1‑135 (MQ=255)
                                                                              cGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGttc                                                        <  1:36841/134‑1 (MQ=255)
                                                                               gCTGTACCAGGCCCATTACCTGGACCTGATGCGCAGGCGCCTAGGGCTGCCCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACCGCCTGgtctt                                                      >  3:201021/1‑135 (MQ=255)
                                                                                    aCCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCGTGTTCTTCCg                                                   <  3:46179/133‑1 (MQ=255)
                                                                                      cAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATCGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCGTCCGCaaa                                               >  2:41559/1‑135 (MQ=255)
                                                                                      cAGGCCCATTACCTGGACATGATGCGCCGGCGCCTAGAGCTGACCACCGCCGAAGATGAAGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCAGCaaa                                               >  2:41553/1‑135 (MQ=255)
                                                                                           ccATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATTATGACATGGTGCTGGTCGAGCGCctg                                                                                                  >  1:226184/1‑79 (MQ=255)
                                                                                           ccATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCctg                                                                                                  <  2:226184/79‑1 (MQ=255)
                                                                                           ccATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTggg                                          >  1:65485/1‑135 (MQ=255)
                                                                                                aCCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGa                                       >  6:433573/1‑133 (MQ=255)
                                                                                                aCCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGa                                       >  5:206594/1‑133 (MQ=255)
                                                                                                aCCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGa                                       >  2:58869/1‑133 (MQ=255)
                                                                                                aCCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGAGCATGCAACGCGGTGGCGTGGACGACAGCCTGATCTTCCGCAAACTGGGGGa                                       >  6:332743/1‑133 (MQ=255)
                                                                                                aCCTGGACCTGATGCGCCGGCGCCTAGGGCTGACAACCGCCGAAGATGATGACATGGTGCTGGTAGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGAGTGGACTACAGCCAGTTCTTACGCACACTGGGTGagc                                     >  1:304038/1‑135 (MQ=255)
                                                                                                       ccTGATGCGCCGGCGCCTAGGGCTGCCCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCg                                >  7:296162/1‑133 (MQ=255)
                                                                                                        cTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTa                             >  7:328187/1‑135 (MQ=255)
                                                                                                        cTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTa                             >  2:5291/1‑135 (MQ=255)
                                                                                                        cTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATCACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTa                             >  2:13379/1‑135 (MQ=255)
                                                                                                         tGATGCGCCGGCGCCTCGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAg                            <  4:427904/135‑1 (MQ=255)
                                                                                                         tGATGCGCCGGCGCCGAGGGCGGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAg                            <  1:207937/135‑1 (MQ=255)
                                                                                                             gcgcCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGCCATGGTGCTGGTCGAGCGCCTGCTGCAGAGCATGCAACGCGGTGGCGAGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGa                        >  4:321420/1‑135 (MQ=255)
                                                                                                             gcgcCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGa                        >  4:115368/1‑135 (MQ=255)
                                                                                                              cgcCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACCGCCTGGTCTTCCGCAAACGGGGTGAGCAGCCGGTAGCCGAg                       >  1:113960/1‑135 (MQ=255)
                                                                                                              cgcCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAg                       >  2:88422/1‑135 (MQ=255)
                                                                                                              cgcCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAg                       >  4:15017/1‑135 (MQ=255)
                                                                                                              cgcCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTATTACGCAAACTGGGTGAGCAGCCGGTAGCCgcg                       >  7:167894/1‑133 (MQ=255)
                                                                                                               gcCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAg                       <  6:206594/134‑1 (MQ=255)
                                                                                                                    cgcCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAGGCGCTg                 <  7:44189/135‑1 (MQ=255)
                                                                                                                     gcCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAGGCGCt                  >  8:12203/1‑133 (MQ=255)
                                                                                                                        tcgggCTGACCACCGCCGAAGAGGAGGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGGGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAGGCGCTGAAgg             <  6:44878/133‑1 (MQ=255)
                                                                                                                        tAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAGGCGCTGAAgg             <  5:211948/135‑1 (MQ=255)
                                                                                                                                 ccaccGCCTAAGATGAGGACATGGTGCTGGGCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAGGCGCTGAAGGTGgcgcg      <  4:96208/133‑1 (MQ=255)
                                                                                                                                 ccaccGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAGGCGCTGAAGGTGgcgcg      <  3:15017/133‑1 (MQ=255)
                                                                                                                                  caccGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAGGCGCTGAAGGTGGCGCGCGa   <  5:168882/135‑1 (MQ=255)
                                                                                                                                   accGCCGAAGATGATGACATGGTGCTGGTCGAGCTCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACGACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAGGCGCTGAAGGTGGCGCGCGa   >  3:448340/1‑134 (MQ=255)
                                                                                                                                   accGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAGGCGCTGAAGGTGGCGCGCGa   >  7:209790/1‑134 (MQ=255)
                                                                                                                                   accGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAGGCGCTGAAGGTGGCGCGCGa   >  6:302695/1‑134 (MQ=255)
                                                                                                                                   accGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAGGCGCTGAAGGTGGCGCGCGAt  >  4:235773/1‑135 (MQ=255)
                                                                                                                                     |                                                                                                                                    
ACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGCAGCCGGTAGCCGAGGCGCTGAAGGTGGCGCGCGAT  >  ACB122/5778548‑5778813

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: