Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,342,221 C→G 10.1% R104G (CGG→GGG)  melA → alpha‑galactosidase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,342,2210CG10.1% 227.1 / 2.4 88R104G (CGG→GGG) melAalpha‑galactosidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (15/61);  new base G (7/3);  total (24/64)
Fisher's exact test for biased strand distribution p-value = 2.20e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

TTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATCCGCATATCAAACAGGTCGGGCTGTGCCATTCGGTGCAGGGAACGGCGGAA  >  NC_000913/4341987‑4342454
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
ttCCATCGCGGGGCGCTGAAAACGGCGCATATTCCCCTGATGGACATTGACCCCACCCGCCGGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCa                                                                                                                                                                                                                                               <  1:1203869/239‑1 (MQ=255)
ttCCACCGGGGGGCGTTGAAACCGGCGCATATGCCCCTGAGGGACTTTGACCCCCCCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGTGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCa                                                                                                                                                                                                                                               <  1:324125/239‑1 (MQ=255)
          aGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGCCGTGGTGGCCTTTCAGATTGGCGGGTATTAACCCTTCCCCGGTGGCTGGTTTCTGGGTCTTTAAGGGCGTGGttggggg                                                                                                                                                                                                                                    >  2:1113984/1‑233 (MQ=255)
             cgcTGAAAACGGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGGTATTAAACTTGGACCGTGGCTGGTTTTGTGGGCTTTTAGGGGCGTgggcgtggacaa                                                                                                                                                                                                                                 >  1:1070249/1‑229 (MQ=255)
                 gAAAGCGGCGCATTTCGCCCTGATGCACATTTGCCCCACCCGCCTGGAAGAGTCGCATATTGGGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCa                                                                                                                                                                                                                             <  1:729301/239‑1 (MQ=255)
                         ggcATATTGCCCTGATGGACATGGCCCCCACCCGCCTGGAAGAGTCGCATATTGGGGTGCGTAAGCTGAGGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGa                                                                                                                                                                                                                     <  1:1113984/237‑1 (MQ=255)
                         caaatttcgCCCTGTGGGCCATGGCCCCCCCCCGCCTGGAAGAGTCGAATTTTGGGGTGCGTAAGCTGAGGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGa                                                                                                                                                                                                                     <  1:755857/230‑1 (MQ=255)
                                  ctctgttgtccttttcccccccccccCGGGAAGGGTCGATTATTGGGGTGGGAAGAGTGATGAGTTCAGCGGGGGCCAGGGGTAAAATCACCTGCCACCCCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTgggg                                                                                                                                                                                                         <  2:178212/223‑1 (MQ=255)
                                  cccTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGGTATTAACCTTGCACGGGGACTGTGTTTTTGGTCTGTAAGGGGCATGGTTTGGaaaccgtctctttttaacatatccg                                                                                                                                                                                                             >  1:852015/1‑215 (MQ=255)
                                  acgggaggcaattggcccccACCCGCCAGGAAGGGTCGCATATTGGGGGTCGTAAGCTGAAGGATTAAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTgggg                                                                                                                                                                                                         <  1:945804/223‑1 (MQ=255)
                                   ccagggtgaattgcaCCCCCCCCGCCTTGAAGGGTCCCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTgggg                                                                                                                                                                                                         <  1:1072343/224‑1 (MQ=255)
                                    cTGATGGACATTGCCCCCCCCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGcc                                                                                                                                                                                                       <  1:116962/239‑1 (MQ=255)
                                      tagggaattggACCCCCCCGGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCgg                                                                                                                                                                                                     <  1:30145/229‑1 (MQ=255)
                                      gtgggcaatggcccccccccGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCgg                                                                                                                                                                                                     <  1:725507/227‑1 (MQ=255)
                                       aTGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTGTTGAAACCTGTCACGGTGGGTGTTTTTGAGGGCTGAAGGGGCACTGGTCTGGAAAAACACATGCCGCTCCATTGGGGGCggg                                                                                                                                                                                                    >  2:495273/1‑239 (MQ=255)
                                       aTGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGGTATGTACCTTGCACGGGGACTGGTTTTTAGGGGTGTAAGGGGCCTGGGGTGGGACAAAACAATGTCCGTACGGTGGGGGCggg                                                                                                                                                                                                    >  1:806602/1‑239 (MQ=255)
                                        accctattccccccccccGCCTGGAAGGTGCGCTTTTTGTGGGGCGAAAGTGGATGGATTCAGCAGGGCCCAGCGGCAAAATCACCTGCCACCCCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCATGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGc                                                                                                                                                                                                   <  2:528337/229‑1 (MQ=255)
                                         ggaactggcccccaccgcgcaggAAGATTCGCATATTGGGGTGGGAAAGGTGAGCGATGTAGGAGGGGCCAGAGGCAAAATCCCCTGCCACCCCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGTTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGc                                                                                                                                                                                                   <  1:1191441/217‑1 (MQ=255)
                                         gcactttcaCCCCACCCGCGTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGCTTCGACAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGc                                                                                                                                                                                                   <  1:195207/230‑1 (MQ=255)
                                         cccctttccccccccccccCTGGAAGGGTCGATTTTTGTGGTGGGTAAGCTGGTGGATTGAGCAGGGGCCGGCGGAAAAATCCCCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCggg                                                                                                                                                                                                    <  2:761895/228‑1 (MQ=255)
                                           acttggcccccccccGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCg                                                                                                                                                                                                  <  1:181204/230‑1 (MQ=255)
                                                    ccccccGCCTGGAAGGGTCGCATATTGTGGTGGGAAAGCGTATGGATTAAGCAGGGGCCAGGGGAAAAATCACCTGGCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgcgc                                                                                                                                                                                       <  1:1366150/236‑1 (MQ=255)
                                                     cACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGGGGCATTTCAGATTGGCGGTGATGAACCTTGTCCGAGGTATGATTTTGGGGTTTGTGATGGGCGAGGGTGGGGAAAAAAAATTTCCCGAACGTTGGGGCCgggggggaggatggggc                                                                                                                                                                                      >  2:181204/1‑228 (MQ=255)
                                                      aCCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATTTACCCTTGCCGGGTACGGAGTTCTTGGGGTGTAAGCAGCATGGTCTGGGACCAACCACTGCCGGACACGTGTGGGCCGGGGGGGATAATgcggg                                                                                                                                                                                      >  2:1090042/1‑237 (MQ=255)
                                                      aCCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGGGGCTTGTTTTGAGGGCTGTAAGCGGCCTGGGCTTGGACAAACCCTTTCCGATACGTTTGGGgcgggcggggattttccct                                                                                                                                                                                      >  1:761895/1‑227 (MQ=255)
                                                       cccGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTTAACGGGGACTGATTTTGCGGGCTTTTAGCGGCGTGGGCTGGGACAAACCAATGCCCATACATTTGGTGCCGGCGGGGTTATgcgctgt                                                                                                                                                                                    >  1:528337/1‑236 (MQ=255)
                                                                     tCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATTAACCTTTCACGGTGACTGGTTTCGAGGGCTTTAATGGGCATGGGCTGGAACAAACCATTTCCCAGACGTTTGGGCCGGGCCGGGTTATGCCCGCCCTCCGGACAATTc                                                                                                                                                                       >  1:416321/1‑238 (MQ=255)
                                                                               tggggCGTAAGCTGATGGATTCAGCAGGGCCCAGCGGCAAAATCCCCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgcgcg                                                                                                                                                                                      <  1:727609/213‑1 (MQ=255)
                                                                                       aaGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATTAAACTTTCTCGGGGCCTGGTTTTTGGGGCTGTAGAGGGGCGGGGTTGGTACAAAACACTTCCGATACAGTTGGGGCGGGGGGGGTTTTGAGGGCGCCCCCTACCAATCCCCACTCGGTGGAAAAttt                                                                                                                                                   >  2:754778/1‑239 (MQ=255)
                                                                                           tGATGGATTCAGCAGGGGCCAGTGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAg                                                                                                                                              <  1:145255/239‑1 (MQ=255)
                                                                                           tGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGa                                                                                                                                                                                                                     <  1:62460/172‑1 (MQ=255)
                                                                                           cgagggtttcggggggggCCAGCGGAAAAATCCCCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCg                                                                                                                                                <  2:1070249/223‑1 (MQ=255)
                                                                                            gagggATTCAGGAGGGGCCAGCGGCAAAACCACCTGCCAACCCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTGGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAg                                                                                                                                              <  1:540842/235‑1 (MQ=255)
                                                                                              tGGATTCAGCAGGGGCCAGGGGAAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGa                                                                                                                                            <  1:484740/238‑1 (MQ=255)
                                                                                              tGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGGGACTGATTTTGAGGTCTGTAAGCGGCACGGGCTGGAACCAACCAATGCCGATACGTTGGGGCCGGGCGGGATTTTTCGCGCGCTCCGTACCAATCCCCCTCTGTgggtcatttttgagggc                                                                                                                                           >  1:483398/1‑224 (MQ=255)
                                                                                              gggATTCAGCAGGGGCCAGCGGAAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGa                                                                                                                                            <  1:1151808/237‑1 (MQ=255)
                                                                                                    gaggagGGGCCGGGGGCAAAATCACCTGCCACCCCCAACAGAAAGAAGGCTTAGAGGATGCCGATTTTGGCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACg                                                                                                                                     <  1:1291860/235‑1 (MQ=255)
                                                                                                    cgccgggggCCAGGGGCAAAATCCCTTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACg                                                                                                                                     <  2:380719/234‑1 (MQ=255)
                                                                                                    cagcagGGTCCAGCGGAAAAATCCCCTGCCAATCCCAACAAAAATAAGCCTTAGAGGACTCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTTGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACg                                                                                                                                     <  1:235013/239‑1 (MQ=255)
                                                                                                    cagcagGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATTAAACTTGCACGGTGACTGATTTTGAGGGCCGTAATCGGCATGGTCTGGAACAAACCATTGCCGATAAGTTGGGGGCGGGCGGGGTTTTTCGCGCCCTACGCACCAATCCCCACCTTTTGGAAAAATTGTGGGACATGACg                                                                                                                                     >  1:438944/1‑239 (MQ=255)
                                                                                                                       aaTCCCCTGCCACACCCAAAAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAATCCATTGCCGATCCGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGcca                                                                                                                  <  1:1033267/239‑1 (MQ=255)
                                                                                                                         tCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGccacc                                                                                                                <  1:96608/239‑1 (MQ=255)
                                                                                                                            cctgcccaccccAACAGAAAGAGGCCTTAGAGGATGCCGATTTTGTCGTGGGGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATg                                                                                                             <  2:1096243/230‑1 (MQ=255)
                                                                                                                              tGCCACACCCAACAGGAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCt                                                                                                           <  1:880559/239‑1 (MQ=255)
                                                                                                                                ccACACCCAACAGAAAGAATCCTTAGAGGATGCCGTTTTTGTCGGGGTGGCATTTCAGATTTGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCa                                                                                                         <  1:535344/239‑1 (MQ=255)
                                                                                                                                 cacaCCCAACAGAAAGGAGCCTTTGAGGATGCCGTTTTTGGCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCaa                                                                                                        <  2:416321/239‑1 (MQ=255)
                                                                                                                                  ccccccAAAAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGGGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAAc                                                                                                       <  2:1303690/236‑1 (MQ=255)
                                                                                                                                   cccccAACAGAAAGGAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGACCCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAGCCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGcca                                                                                                                  <  1:1013918/225‑1 (MQ=255)
                                                                                                                                       cAACAGAAAGAAGCCTTAGGGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGt                                                                                                  <  1:1356245/239‑1 (MQ=255)
                                                                                                                                       cAACAGAAAGAAGCCTTAGAGGATGCCGATTGTGTAGTGGGCGCATTTCAGAATGGCGGTTATGAACCTTGCACGGTGACTGAGTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAATCCATTGCCGATACGTTGGGGCCGGGCTGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCTAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGt                                                                                                  <  1:1286999/239‑1 (MQ=255)
                                                                                                                                        aaCAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGcca                                                                                                                  <  1:1121284/222‑1 (MQ=255)
                                                                                                                                              aagaagCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTGCACGGGGGCTGGGTTTGAGGGCTGTTAGAGGCCTGGTCTGGGAGCAAACACTTCCCAGACGTTTGGGGCGGGGCGGTGTTTGCGCGGGCTGCGTCCCATTCCGCCTCGTCGGCGAATTTGTGGGGGCACGACGGAAGTAGGCCCCGCCGccaccacccccaaataaagtaaaaca                                                                                           >  2:573389/1‑219 (MQ=255)
                                                                                                                                               agaagCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCaa                                                                                         <  1:365639/239‑1 (MQ=255)
                                                                                                                                                  agCCTTGGATGTTGCCGTTTTTGTCGTGGTGGCATTTCAGATTGGCGGTTTTGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATgg                                                                                      <  2:1315480/239‑1 (MQ=255)
                                                                                                                                                  agCCTTAGAGGTTGCCGTTTTTGTCGGTGTGGCATTTCGAGTTGGCGGGTATGAACCTTGCACGGTGACTGACTTGGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCACGCTACGTTCCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATgg                                                                                      <  1:147277/239‑1 (MQ=255)
                                                                                                                                                   gtccttgagGGTGCCGTTTTTGTCGGGGTGTCAGTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATgg                                                                                      <  2:1286266/232‑1 (MQ=255)
                                                                                                                                                   gccttgtgggtggcccgttttGTCGTGGGGGGATTTGAGATTGGCGGTTTTGAACCTTGCCCGGTTACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCCGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATg                                                                                       <  2:545615/220‑1 (MQ=255)
                                                                                                                                                    ccTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACa                                                                                                                                          <  1:115430/186‑1 (MQ=255)
                                                                                                                                                       tAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTGCACGGGGACTGATTTCGAGGTCTGTAAGCGGCATGGGCTTGAACAAAACATTGCCGGTACGTTTGGGCCGGGCGGTATTATGCGCGCGCCACGTACCATTCCGCCTCTGTTGCAAAATTTCGGGGACATGACGGGAGTCTGCCCCCGTGCCACCCTTCCCAACTATGTTAAACCCAAGGCCGag                                                                                  >  1:324463/1‑238 (MQ=255)
                                                                                                                                                         gagGATGCCGATTTTGGCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGaa                                                                               <  1:754778/239‑1 (MQ=255)
                                                                                                                                                          agGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAg                                                                                                                                              <  1:1270608/176‑1 (MQ=255)
                                                                                                                                                              ttccgttttgtgcttggtggATTTTCGGATTGGGGGTTATGAACCTTGCACGGGCACGGATTTAGAGGTCGGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATAcc                                                                           <  2:806602/224‑1 (MQ=255)
                                                                                                                                                                    ttttGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCggg                                                                                                                                                                                                    <  1:24020/114‑1 (MQ=255)
                                                                                                                                                                    ttttGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTggg                                                                       <  1:1103413/236‑1 (MQ=255)
                                                                                                                                                                       tGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGccacc                                                                                                                <  1:338076/193‑1 (MQ=255)
                                                                                                                                                                                    tttAAGTTTGGGGGTTTTGGTCCTTGCACGGTTGTTGATTTGGGGGGTTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATCCGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATcc                                                   <  2:282867/239‑1 (MQ=255)
                                                                                                                                                                                                    aTGAACCTTGCACGGTGTCTGATTTCGAGGTCTGTAAGCGTCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCAGGGCGTTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATCCGCATATCAAACAGGt                                    <  1:1174815/238‑1 (MQ=255)
                                                                                                                                                                                                      tttgcaccttgcccgttgtattgtttCGGGGTCTTGAAGGGGCATGGTTTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATCCGCATATCAAACAGGt                                    <  2:1219065/216‑1 (MQ=255)
                                                                                                                                                                                                          gactttGCACGGTTACTGTTTTCGGGGTATGTAATGGGCATGGTCTGGAACAAACCATTGCCGATCCGTTGGGGCCGGGCGGTATTATGCGCGCGTTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATCCGCATATCAAACAGGTCggg                                <  2:438944/235‑1 (MQ=255)
                                                                                                                                                                                                          aCCTTGCATGGTGACTGATTTCGGGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATCCGCATATCAAACAGGTCggg                                <  1:8312/238‑1 (MQ=255)
                                                                                                                                                                                                          aCCTTGCACGGTGACTGATTGCGAGGTCTGTAAGGGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATCCGCATATCAAACAGGTCggg                                <  1:221480/238‑1 (MQ=255)
                                                                                                                                                                                                             cTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGGGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCACATTTGCGCGGACATGACGGAAGTCTGCCCCGATGCCAACATGCTCAACTATGTAAACCCCATGGAGATGAAAACCTGGCGCAGTCTTGCCCCCTATCCGCAAACAAAAACGGTCGGGcgg                             >  2:558866/1‑237 (MQ=255)
                                                                                                                                                                                                               tGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTGCGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTTCCCCGATGCCACCATCCTCAACTATGTTAACCCAATGGCGATGAATACCTGGCGGCTGAATGCCCGCTATCCGCATATCAAACAGGACgggggggg                           >  2:1213810/1‑234 (MQ=255)
                                                                                                                                                                                                               tGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACa                                                                                                                                          <  1:381952/130‑1 (MQ=255)
                                                                                                                                                                                                               tGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACa                                                                                                                                          >  2:381952/1‑130 (MQ=255)
                                                                                                                                                                                                               tGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCa                                                                                          >  2:962903/1‑177 (MQ=255)
                                                                                                                                                                                                               tGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCa                                                                                          <  1:962903/177‑1 (MQ=255)
                                                                                                                                                                                                               tGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGATTACCTGGGCGATGTATGCCCGCTATCCGCCTATCCACAAGGTCGGCCtgt                            >  1:1042697/1‑238 (MQ=255)
                                                                                                                                                                                                                 cACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGCTGAATACCTGGGCGATGAATGCCCGCTATCCCAATATCAAAAAGGTGGGGGTGGGca                         >  2:403364/1‑238 (MQ=255)
                                                                                                                                                                                                                    ggTGACTGTTTTTGAGGGCTGTAAGCGGCATGGTTGGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTAACATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATCCGCATATCAAACAGGTCGGGCTGTGCCAtt                      <  2:324463/239‑1 (MQ=255)
                                                                                                                                                                                                                    ggTGACTGATTTCGAGGTCTGGAAGCGGCATGGTCTGGAAAAAACCATTGCCGATACGTTGGGGCCTGGCGGTATTAATCGCGCGCTTCGTACCATTCCGCATCTGTGGCAAATTGGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATCCGCAAATCAAACAGGTCGGGCTGTGCCAtt                      <  1:667706/239‑1 (MQ=255)
                                                                                                                                                                                                                     gggACGGATGTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAGCCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATCCGCATATCAAACAGGTCGGGCTGTGCCAtt                      <  1:545612/236‑1 (MQ=255)
                                                                                                                                                                                                                      tGACTGATTTCGAGGTCTGTAAGCGGCATGGTATGGAACAAACCATTGCCGATTCGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGTCATGACGTGAGTCTTCCCCAATGCCACCATGCTCAACTATGTTAACCCAAAGGCGATGAATACCTGTGCGATGTATGCCCGCTATCCCCACATCAAAACGGTCGGGcgtggccaatcg                    >  2:667706/1‑228 (MQ=255)
                                                                                                                                                                                                                                ggAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATCCGCATATCAAACAGGTCGGGCTGTGCCATTCGGTGCAGGGaa          <  1:495273/238‑1 (MQ=255)
                                                                                                                                                                                                                                cGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATCCGCATATCAAACAGGTCGGGCTGTGCCATTCGGTGCAGGGaa          <  1:403364/239‑1 (MQ=255)
                                                                                                                                                                                                                                  aGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTTTGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATCCGCATATCAAACAGGTCGGGCTGTGCCATTCGGTGCAGGGAAcg        >  1:485893/1‑239 (MQ=255)
                                                                                                                                                                                                                                       tGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGGCATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGAGGTATGCCCGCTATCCGCATATCAAACAGGTCGGGCGGTGCAATTCGGTGCAGGGAAGGGGGGa   >  1:646140/1‑239 (MQ=255)
                                                                                                                                                                                                                                         gaaGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATTCCTGGGCGATGTATGCCCGCTATCCGCATATCAAACAGGTCGGGCTGTGCCATTCGGTGCAGGGAACGGCGGaa  <  1:1249289/237‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
TTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCCCCGATGCCACCATGCTCAACTATGTTAACCCAATGGCGATGAATACCTGGGCGATGTATGCCCGCTATCCGCATATCAAACAGGTCGGGCTGTGCCATTCGGTGCAGGGAACGGCGGAA  >  NC_000913/4341987‑4342454

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: