Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,368,876 A→G 12.9% E71G (GAG→GGG)  fxsA → suppressor of F exclusion of phage T7

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,368,8760AG12.9% 153.2 / 2.5 81E71G (GAG→GGG) fxsAsuppressor of F exclusion of phage T7
Reads supporting (aligned to +/- strand):  ref base A (27/40);  new base G (10/2);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 9.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGCTGCTATTACAGGAGAAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCT  >  NC_000913/4368639‑4369113
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                        
cccgggttttttttcagggaaaaCCTTTGGGTTGGTTTCCTTTTATTTCCATTTTCCTTTATGGCTATATTGAGATTGCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGgag                                                                                                                                                                                                                                                         <  2:539129/219‑1 (MQ=255)
acggggtgctTTTGCAGGGAAAACCTTGGGTCTTGTTCACTTTTATTGCCATTTTCCTTTATGGCTAAAATGAGATTTCAAGCTTTATGCAGGTTGCCCGGGTATTGGGGGTATTGCTGACCCTCGGGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATGTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGgag                                                                                                                                                                                                                                                         <  2:71559/233‑1 (MQ=255)
aCCGGCTGTTTTTACAGGAAAAACCTTTGCGCTTGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGgag                                                                                                                                                                                                                                                         <  2:128114/239‑1 (MQ=255)
aCCGGCTGCTATTACAGGAGAAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGgag                                                                                                                                                                                                                                                         <  1:107841/239‑1 (MQ=255)
aCCGGCTGCTATTACAGGAGAAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGgag                                                                                                                                                                                                                                                         <  1:664877/239‑1 (MQ=255)
     ctgctATTACAGGAGAAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCCACGGCGGAgatgat                                                                                                                                                                                                                                                    >  1:927414/1‑239 (MQ=255)
                  gaaacccTTTTCTCTGTTTATCTTTTTTTGCCATTTGCCTTTATGGTTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑tt                                                                                                                                                                                                                                        <  2:250599/231‑1 (MQ=255)
                  agaAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCTGCGGGTGTAAACCCACCGGCGGAGATGATTAAAAAG‑G‑tttt                                                                                                                                                                                                                                      >  1:139897/1‑237 (MQ=255)
                   gaAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCCCTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCAGCGGCGGAGATGATTTAAAGG‑G‑tttt                                                                                                                                                                                                                                      >  1:1083201/1‑236 (MQ=255)
                     aaCCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGCGGCGGGGGTGGTTAAAAGG‑G‑GTTCGCTg                                                                                                                                                                                                                                  >  1:175878/1‑239 (MQ=255)
                           tGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCCGCGGCGGAGATGATTAAAAGG‑G‑TTTCCCTGATCAtt                                                                                                                                                                                                                            >  1:1038391/1‑239 (MQ=255)
                                  gTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGACAACCCAGCGGCGGAGATGGTTTAAAGA‑G‑GTTTGCTGCTCActgct‑gcgt                                                                                                                                                                                                                    >  1:548389/1‑236 (MQ=255)
                                  gTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGCGGCGGGGATGATTAAAAGT‑G‑TTTCGCTGATCATTgtgtgtt                                                                                                                                                                                                                     >  1:956169/1‑233 (MQ=255)
                                    ttttttttttTGCCTTTTTTCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTtg                                                                                                                                                                                                                  <  2:927414/235‑1 (MQ=255)
                                    ttctttttATGGACATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTTTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTtg                                                                                                                                                                                                                  <  1:228827/235‑1 (MQ=255)
                                     tccTTTGAGTGCCATTTTCCTTTATGGTGATATTGAGATTTCAATCTTTAGTCAGGTTGCCCATGGATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTtgc                                                                                                                                                                                                                 <  1:822031/238‑1 (MQ=255)
                                       cTTTTTTTTCCTTTTTTCTTTATGTTTATATTGAGATTTCAATTTTTATTCAGGTTGCCCCTGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTtgc‑t‑g                                                                                                                                                                                                             <  2:899163/239‑1 (MQ=255)
                                            aTTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCCTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCCGCGGCGGGGATGATTAAAAAA‑G‑TTTCTCTGGTCATTTCT‑GGGTtgtggggttt                                                                                                                                                                                                          >  2:55601/1‑231 (MQ=255)
                                                          ttATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCACCGGCGGGGGTGATTTAAAGT‑G‑GTTCGCTGTTAATTGCT‑GGTTTGC‑T‑GGGTTTTTTACC‑GGGC‑tttt                                                                                                                                                                                        >  1:690619/1‑239 (MQ=255)
                                                          ttATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAACGGCGGCGTTGATTAAACAT‑G‑TTTTACTGATCATTGGT‑TGGTTTC‑T‑ACTTTTTTTTCC‑GCGC‑tttt                                                                                                                                                                                        >  1:946790/1‑239 (MQ=255)
                                                             tGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGGAAAACCAGCGGCGGAGGTTGTTAAAAAA‑GTTTTCGC‑GAGCATTTTC‑GGTTTTG‑T‑GCTTTTTTTAAC‑GAGG‑gttttt                                                                                                                                                                                      >  2:653421/1‑238 (MQ=255)
                                                                 tataTTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGGAAAACCAGCGGCGGGGGGTGTTAAAAAT‑G‑TTTCGCTGGTAATTTCTGGGTTTGG‑T‑CCTTTT‑TTACA‑GGGA‑ttttttttacc                                                                                                                                                                                 >  2:836124/1‑234 (MQ=255)
                                                                  atatTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCGAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGAATGGTTAAAAGT‑G‑TTTCGCGCTTAATTGGC‑GGTTTTCGT‑CCTTTTTTTCCC‑GGGT‑ttttttccg                                                                                                                                                                                   >  1:82287/1‑234 (MQ=255)
                                                                       ttgttttatattttttttCAGTTTGCCCATTTATTGGGGGTTTTTCTGACCCTCGAGCTGGTTTTATTCACGTCCGTTATCGGTATGTCACTGGTAAGTAACCAGGGCTTTAAGAATTTCGGGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGGGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCt                                                                                                                                                                          <  2:506195/227‑1 (MQ=255)
                                                                       gcgaTTTCAATTTTTTTTCAGGTTGGCCATGTATTGGGGGTTTTGCTGACCCCCGGTGAGGTTATATTCACGGCAGTTAGCGGTATGTCACGGGTACGTAACCGTGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCTGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTTC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCt                                                                                                                                                                          <  2:817446/237‑1 (MQ=255)
                                                                        agaTTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTCAAAGG‑G‑TTTCGCTGGTCCTTGCT‑GGGTTGC‑T‑GCTTTTATTTCC‑GGGC‑TTTTTTTAC‑GACTTCCTc                                                                                                                                                                         >  1:121873/1‑239 (MQ=255)
                                                                        agaTTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGCCGGCGGGGGAAAACCCAGCGGCGGGGGTGGATTAAAGG‑G‑TTTGCCTGGGCAATTGT‑GGGTTGC‑T‑CGTTTTATTTAC‑GAGC‑ttttttttcggcattctc                                                                                                                                                                          >  2:358057/1‑227 (MQ=255)
                                                                          gtttAAATTTTTAATCAGGTTGGTCATGTAGTGGGGGTATTGGTGACCCTCGTGCAGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATGTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCg                                                                                                                                                                        <  1:407709/237‑1 (MQ=255)
                                                                          aTTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTTAAAGA‑G‑TTTTCCTGGTCATTTCT‑GGGTTGC‑T‑GCTTTTATTAAC‑TGGG‑TTTTTTACC‑GGCTTCCTCgg                                                                                                                                                                       >  1:573748/1‑239 (MQ=255)
                                                                           tttCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGGGGGTTTAAAAAA‑G‑GTTTTCTGGACATTGTT‑GGGTTTG‑T‑TCTTTTTTTACC‑GGGG‑ttttttttcaccttcccccgg                                                                                                                                                                       >  2:631343/1‑224 (MQ=255)
                                                                                           ggTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCAGCGGCGGAGAAGGTTTAAAAG‑G‑TTTTGCTGGATGTTTCT‑TGTTTTT‑T‑GCTTTTTTTTCCGGGGG‑TTTTTTCCC‑GCCTCTCCCGGTTTTCtttttttttt                                                                                                                                                        >  2:356976/1‑238 (MQ=255)
                                                                                              tGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTgccgc                                                                                                                                                   <  1:712138/239‑1 (MQ=255)
                                                                                               gCCCATTGATTGGGGGTTTTGCTGACCCTCGGGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTgccgc                                                                                                                                                   <  1:584382/238‑1 (MQ=255)
                                                                                                     gTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCAACGGCGGAGATGGTTAAAAGG‑G‑GTTCGCTGGTCATTGCT‑GGTTTGC‑T‑GGTTTTATTACC‑GGGC‑TTTTTTACC‑GAATTCCTCGGGCTTTTAATTTTATTTCCGCCCGTGCa                                                                                                                                            >  1:950581/1‑239 (MQ=255)
                                                                                                      tATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTtgc‑t                                                                                                                                                                                                               <  1:319303/175‑1 (MQ=255)
                                                                                                                aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTGAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGGTTAAAAGA‑G‑GTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCCTTTTTTACC‑GGGC‑TTTTTTTCC‑CACTTACTCCGGCTTTCACTTTTTTTGCCCCCGGGGGaaaaaaaaacaa                                                                                                                                 >  2:415928/1‑231 (MQ=255)
                                                                                                                aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGGTTAAAAAT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTTCCGCCGGTGCAAAAACa                                                                                                                                      >  1:857576/1‑234 (MQ=255)
                                                                                                                aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGGAAACCCAGCGGCGGGGGTGATTTAAAAG‑G‑TTTCGCTGGTCTATGCT‑TGGTTGC‑T‑GCCTTTTTTACC‑CGGG‑GTTTTTTTC‑GCCTTTCTCGGGTTTTTTCTTTTTTTTTCCCCCGGGCaaaaaaaactta                                                                                                                                 >  2:374491/1‑231 (MQ=255)
                                                                                                                aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAACGGCGGAGATGATTAAAAATGG‑TTTCCCT‑ATTACTGCC‑GGTTTGT‑T‑GCTTTTTTTACC‑GGGCGTTTTTTCCC‑GCTTTCCCCGTTTTTCT‑CTTTTTTTTCCCTCCGTGGAAAAACATCTgc                                                                                                                                 >  2:585509/1‑238 (MQ=255)
                                                                                                                 tttcTGTCCCTCGTGCTGGTTTTATTCACGTGAGTTTCCGGTATGTCACTGGTACGGAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑Ac                                                                                                                               <  2:225421/236‑1 (MQ=255)
                                                                                                                  tgtttcccctgtggCTGTTTTTATTCCCGTTAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACa                                                                                                                              <  2:573748/226‑1 (MQ=255)
                                                                                                                  tGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAAACCAACGGGGGGAGTGATTAAAAAG‑G‑GTTTGCTGGTCATTTGT‑GGTTTTC‑TCGTTTTTTTTCCC‑GGGT‑TTTTTTACC‑GCCTCCCCCGG‑CTTTTTCTTTTTTTGCCCCCCGGGAaaaaaaacccgcaa                                                                                                                               >  2:876656/1‑229 (MQ=255)
                                                                                                                         ccTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGtcttc                                                                                                                                                                  <  1:571302/197‑1 (MQ=255)
                                                                                                                                   gTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAAT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACCTCCTCGGGCTTCTACTTTTATTTCCGCCGGTGCAAAAACATCTG‑GCAGTG‑GAGGTGATTCCg                                                                                                              >  1:89972/1‑237 (MQ=255)
                                                                                                                                       ttttcccggttgtttttggTTTGTCACTGGTACTTAACCAGGGGTTTAAGAATTTCGGGCTGATGCAGTAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTg                                                                                                        <  2:496011/222‑1 (MQ=255)
                                                                                                                                         ttCCCGTCATTTATCGGTATGTCACTGGTACGGAACCTGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAAGGGCGGCGGGTGAAAACCCAGGGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTgcg                                                                                                      <  2:1005252/239‑1 (MQ=255)
                                                                                                                                          gacaggcAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATGTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTgcg                                                                                                      <  1:243733/232‑1 (MQ=255)
                                                                                                                                           cACTTCAGTTATCGGTATGTCACTGGTACGTAACCGGGGCTTTAAGAATTTCGTGCTGATGCAGCCAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTGCC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTgcg                                                                                                      <  1:158367/237‑1 (MQ=255)
                                                                                                                                             cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGGTTAAAAGG‑G‑TTTCGCTGATCATTTCT‑GGGTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTTTACTTTTTTTTCCGCCGGTGCAAAAGACTCTT‑ACCGTG‑AAGttttgtcccccttttggcttt                                                                                                  >  1:115720/1‑220 (MQ=255)
                                                                                                                                             cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCt                                                                                                                                                                          <  2:834789/169‑1 (MQ=255)
                                                                                                                                             cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCt                                                                                                                                                                          >  1:834789/1‑169 (MQ=255)
                                                                                                                                             cGTCAGTTAGCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCttt                                                                                                  <  1:374491/239‑1 (MQ=255)
                                                                                                                                              gTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCCCATTTTCGCtttt                                                                                                 >  1:397494/1‑239 (MQ=255)
                                                                                                                                              gTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGAGGCGGAGATGATTAAAAAG‑G‑TTTGGCCGGTCATTGCT‑GGTTTGC‑T‑GCCTTTATTACC‑CGGC‑TTTTTTACCTGCCTCCCCCGGT‑TTCTACATTTTTTTACGCCCGTGCAAAAACCTCTGTACTG‑G‑AAGTTGATGCCCGATTTTCGCtttt                                                                                                 >  2:124030/1‑239 (MQ=255)
                                                                                                                                                 aGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACCGTG‑AAGTTGATGCCGCCTTTTCGCtttttt                                                                                               >  1:360739/1‑236 (MQ=255)
                                                                                                                                                               cACTGTTACGTAACCGGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTggcgg                                                                                 <  1:653421/238‑1 (MQ=255)
                                                                                                                                                               cACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTggcggc                                                                                <  1:867587/239‑1 (MQ=255)
                                                                                                                                                                aCTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCCTTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑CACTTCCTCCGGCTTCTACTTTTTTTTCCGCCGGGGCAAAAACATCTG‑ACCGTGCAAGTT‑ATCCCCCATTTTTGCCTTTCTCTCCAGCCTGggggc‑g                                                                              >  2:867587/1‑239 (MQ=255)
                                                                                                                                                                    gTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTTCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTTATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGttt                                                                          >  1:875630/1‑239 (MQ=255)
                                                                                                                                                                          aaCCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCg                                                                                                                                                                                                                                     >  1:916290/1‑87 (MQ=255)
                                                                                                                                                                          aaCCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCg                                                                                                                                                                                                                                     <  2:916290/87‑1 (MQ=255)
                                                                                                                                                                            ccAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑TCAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCCTGCCTGGCGGC‑GGGTTTAGCGCCg                                                                  >  1:460613/1‑239 (MQ=255)
                                                                                                                                                                              aGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTTCTTTTATTTCCGCCGGTGCAAAAACATCTG‑ACAGGG‑AAGTTTATTCCCCATTTGCGCTTTTTTTCCCTGCCTGGCCGC‑GGTTTTTGTGCCggg                                                                >  1:986935/1‑239 (MQ=255)
                                                                                                                                                                                ggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGGTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGGGCCAAAGCCTCTG‑ACAGTG‑AAGTTGATGCCGCATTTTCGCTTTTCTCCCATGCCCGGCGGC‑GGTTTTAGAGCCcggg                                                               >  1:214990/1‑237 (MQ=255)
                                                                                                                                                                                    ttAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGGCGCCGGTGCCAAAGCATCTG‑ACAGGG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGGTTTTAGGCCCgggccggg                                                           >  1:961343/1‑234 (MQ=255)
                                                                                                                                                                                    ttAAGAATTTCGGTCTTATGCAGCAAAAAATGGCGGCGGGTTAAACCCCCGCGGGGGAGATGATTAAAAGT‑T‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGt                                                          <  2:857576/239‑1 (MQ=255)
                                                                                                                                                                                    ttAAAAATTTCGGGGTGGTGGAGAAAAAAATGGGGGGGGGTGAAAACCCAGGGGCGGAGATGGTTAAAAGT‑G‑TTTGGCTGATCATTGCT‑GGTTTGC‑T‑TCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGt                                                          <  2:697582/239‑1 (MQ=255)
                                                                                                                                                                                      aaGAATTTCGTGCTGATGCAACAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTgg                                                        <  1:884267/239‑1 (MQ=255)
                                                                                                                                                                                           tttCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGG‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTTCC‑GACTTCCCCGGGCTTCTACTTTTATTGCCGCCGGGGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTTCGCTTTTCTCGCATGCCTGGCGGC‑GGGTTTATGGCCGGGGCCggggggggga                                                   >  1:12919/1‑237 (MQ=255)
                                                                                                                                                                                               gTGCTGATGCAACAAAAAATGGCGGCGGGTGAAAACCCAGGGGGGGGGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATAc                                               <  2:1038391/239‑1 (MQ=255)
                                                                                                                                                                                                    tgtggagcaAAAAATGGGGGGGGGTGAAAACCCAGGGGCGGAAATGATTAAAAGT‑T‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTg                                          <  2:950581/234‑1 (MQ=255)
                                                                                                                                                                                                    ggggcagcaAAAAGTGGCGGGGGTTGAAAACCCGGCGGCGGAGATGATTAAAAGT‑G‑TTTCGTTGTTATTTGCT‑GGTTTGC‑T‑TCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTg                                          <  2:946790/236‑1 (MQ=255)
                                                                                                                                                                                                     gaggagaaaaaaaTGGGGGGGGGTGAAAACCCGGCGGCGGAGATTATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGTTGGCGGTAATACTTTTg                                          <  2:175878/231‑1 (MQ=255)
                                                                                                                                                                                                     aTGCCGGAAAAAATGGCGGCGGGTGAAATCCCGGCGGGGGGGATGATTAAAAGT‑G‑TTTCGCTTATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTCCTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGa                                         <  2:177712/239‑1 (MQ=255)
                                                                                                                                                                                                      agcagcaAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTCCC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGa                                         <  1:585509/237‑1 (MQ=255)
                                                                                                                                                                                                            aaaaaaTGGGGGGGGGTGAAAACCCAGCGGGGGAGATGATTAAAAGT‑G‑TTTCGCTTATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAg                                  <  2:89972/239‑1 (MQ=255)
                                                                                                                                                                                                            aaaaaaTGGCGGCGGGTGAAAACCCAGTGGCGGATATGATTAAAAGT‑T‑TTTCGCTGATCATTTCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGa                                   <  1:358057/238‑1 (MQ=255)
                                                                                                                                                                                                               aaaTGGCGGCGGGTGAAAACCCAGCGGGGGATATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTAc                               <  1:474854/239‑1 (MQ=255)
                                                                                                                                                                                                                    aggaagggTGAAAACCCAGCGGCGGAAATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                          <  1:58712/234‑1 (MQ=255)
                                                                                                                                                                                                                     gggcggGTGAAAACCCAGGGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCGCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                          <  1:570542/237‑1 (MQ=255)
                                                                                                                                                                                                                     cggcggGTGAAACCCCAGCGGGGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GTTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                          <  1:481358/238‑1 (MQ=255)
                                                                                                                                                                                                                     cggcggGTGAAAACCCGGCGGGGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTCCC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                          <  2:139897/238‑1 (MQ=255)
                                                                                                                                                                                                                     cggcggGTGAAAACCCAGGGGCGGTGATGATTAAAAGT‑T‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                          <  1:356976/238‑1 (MQ=255)
                                                                                                                                                                                                                       gcggGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGaa                        <  1:55601/238‑1 (MQ=255)
                                                                                                                                                                                                                                      ggtgcggAAATGATTAAAAGT‑G‑TTTCGCTGATCATTTCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GATTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGa        <  1:736605/236‑1 (MQ=255)
                                                                                                                                                                                                                                      gcggggGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGa        <  1:631343/239‑1 (MQ=255)
                                                                                                                                                                                                                                            gagaTGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGGGGGCGGAAAACTTTTTATGGTGAGTACCAGCGAAAGGATGGTGAggcgggccgcct  >  1:647410/1‑228 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                        
ACCGGCTGCTATTACAGGAGAAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGT‑G‑TTTCGCTGATCATTGCT‑GGTTTGC‑T‑GCTTTTATTACC‑GGGC‑TTTTTTACC‑GACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTG‑ACAGTG‑AAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCT  >  NC_000913/4368639‑4369113

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