Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB143 445,359 A→G H76R (CAC→CGC)  argD_1 → Acetylornithine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143445,3590AG97.0% 238.4 / ‑2.0 67H76R (CAC→CGC) argD_1Acetylornithine aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (38/27);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

CCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCACAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTACCGGCAGTGACCAGGC  >  ACB143/445226‑445491
                                                                                                                                     |                                                                                                                                    
ccTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCcgca                                                                                                                                    >  5:35998/1‑136 (MQ=255)
ccTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGAGGTTGTGGGACAACGAAGGGAATGCCTACCTGGACTTCCCTCAGGGCGGTGCGGTCAACAGCTTGGGCcgca                                                                                                                                    >  4:200017/1‑136 (MQ=255)
    aGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTTTTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGcc                                                                                                                                 <  7:105051/135‑1 (MQ=255)
    aGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCGCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGcc                                                                                                                                 <  1:130282/135‑1 (MQ=255)
       gcCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAg                                                                                                                             >  1:51233/1‑136 (MQ=255)
        ccTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGc                                                                                                                            <  3:332797/136‑1 (MQ=255)
        ccTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGc                                                                                                                            <  2:70516/136‑1 (MQ=255)
        ccTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGc                                                                                                                            <  3:101165/136‑1 (MQ=255)
            aTGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGc                                                                                                                        >  8:48326/1‑136 (MQ=255)
              gCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGcc                                                                                                                                 >  5:368676/1‑125 (MQ=255)
              gCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGcc                                                                                                                                 <  6:368676/125‑1 (MQ=255)
               cccGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTgg                                                                                                                     <  1:244633/136‑1 (MQ=255)
                  ggCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTCCCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCCGCGTGCTGGTc                                                                                                                   >  4:284293/1‑135 (MQ=255)
                      ccGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGc                                                                                                              >  8:225931/1‑136 (MQ=255)
                             tGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTTCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGc                                                                                                       <  5:311422/136‑1 (MQ=255)
                             tGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACATGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTggg                                                                                                        >  6:211792/1‑135 (MQ=255)
                                    cAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAgg                                                                                                <  6:337765/136‑1 (MQ=255)
                                    cAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAgg                                                                                                <  4:64536/136‑1 (MQ=255)
                                    cAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTCGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAgg                                                                                                >  5:56319/1‑136 (MQ=255)
                                             gtgcgtgGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCg                                                                                        <  7:1072/135‑1 (MQ=255)
                                               gcgtgGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGtt                                                                                      >  2:30228/1‑135 (MQ=255)
                                                    gCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCaa                                                                                >  2:304777/1‑136 (MQ=255)
                                                         ggCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCcgg                                                                           <  2:20931/136‑1 (MQ=255)
                                                            tCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCACCCcgggct                                                                        >  5:57242/1‑136 (MQ=255)
                                                                 gTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCcgggctcgggc                                                                   >  6:244311/1‑136 (MQ=255)
                                                                 gTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCcgggctcgggc                                                                   >  6:118823/1‑136 (MQ=255)
                                                                 gTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCcgggctcgggc                                                                   >  1:304099/1‑136 (MQ=255)
                                                                  ttGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCcgggctcgggct                                                                  >  8:120640/1‑136 (MQ=255)
                                                                    gtgGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCtt                                                                 <  5:191115/135‑1 (MQ=255)
                                                                      gggACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCca                                                              >  3:330496/1‑136 (MQ=255)
                                                                       ggACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCca                                                              >  5:168016/1‑135 (MQ=255)
                                                                           aaCGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAAcc                                                         >  1:279563/1‑136 (MQ=255)
                                                                            aCGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCg                                                        <  4:203524/136‑1 (MQ=255)
                                                                               aaGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTgg                                                     >  4:284964/1‑136 (MQ=255)
                                                                               aaGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTgg                                                     >  8:116734/1‑136 (MQ=255)
                                                                                aGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGc                                                    >  4:134582/1‑136 (MQ=255)
                                                                                aGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGc                                                    >  8:280316/1‑136 (MQ=255)
                                                                                 gggCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCt                                                   <  3:284964/136‑1 (MQ=255)
                                                                                 gggCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCt                                                   <  7:197999/136‑1 (MQ=255)
                                                                                          tACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGc                                          >  3:171682/1‑136 (MQ=255)
                                                                                               ggACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGt                                      >  3:376127/1‑135 (MQ=255)
                                                                                               ggACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGAGGCTTGCTGGGGCTGGt                                      >  3:242770/1‑135 (MQ=255)
                                                                                                gACCTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCa                                    <  5:270581/136‑1 (MQ=255)
                                                                                                  cTTCACTCCGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCCGGTCAAAGCGTTGGTCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAAc                                  <  3:276463/136‑1 (MQ=255)
                                                                                                    tCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCg                                >  8:199115/1‑136 (MQ=255)
                                                                                                     cACTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCgg                               >  5:218735/1‑136 (MQ=255)
                                                                                                      aCTCAGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCgg                               <  1:231299/135‑1 (MQ=255)
                                                                                                        tCACGGCGGTGCGGTCAACAGCTAGGGCCGCAGCCCCAGCGTGATCGTCAAAGCGTTGGGCAACCAGGATCACGCGTTGATCAACCCGGGCTCGGGATTCCACAACCGTGGCTTGCTGGCGCTGGTCAACCGGcag                            >  4:363867/1‑134 (MQ=255)
                                                                                                         cAGGGCGGTGCGGTCAACAGCTTGGGCCACAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGtt                                                                                      <  8:240520/77‑1 (MQ=255)
                                                                                                         cAGGGCGGTGCGGTCAACAGCTTGGGCCACAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGtt                                                                                      >  7:240520/1‑77 (MQ=255)
                                                                                                          aGGGCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCtgtg                          >  6:335894/1‑136 (MQ=255)
                                                                                                           gggCGGTGCGGTCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTgtgt                         <  1:304777/136‑1 (MQ=255)
                                                                                                               ggtgcggtCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCa                       >  6:362865/1‑134 (MQ=255)
                                                                                                                   cggtCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTAc                 >  4:193791/1‑136 (MQ=255)
                                                                                                                    ggtCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTAcc                >  5:261397/1‑136 (MQ=255)
                                                                                                                    ggtCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTAcc                >  8:15554/1‑136 (MQ=255)
                                                                                                                    ggtCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTAcc                >  6:180311/1‑136 (MQ=255)
                                                                                                                     gtCAACCGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAACGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTACCg               <  5:111481/136‑1 (MQ=255)
                                                                                                                     gtCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTACCg               <  5:22406/136‑1 (MQ=255)
                                                                                                                     gtCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTACCg               <  2:28433/136‑1 (MQ=255)
                                                                                                                     gtCAACAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTACCg               <  1:75965/136‑1 (MQ=255)
                                                                                                                         aCAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCACCCGGCTGTGTCAGAGTACCGGCAg           >  1:31576/1‑136 (MQ=255)
                                                                                                                          catgttGGGCCGCAGCCCCAGCGTGCTGGTCAAAGGGTTGTGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTACCGGCAGt          <  3:158133/132‑1 (MQ=255)
                                                                                                                          cAGCTTGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTACCGGCAGt          <  2:279563/136‑1 (MQ=255)
                                                                                                                              ttGGGCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTACCGGCAGTGAc       >  4:161389/1‑135 (MQ=255)
                                                                                                                                gggCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTACCGGCAGTGACCAg    >  6:317854/1‑136 (MQ=255)
                                                                                                                                gggCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTACCGGCAGTGACCAg    >  6:176001/1‑136 (MQ=255)
                                                                                                                                 ggCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTACCGGCAGTGACCAgg   <  6:168016/136‑1 (MQ=255)
                                                                                                                                  gCCGCAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTACCGGCAGTGACCAGGc  >  4:6445/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                                    
CCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGACAACGAAGGGCATGCCTACCTGGACTTCACTCAGGGCGGTGCGGTCAACAGCTTGGGCCACAGCCCCAGCGTGCTGGTCAAAGCGTTGGGCAACCAGGCTCAGGCGTTGATCAACCCGGGCTCGGGCTTCCACAACCGTGGCTTGCTGGGGCTGGTCAACCGGCTGTGTCAGAGTACCGGCAGTGACCAGGC  >  ACB143/445226‑445491

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