Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,341,972 G→A 8.9% P242L (CCG→CTG)  tsdA ← Thiosulfate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,341,9720GA 8.9% 208.8 / 5.7 68P242L (CCG→CTG) tsdAThiosulfate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base A (3/3);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCCCATCATTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCGGCGGCATTACATACTCGCCTGCGACCTTTTGCCCTTCAC  >  ACB143/4341839‑4342107
                                                                                                                                     |                                                                                                                                       
aaCCCCATCATTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCCCGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGc                                                                                                                                       >  3:349189/1‑136 (MQ=255)
  ccccATCATTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACGAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCtt                                                                                                                                     <  5:92988/136‑1 (MQ=255)
  ccccATCATTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCtt                                                                                                                                     <  5:304348/136‑1 (MQ=255)
   cccATCATTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGGGGTTGTTACACCCCGCCACGTCCCAGGCCAGCGGATAGCTGAGGCTGCCCGGCTTg                                                                                                                                    >  1:179299/1‑136 (MQ=255)
   cccATCATTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTg                                                                                                                                    >  1:108796/1‑136 (MQ=255)
   cccATCATTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTg                                                                                                                                    >  7:254445/1‑136 (MQ=255)
   cccATCATTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTg                                                                                                                                    >  5:392571/1‑136 (MQ=255)
    ccATCATTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGc                                                                                                                                   <  6:384768/136‑1 (MQ=255)
     catcatTGGCATGAAACTTCCCCCGCGTCTTCTCGATGGATCCCTCTCCCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTGGCTCCCCGCCGCCACGGCCCCGGCCTGCTGAGCCCTGCGGCTGCCCAGCTTGcc                                                                                                                                  >  8:59094/1‑136 (MQ=255)
     catcatTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCCCGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGcc                                                                                                                                  >  2:322403/1‑136 (MQ=255)
     catcatTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATACACGCCGCCAAGTCACATGCCTGGTGATAGCTGAGGCTGCCCAGCTTGcc                                                                                                                                  >  8:63917/1‑136 (MQ=255)
     catcatTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCCTGGCGGTTGAACCACCCCGCCACCTCCCACGCCTGCAGATACCTGAGGCGGCCCACCTTGcc                                                                                                                                  >  8:63916/1‑136 (MQ=255)
       tcatTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCa                                                                                                                                <  7:370323/136‑1 (MQ=255)
        catTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCCTGGCGGTTGATCCACGCCGCCCCGTCCCCGGCCTGCTGCTCGCGGAGGCTGCCCGGCTTGCCCAg                                                                                                                               >  6:357503/1‑136 (MQ=255)
         atTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGc                                                                                                                              <  2:210621/136‑1 (MQ=255)
         atTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCGCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGc                                                                                                                              <  2:370389/136‑1 (MQ=255)
          ttGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCg                                                                                                                             >  3:118657/1‑136 (MQ=255)
               aTGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCCCGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCAtg                                                                                                                        >  2:324213/1‑136 (MQ=255)
               aTGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATACCCGCCGCCACCTCCCAGGCCAGCTGATCGCGGAGGCTGCCCGGCTTGCCCAGCGGCAtg                                                                                                                        >  4:220785/1‑136 (MQ=255)
               aTGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCGAGTGGCGGGTGGTACCCGCCGCCCCGGCCCAGGCCTGCTGTTCGCGGAGGCTGCCCGGCTTGCCCCGCGGCAtg                                                                                                                        >  7:136545/1‑136 (MQ=255)
                               gTCTTCTCGATGGATCCCTCTACAAACCGCGGGTCCTGTGGCCGCTCGGGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAacg                                                                                                        <  8:141548/136‑1 (MQ=255)
                                 cttctCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAacgac                                                                                                      >  3:224175/1‑136 (MQ=255)
                                  ttctCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAacgacg                                                                                                     <  8:166990/136‑1 (MQ=255)
                                   tctcGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAacgacg                                                                                                     <  3:212233/135‑1 (MQ=255)
                                    ctcGATTGATCCCTATACAACCCGCGGTTCCTGTGGCCGCTCGTGGCGTTTGATCCACGCCGCCACGTCCAAGGCCTGCTGAGCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTCATGAacgacg                                                                                                     <  4:65313/134‑1 (MQ=255)
                                    ctcGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAacgacg                                                                                                     <  8:293436/134‑1 (MQ=255)
                                    ctcGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAacgacg                                                                                                     <  4:258330/134‑1 (MQ=255)
                                    ctcGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAacgacg                                                                                                     <  5:123522/134‑1 (MQ=255)
                                    ctcGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAacgacg                                                                                                     <  7:306221/134‑1 (MQ=255)
                                      cGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGAcgccgc                                                                                                 >  3:369370/1‑136 (MQ=255)
                                        aTGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCgg                                                                                               >  5:339170/1‑136 (MQ=255)
                                             tCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTg                                                                                          <  4:109372/136‑1 (MQ=255)
                                             tCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTg                                                                                          <  7:22013/136‑1 (MQ=255)
                                             tCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTg                                                                                          <  2:307677/136‑1 (MQ=255)
                                             tCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTg                                                                                          <  3:284602/136‑1 (MQ=255)
                                             tCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTg                                                                                          <  2:72348/136‑1 (MQ=255)
                                             tCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTg                                                                                          <  3:291214/136‑1 (MQ=255)
                                               ccTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGAt                                                                                        >  8:277815/1‑136 (MQ=255)
                                                    aCCAACCGCGGGTCCTGTGGCCGCTGGTGGCGTTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGt                                                                                   <  6:50721/136‑1 (MQ=255)
                                                          cgcgGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGc                                                                             <  4:246336/136‑1 (MQ=255)
                                                           gcgGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGcc                                                                            >  2:28079/1‑136 (MQ=255)
                                                           gcgGGTCCTGTGGCCGCTCCTTGCCCTTGAAACACCCCGCCACGTCCCACGCCAGCTGAACGCTGAGGCTGCCCCGCTTTCCCAGCGGCAATTTTTAAATTATGAAACACCCCCCCCTAATGATTCGGTGCATGc                                                                             >  2:163786/1‑135 (MQ=255)
                                                             gggTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAg                                                                          >  4:197673/1‑136 (MQ=255)
                                                               gTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGcc                                                                        <  8:16544/136‑1 (MQ=255)
                                                                  cTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGccccc                                                                     <  7:384535/136‑1 (MQ=255)
                                                                  cTGTGGCCGCTCGTGGCGGTGGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGccccc                                                                     <  1:127513/136‑1 (MQ=255)
                                                                     tgGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGcccccc                                                                    <  3:287848/134‑1 (MQ=255)
                                                                        ccGCTCGTGGTGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCGGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTg                                                               <  5:292417/136‑1 (MQ=255)
                                                                          gCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTTCTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCCGCCCCCCAGTTGTa                                                             >  2:318880/1‑136 (MQ=255)
                                                                            tCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAg                                                            >  2:34519/1‑135 (MQ=255)
                                                                                ggCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATc                                                       >  7:35770/1‑136 (MQ=255)
                                                                                              gccgccACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAgcg                                         <  5:374127/136‑1 (MQ=255)
                                                                                              gccgccACGTCCCAGGCCTGCTGATCGATGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAgcg                                         <  6:129286/136‑1 (MQ=255)
                                                                                                 gccACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAgcggcg                                      >  6:130826/1‑136 (MQ=255)
                                                                                                   cACGTCCCAGGCCTGCTGTTGGCTGAGGCTGCCCAGCTTGCGGAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCggcggc                                    <  7:63916/136‑1 (MQ=255)
                                                                                                   cACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCAGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCggcggc                                    <  7:63917/136‑1 (MQ=255)
                                                                                                   cACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCAGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGAGGCGGGGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAACGCggcggc                                    <  7:59094/136‑1 (MQ=255)
                                                                                                     cGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTg                                                                                          >  5:221804/1‑80 (MQ=255)
                                                                                                     cGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTg                                                                                          <  6:221804/80‑1 (MQ=255)
                                                                                                      gTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCGGCGGCAt                                  >  8:126806/1‑135 (MQ=255)
                                                                                                          cAGGCCTGCTGATCGCGTAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCGGCGGCATTACAt                             <  8:136545/136‑1 (MQ=255)
                                                                                                           aGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACTACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCTGCGGCATTACATa                            >  4:202611/1‑136 (MQ=255)
                                                                                                                    gATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCGGCGGCATTACATACTCGCCTGc                   >  5:110626/1‑136 (MQ=255)
                                                                                                                    gATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCGGCGGCATTACATACTCGCCTGc                   >  3:26531/1‑136 (MQ=255)
                                                                                                                     aTCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCGGCGGCATTACATACTCGCCTGc                   >  3:132775/1‑135 (MQ=255)
                                                                                                                     aTCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCGGCGGCATTACATACTCGCCTGCg                  >  6:210725/1‑136 (MQ=255)
                                                                                                                                 gCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCGGCGGCATTACATACTCGCCTGCGACCTTTTGCCCt      <  5:310455/136‑1 (MQ=255)
                                                                                                                                  cccGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCGGCGGCATTACATACTCGCCTGCGACCTTTTGCCCt      <  1:129578/135‑1 (MQ=255)
                                                                                                                                     ggCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCGGCGGCATTACATACTCGCCTGCGACCTTTTGCCCTTCAc  >  6:361693/1‑136 (MQ=255)
                                                                                                                                     ggCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCACAAAGCGGCGGCATTACATACTCGCCTGCGACCTTTTGCCCTTCAc  >  8:340854/1‑136 (MQ=255)
                                                                                                                                     ggCTTGCCCAGCGGCATGTTGTACTTGATGAAAGACGCCGAGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCACAAAGCGGCGGCATTACATACTCGCCTGACACCTTTTGCCATTCAc  >  8:340850/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                                       
AACCCCATCATTGGCATGAAACTTCACCCGCGTCTTCTCGATGGATCCCTCTACCAACCGCGGGTCCTGTGGCCGCTCGTGGCGGTTGATCCACGCCGCCACGTCCCAGGCCTGCTGATCGCTGAGGCTGCCCGGCTTGCCCAGCGGCATGTTGTACTTGATGAACGACGCCGCGGTATTGATGCGGTGCATGCCAGCCCCCCAGTTGTAGGAATCCTTGCCCCAAAGCGGCGGCATTACATACTCGCCTGCGACCTTTTGCCCTTCAC  >  ACB143/4341839‑4342107

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: