Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,877,673 A→G 12.6% intergenic (+57/+894) iap → / ← ygbF aminopeptidase in alkaline phosphatase isozyme conversion/CRISPR adaptation ssRNA endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,877,6730AG12.6% 217.2 / 2.2 95intergenic (+57/+894)iap/ygbFaminopeptidase in alkaline phosphatase isozyme conversion/CRISPR adaptation ssRNA endonuclease
Reads supporting (aligned to +/- strand):  ref base A (26/53);  new base G (11/2);  total (37/58)
Fisher's exact test for biased strand distribution p-value = 6.16e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAA  >  NC_000913/2877447‑2877907
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
aGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCatcatc                                                                                                                                                                                                                                             >  1:179064/1‑238 (MQ=255)
    aCGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGACTCAAAGTGGAATCAATCATCGCAt                                                                                                                                                                                                                                         >  1:603691/1‑238 (MQ=255)
     cgcgCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCAAAGGCGTCTTGATGGGGTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGAAACCAAGAGAAAAGCATCAACGctacc                                                                                                                                                                                                                                       >  2:647787/1‑235 (MQ=255)
           aaaaCCCCTGCCTGCCGGGCGGGAAATAGCTGGCAGAGCGGAAGACTGGACAATCCCCAACATATTGATAAGGCTCTTCCTGGACGAATAGAACGTCGCTGCCGGGACGTTATGCGGATAATGCTACCTCTGGGGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTgcg                                                                                                                                                                                                                                 <  2:451960/239‑1 (MQ=255)
            aaaCCCCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTgcgc                                                                                                                                                                                                                               <  1:1329629/239‑1 (MQ=255)
            aaaCACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAAAGGGGGGGGGGATCCAAAGGAGATCTAACATCACACCCAgggagc                                                                                                                                                                                                                                >  2:1329629/1‑234 (MQ=255)
             actcactgctctccgcgggggAAATAGCAGGCTTGAGGGAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGAAAGATTGGAACGTGGATGCCGGGGCGTTATGCGGAGAATGATACCTCTGGGGGAGTGGTTGGCGAAGGCGTTTTGATGGGTTTGAAAATGGGCGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCATGTGATATCCATCATCGCATCCAGTGCGcc                                                                                                                                                                                                                              <  2:202866/221‑1 (MQ=255)
             acccccTGCCTTCCCGGGGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGcc                                                                                                                                                                                                                              <  2:1306777/235‑1 (MQ=255)
             aaCCCCTGCCTTCCCGGGGGGAAATAGCTGGAATGACGAAAGCCTGGATAATCCCCAACATATTGAAAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGcc                                                                                                                                                                                                                              <  2:801056/239‑1 (MQ=255)
                          ccGGGGGGAAATATCTGGCATAACGTAAGACTGGATAATCCCCAAAATATTGATAAGGCTCTTCCTGGAAGAATAGAGCGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCg                                                                                                                                                                                                                 <  2:345186/239‑1 (MQ=255)
                                 gAAATGGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATg                                                                                                                                                                                                          <  1:1350406/238‑1 (MQ=255)
                                 gAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATg                                                                                                                                                                                                          <  1:764096/238‑1 (MQ=255)
                                 gAAATAGCTGGAATGACGTAAGGCTGGATAATCCCCAACATATGGATAAGGCTCTTCCTGGAAAAATAGAACGTCGCTGCCGGGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATg                                                                                                                                                                                                          <  2:459206/238‑1 (MQ=255)
                                   aaTAGCTGGCAGGGCGTAAGCCTGGATAATCCACAACATATTGATAAGGCTCTTCCGGGAAGAATAGAACGTCGGTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCg                                                                                                                                                                                                        <  1:539041/238‑1 (MQ=255)
                                       gctgtctggACGGAAGGCGGGATAAGCCCCAACAAATTGAGAAGGCTCTTCCGGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCgggg                                                                                                                                                                                                     <  1:1214976/229‑1 (MQ=255)
                                       gCTGGCTTGCCGTAAGACTGGATTATCCCCAACATTTTGATAAGGCTTTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGaa                                                                                                                                                                                                   <  2:573150/239‑1 (MQ=255)
                                        ggtgggtataagacgcaggtgtAATCCCCAACCTATTGAAAAGGCTCTTCCGGGAAGAATAGTACGGCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGaa                                                                                                                                                                                                   <  1:1360747/217‑1 (MQ=255)
                                         tGGGATGTCTTAAGACTTGATAACCCCCAACATATTGAAAAGGCTCTTCCTGGAAGAATAGAACGTCTCTGCCGTGACGTTAGGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAaca                                                                                                                                                                                                 <  2:462900/239‑1 (MQ=255)
                                          ggCATGACGTAAGACTGGATAATCCCCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAacac                                                                                                                                                                                                <  1:1160460/239‑1 (MQ=255)
                                           gCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGACTCCAAGTGGTATCCATCATCGCCTCCGATGCGCCCCGTTTATCCCCGCCGATGCgggggaaccc                                                                                                                                                                                               >  1:378326/1‑233 (MQ=255)
                                            cATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTGGCTGCCGGGACGTTATGCGGATAATGCAACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGTGACTGGGAGTACTACCGCAGCAGCGGGGGAAGTGCAAAGGGGAGACATCACCGACCCCACGGCGCCCGGATTATcccccccgaacagcggaaccaca                                                                                                                                                                                               >  2:951850/1‑220 (MQ=255)
                                                 cGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGACCGGCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGTTTTGAAAATGGTAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATg                                                                                                                                                                                                          <  1:99618/222‑1 (MQ=255)
                                                     aGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGGGATATCCACCACCCCATCCCGTGGGCCCGGTTTATCCCCCCCCGTGCgggggaaacaaacgaccggg                                                                                                                                                                                     >  2:578867/1‑223 (MQ=255)
                                                           gATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGAATCCAAAGGGTATTCACCCTCCCATCCCGGGCGCCCGGGGTATCCCCGCTGAGGCGGGGaaaaacacccgcagagggggaa                                                                                                                                                                               >  2:19627/1‑219 (MQ=255)
                                                                tCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCACAGGATATTCATCATCCCCTTCCGGGGGCCCCGGTTTTCCCCCCGTGAGGGGGGGAAAACCacggccgggggggaaaacc                                                                                                                                                                            >  2:1242897/1‑220 (MQ=255)
                                                                        atatTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCg                                                                                                                                                                                                                 <  1:1042887/193‑1 (MQ=255)
                                                                          ttttaTAAGGCTTTTCCTGGAAGAATAGACTGTCGGTGGCGGGACGTTATTGGGATAATGCTACCTTGGGTGAAGGAGTTGGGGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGCTGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGc                                                                                                                                                          <  2:830150/235‑1 (MQ=255)
                                                                          atTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGc                                                                                                                                                          <  1:316535/239‑1 (MQ=255)
                                                                          atTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAATGATATCCATCATCGCATCCCGTGCGCCCGGTTTATCCCCGCTGGTGCGGGGGACACCAGCGTCAGGCGTGAAATATCCCCCACGTTg                                                                                                                                                            >  1:1015091/1‑238 (MQ=255)
                                                                              aTAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGACTCCAACGGGAATCCCTCATTGCATCCCATGCGCCCCGGTTTTCCCCCCCGGTGCGGGGGAAACCACCGCCCAGGCGCGaaacccacccccctgtccggg                                                                                                                                                         >  2:160288/1‑221 (MQ=255)
                                                                                  ggCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGACTTCAAAGGAAATCCATCATCCCCTCCAGTGCGCCCCGGTTATCCCCGCCGATGGGGGGGAAACCAACCGCACGCGGGAAAtccccacccccttggcggggta                                                                                                                                                    >  2:26099/1‑218 (MQ=255)
                                                                                  ggCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGACTCCAAGTGGTATCCATCATCGCATCCAGTGCGCCCGGGTTATCCCCGCTGATGCGGGGAACACCAGCGGCAAGCGTGAAATCTCACCGCCGTTGCCGGTTTa                                                                                                                                                 >  1:156530/1‑238 (MQ=255)
                                                                                  ggCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCCTCAGGCGTGGAATCTCACCCCCGTTGCCGGGTTa                                                                                                                                                 >  1:265501/1‑238 (MQ=255)
                                                                                        tCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGACTCCAAATGATATCCCTCATCGCATCCAGTGCGCCCGGTTTTTCCCCACTCGTGCGGGGAAAACCACCGGCGCGGGGGGAACTATACCCCCGGTGCCGGTTTATcccccc                                                                                                                                         >  2:1227613/1‑236 (MQ=255)
                                                                                            ggAAAAATAGAACGGCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGgc                                                                                                                                     <  1:1227613/239‑1 (MQ=255)
                                                                                                 aaTAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGTGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCca                                                                                                                                                                                                                                                   <  1:760433/136‑1 (MQ=255)
                                                                                                 aaTAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCg                                                                                                                                                                                <  1:287568/196‑1 (MQ=255)
                                                                                                   tAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGAATCCAAGAGGTATCCATCATCGCATCCACTGCGGCCCGCTTATTCCCCCTGCTGCGGGGAGAACCACCGGCCCGGGGGGAAACTCACCCTCGCTTGCGGtttttacccccccggccggggg                                                                                                                               >  2:103059/1‑220 (MQ=255)
                                                                                                   tAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGGTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGGGGTAAATCTCACCGTCGTTTCCCGTTTATtccctccggccgggggt                                                                                                                              >  1:1248454/1‑225 (MQ=255)
                                                                                                   tAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGGGCGCCCGGTTTATCCCCGCCGGTGCGGGGAACACCAGACGCCGGCGGGAAATTTCCACGCCGTTGCCGGTTTATCCCCCCcggggcggggg                                                                                                                             >  2:416179/1‑238 (MQ=255)
                                                                                                      aaCGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCgggg                                                                                                                              <  1:302296/235‑1 (MQ=255)
                                                                                                        cgtccgtGCCGTGACGTTTTGGCGGTTATGCTTCCTCGGGTGAAGGAGTTGGCGAAGGAGTTTTGGTGGGTTTGAAAATGGGAGCTGGGAGTGCTACCGCAGGGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAActct                                                                                                                        <  2:378326/233‑1 (MQ=255)
                                                                                                        cggcGCTGCCGTGACGTTATGCGGATAATGCTACCTCGGGTGAAGGAGTTGGCAAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAActct                                                                                                                        <  1:647787/236‑1 (MQ=255)
                                                                                                        cGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGAATCCCAATGAAATCCATCATCGCTCCCACAGCGCCCGGTTTATCCCCCCCGTTGCGGGGGACACCCCCGTCAGGCGGGGAAATCCCCCCCACTTGCCGGGTTATccccctcggcccgggggaccc                                                                                                                          >  2:590509/1‑220 (MQ=255)
                                                                                                           cGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGACTCCAAGTGGTATCCCACACCGCAGCCCCGGCGCCCCGTTTATTCCCCCCGATGGGGGGGAACCCCCGCGCAGGCGCGGAATTCCACCGCCGCTGCCGGGTTTACCCCGCCTGCCGGGGGGAAccccgg                                                                                                                      >  2:143536/1‑234 (MQ=255)
                                                                                                           cGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCCAGTGGTATCCATCATACCATCCCGTGCGCCCCGTTTTTTCACCCCGACTCGGGGGACCCCCCAGCCAGGCGGGGAAATCCACCCTCCTTGCCGGGTTTTTCCTCCTGGCGCGGGGGACACTAgg                                                                                                                     >  2:1035792/1‑239 (MQ=255)
                                                                                                           cGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGGGATTCATCATCGCATCCCAGGCGCCCGGGTTATCCCCGCTGACTCGGGGAACACCAACACCAGGCGGCGAAATCCCCCCTCGTTGCCGGTTTTACCCCCGCGGGCGGGGGaaaacccgg                                                                                                                      >  2:68066/1‑233 (MQ=255)
                                                                                                           cGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGGGATATCCCTCATCCCCTCCCGTGCCCCCGGGTTATCCCCCCTGAGACGGGGACAACCCGCGTCAGGCGGGGAATATCACCCCCGTTTCCGGTTTTTCCCCGCGGGTGCgggggcatcccgg                                                                                                                     >  2:1040314/1‑230 (MQ=255)
                                                                                                                ccGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGACTCCAAATGATATTCCTCATCGCATCCAGTGCGCCCCGGTTATCCCCGCCGATGCGGGGAAAACCAGCGTCAGGGGGGAAATTCTCCCCGCGTTGCCGGGTTTTCCCCCCGTGCGCgggggaacctccggtttg                                                                                                                >  2:302296/1‑225 (MQ=255)
                                                                                                                 cgggaaggtttGCGGTTAATTCTACCTCTGGTGAAGGGTGTGGCGAAGGCGGCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCGACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCCGCTGGCGCGGGGAACTCTCGGTTCAgg                                                                                                               <  2:200676/229‑1 (MQ=255)
                                                                                                                 cGTGACCTTATGGGGTTTATGCTTCCTCTGGTGAAGGAGGTGGCGAAGGCGGCTTGGGGGGTCTGAAAATGGGAGCTGGGAGTTCTACCGGAGGGGCGGGGGAACCCCAAGTGATATCCATCATCGCCTCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAgg                                                                                                               <  2:602164/239‑1 (MQ=255)
                                                                                                                   tGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGACTCCAAATGATATCCATCATCACATCCCGTGCGCCCCTGTTTTTCCCCCCGCTGTGGAGAAAACCAGCACCAAGCGGGAAATATCCCCGCAGTTTCCGTTTTTTCCCCTTGTGGCGGGGGGAAccccgggtcgggga                                                                                                              >  2:923546/1‑226 (MQ=255)
                                                                                                                         tttggggATAATGCTTCCTGTGGTGAAGGAGGTGGCGAAGGCGGCTGGATGGGGTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCa                                                                                                        <  2:394110/233‑1 (MQ=255)
                                                                                                                         tATTCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGATTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCa                                                                                                        <  1:923546/238‑1 (MQ=255)
                                                                                                                         tATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCaa                                                                                                       <  1:416179/239‑1 (MQ=255)
                                                                                                                          aTGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGAACCCAAGTGATATCCATCATCGCATCCAGTGCCCCCGGTTTTTCCCCGCTGATGCGGGGAACACCCCCGGCCGGCGTGAAAATCCACCCTCGTTGCCCGGTTTTCCCCTCTTGCCGCGGGGACCCCCGGGTTCGGGCGtggaaa                                                                                                       >  2:1151015/1‑234 (MQ=255)
                                                                                                                          aTGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTTCCGGTTTTTTCCTGCTGGCGCGGGGGACTCTCGGTTAAGGCGTTGCaa                                                                                                       >  1:845483/1‑238 (MQ=255)
                                                                                                                              ggATAATGCTACCTTTGGTGAAGGAGTTGGCGAAGGGTTCTTGAGGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCt                                                                                                   <  2:179064/238‑1 (MQ=255)
                                                                                                                               gATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGc                                                                                                                                                          <  1:728219/186‑1 (MQ=255)
                                                                                                                                aTAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGACTCCAAGTGATATCCCTCATCGCCTCCAGGTCGCCCGCGTTTTCCCCCCTGATCCGGGGGACACCACGCTCCGGCGGTAAACTTTCCCCCCCTTGCCGGTTTATTCCCTCTGCCGCGGGGGACTCCCCGCTCCGGCGGTGGAAACCCgg                                                                                                 >  2:264968/1‑238 (MQ=255)
                                                                                                                                aTAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGACTCCAAGTGATATCCATCATCGCATCCCGTGCCCCCCGTTTATCCCCGCTGGTTCGGGGAACACCAGCGCCCGGCGGGAAATTTCACCGCCGTTTCCGGGTTTTCCCTTCTGGCGCGGGGGACTCTCTGGTCAGGGGTTGGAAACCCGGc                                                                                                >  2:1150538/1‑239 (MQ=255)
                                                                                                                                aTAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGACTCCAAAGGGTATCCACCATCGCATCCAGGTGGCCCGGGTTTTCCCCCCTGATGCGGGGAAAACCCGCGGCAGGCGGGAAAACTCCCCGTCGCTGCCGGTTTTTTCCCGCTGGCGCGGGGGAACTCCCggcccggcggggccaaaccggg                                                                                                 >  2:1026075/1‑219 (MQ=255)
                                                                                                                                    tGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGAGCCCAGGGGGCCCGGTTTTTCCCCCCTGATGCGGGGAAAACCAACGTCAGCGCGGAAATCTCACCGCCGCTGCCGGGTTTTCCCTCCTGGCGCGGGGAACCCCCCGTTGAGgtgtggggaaaccgcgtcac                                                                                            >  2:929712/1‑223 (MQ=255)
                                                                                                                                       tACCTCTGGTGAAGGATTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCggg                                                                                         <  1:1035792/239‑1 (MQ=255)
                                                                                                                                        tcgtgtGGGGAAGGAGGTGGTGGAGGGGTGTTGATGGGTTTGAAAATGGGAGGTGGGCGTTCCGCCGTGGAGGCGGGTGAACTGCAAGTGCTATCCAGCATCGCCTACAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTGCCggg                                                                                         <  2:1318307/233‑1 (MQ=255)
                                                                                                                                              ggTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAAATCCAAGTGATAACCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGGGTGAAATCTCCCCGTCGTTGCCGGGTTATCCCTGCTGGCGCGGGGAACTCTCCGTTCAGGCGTTGCAAACCTGGCTTCCGGGCggggta                                                                                  >  1:1155851/1‑235 (MQ=255)
                                                                                                                                                tGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTAt                                                                                 <  1:1040314/238‑1 (MQ=255)
                                                                                                                                                tGAAGGAGTTGGCGAAGGCGTCTTGAGGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTAt                                                                                 <  1:590509/238‑1 (MQ=255)
                                                                                                                                                 gAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATc                                                                                <  1:776099/238‑1 (MQ=255)
                                                                                                                                                          ggCGAGGGGGTTTTGATGGGTTGGAAAATGGGATCTGGGGTTTCTACCGGCGCGGCGGGGGAACTCCAAGGGTTATCCAGCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAAcg                                                                      <  2:135207/239‑1 (MQ=255)
                                                                                                                                                                    tCTTGGTGGGTTTGAAAATTGGAGCTGGGAGTTCGACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTACCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACt                                                             <  1:1211502/238‑1 (MQ=255)
                                                                                                                                                                    tCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGc                                                                                                                                                          <  1:680998/149‑1 (MQ=255)
                                                                                                                                                                     cTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGGACTCCAAGTGATATCCATCATCGCATCCCATGCGCCCCGTTTATCCCCGCTGGTGCGGGGAAAACCCGCGTCAGGCGTGAAATCTCACCGCCCGTGCCGGGTTATCCCCGCTGGGGCGGGGAGCTCTCGGGTCAGGGGTTGCAAACCCGGGCCCCCGGCGGGTTGTCCCCCCCCACGGGGGGGAcgca                                                           >  2:1122401/1‑236 (MQ=255)
                                                                                                                                                                       ttaggggTTTGAAAATGGGAGGGTGGAGTTCTCACGGAGAGGAGGGGGAACTCCAAGTGATATCCAGCATCGGATCCAGGGCGCCGGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGTGTGAAATCTCACCGTCGTTGACGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTa                                                         <  2:1077175/235‑1 (MQ=255)
                                                                                                                                                                       tGATGGGTTTGAAAAGGGGAGCTGGGAGTTCTGCCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTa                                                         <  1:1210686/239‑1 (MQ=255)
                                                                                                                                                                           gggTTTGAAAATGGGAGCTGGGATTTCCACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTcc                                                     <  1:500871/239‑1 (MQ=255)
                                                                                                                                                                           gggTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCCAGTGATATACCTCATCGCATCCCGTGCGCCCCGTTTATCCCCCCCGATGCGGGGGAAACCAGCGTCCTGCGGGGAATCTCACCCTCGTTTCCGGGTTATACCCGCTGTGTCGGGGGACCTCCGGGTCCGGCCGTTCAAAACGCGCTACCCGGCGGGTTTTCCCCCCCACCGCGGGGGAcacctgaccc                                                     >  2:1228291/1‑230 (MQ=255)
                                                                                                                                                                              tttGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAATGATATTCATCATCGCATCCCAAGCGCACGGTTTATACCCGCCGATGCGGGGAACACCAACGGCAAGCGGGAAATATCCCCGTCGGTGCCGGTTTATCCCGTCTGGCGCGGGGGACCCTCCGGTCCGGCGGTGCGAACCAGGCATCCGGGGGGGTACTCCCGCCCAACGAgggggaacaagagctacc                                                   >  2:294838/1‑224 (MQ=255)
                                                                                                                                                                              tttGAAAATGGGAACTGGGAATTCTACCGGAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCcat                                                   <  1:19627/238‑1 (MQ=255)
                                                                                                                                                                                 gaaagggggAGCTGGGAGTTCTACCGGAGGGGCGGGGGAACTCCAAGTGGTATCCATCATCGCAGCCAGGGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCAtt                                               <  2:265501/233‑1 (MQ=255)
                                                                                                                                                                                 gAAAATGGGAGCTGGGAGGTCTACCGCAGAGGCGGGGGAACTCCAAGTTATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCAtt                                               <  1:26099/239‑1 (MQ=255)
                                                                                                                                                                                               ggAGTTCTGCCGCAGGGGCGGGGGAACTCCAAGTGATATCCATCCTCGCAACCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCATCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTg                                  <  1:152419/238‑1 (MQ=255)
                                                                                                                                                                                                 aGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCcat                                                   <  1:482649/219‑1 (MQ=255)
                                                                                                                                                                                                    tCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCAGCGCATACAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACAACAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCGGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTcc                            <  1:1116574/239‑1 (MQ=255)
                                                                                                                                                                                                    tCTACCGCAGAGGCGGGGGAACTCCAAGTGATAGCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTcc                            <  1:103059/239‑1 (MQ=255)
                                                                                                                                                                                                     ctcccGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTcc                            <  1:1122401/235‑1 (MQ=255)
                                                                                                                                                                                                               ggCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGGGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTAt                    <  1:143536/236‑1 (MQ=255)
                                                                                                                                                                                                                     ggAACTCCAAGTGATATCCATCATCGCCTCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTg            <  1:1151015/238‑1 (MQ=255)
                                                                                                                                                                                                                        tcTCCAAGTGTTAACCATCATCGCATCCGGGGCGCCCGGTTTACCCCCGGTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGGCGGTGCCGGTTTACCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGgcg         <  2:1155851/237‑1 (MQ=255)
                                                                                                                                                                                                                        acttctaGTGGTTTCCATCCTCGCTTCCAGTGCGCCCGGTTTATCCCCGCAGATGCGGGGAACACCAGCGCCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGgcg         <  2:982739/232‑1 (MQ=255)
                                                                                                                                                                                                                        aCTCCAAGTGATATCCATCATCGCATCCAGGGCGCCCGGTTTATCCCCGGTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGgcg         <  1:577975/238‑1 (MQ=255)
                                                                                                                                                                                                                          tCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCATGCGTTGCAAACCGGCCTACCGGCGGGCTTATCCCCGCTGACCCGGGGAACTCGGAGCCCTTAATTCCACCCTTGTGTGCACCCCCGGTGTATccccccggggcggg      >  2:757426/1‑229 (MQ=255)
                                                                                                                                                                                                                           caaaGTGATAGCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCggg     <  1:39529/237‑1 (MQ=255)
                                                                                                                                                                                                                              aGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTCTCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCAACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGa   >  1:775961/1‑238 (MQ=255)
                                                                                                                                                                                                                                 gATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGaa  <  1:578867/236‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
AGCAACGCGCAAAAACACCTGCCTTCCCGGCGGGAAATAGCTGGCATGACGTAAGACTGGATAATCACCAACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAA  >  NC_000913/2877447‑2877907

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