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

ACCGGCTGCTATTACAGGAGAAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCT  >  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)
                  gaaacccTTTTCTCTGTTTATCTTTTTTTGCCATTTGCCTTTATGGTTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGtt                                                                                                                                                                                                                                        <  2:250599/231‑1 (MQ=255)
                  agaAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCTGCGGGTGTAAACCCACCGGCGGAGATGATTAAAAAGGtttt                                                                                                                                                                                                                                      >  1:139897/1‑237 (MQ=255)
                   gaAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCCCTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCAGCGGCGGAGATGATTTAAAGGGtttt                                                                                                                                                                                                                                      >  1:1083201/1‑236 (MQ=255)
                     aaCCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGCGGCGGGGGTGGTTAAAAGGGGTTCGCTg                                                                                                                                                                                                                                  >  1:175878/1‑239 (MQ=255)
                           tGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCCGCGGCGGAGATGATTAAAAGGGTTTCCCTGATCAtt                                                                                                                                                                                                                            >  1:1038391/1‑239 (MQ=255)
                                  gTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGACAACCCAGCGGCGGAGATGGTTTAAAGAGGTTTGCTGCTCActgctgcgt                                                                                                                                                                                                                    >  1:548389/1‑236 (MQ=255)
                                  gTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGCGGCGGGGATGATTAAAAGTGTTTCGCTGATCATTgtgtgtt                                                                                                                                                                                                                     >  1:956169/1‑233 (MQ=255)
                                    ttttttttttTGCCTTTTTTCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTtg                                                                                                                                                                                                                  <  2:927414/235‑1 (MQ=255)
                                    ttctttttATGGACATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTTTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTtg                                                                                                                                                                                                                  <  1:228827/235‑1 (MQ=255)
                                     tccTTTGAGTGCCATTTTCCTTTATGGTGATATTGAGATTTCAATCTTTAGTCAGGTTGCCCATGGATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTtgc                                                                                                                                                                                                                 <  1:822031/238‑1 (MQ=255)
                                       cTTTTTTTTCCTTTTTTCTTTATGTTTATATTGAGATTTCAATTTTTATTCAGGTTGCCCCTGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTtgctg                                                                                                                                                                                                             <  2:899163/239‑1 (MQ=255)
                                            aTTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCCTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCCGCGGCGGGGATGATTAAAAAAGTTTCTCTGGTCATTTCTGGGTtgtggggttt                                                                                                                                                                                                          >  2:55601/1‑231 (MQ=255)
                                                          ttATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCACCGGCGGGGGTGATTTAAAGTGGTTCGCTGTTAATTGCTGGTTTGCTGGGTTTTTTACCGGGCtttt                                                                                                                                                                                        >  1:690619/1‑239 (MQ=255)
                                                          ttATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAACGGCGGCGTTGATTAAACATGTTTTACTGATCATTGGTTGGTTTCTACTTTTTTTTCCGCGCtttt                                                                                                                                                                                        >  1:946790/1‑239 (MQ=255)
                                                             tGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGGAAAACCAGCGGCGGAGGTTGTTAAAAAAGTTTTCGCGAGCATTTTCGGTTTTGTGCTTTTTTTAACGAGGgttttt                                                                                                                                                                                      >  2:653421/1‑238 (MQ=255)
                                                                 tataTTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGGAAAACCAGCGGCGGGGGGTGTTAAAAATGTTTCGCTGGTAATTTCTGGGTTTGGTCCTTTTTTACAGGGAttttttttacc                                                                                                                                                                                 >  2:836124/1‑234 (MQ=255)
                                                                  atatTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCGAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGAATGGTTAAAAGTGTTTCGCGCTTAATTGGCGGTTTTCGTCCTTTTTTTCCCGGGTttttttccg                                                                                                                                                                                   >  1:82287/1‑234 (MQ=255)
                                                                       ttgttttatattttttttCAGTTTGCCCATTTATTGGGGGTTTTTCTGACCCTCGAGCTGGTTTTATTCACGTCCGTTATCGGTATGTCACTGGTAAGTAACCAGGGCTTTAAGAATTTCGGGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGGGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCt                                                                                                                                                                          <  2:506195/227‑1 (MQ=255)
                                                                       gcgaTTTCAATTTTTTTTCAGGTTGGCCATGTATTGGGGGTTTTGCTGACCCCCGGTGAGGTTATATTCACGGCAGTTAGCGGTATGTCACGGGTACGTAACCGTGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCTGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTTCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCt                                                                                                                                                                          <  2:817446/237‑1 (MQ=255)
                                                                        agaTTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTCAAAGGGTTTCGCTGGTCCTTGCTGGGTTGCTGCTTTTATTTCCGGGCTTTTTTTACGACTTCCTc                                                                                                                                                                         >  1:121873/1‑239 (MQ=255)
                                                                        agaTTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGCCGGCGGGGGAAAACCCAGCGGCGGGGGTGGATTAAAGGGTTTGCCTGGGCAATTGTGGGTTGCTCGTTTTATTTACGAGCttttttttcggcattctc                                                                                                                                                                          >  2:358057/1‑227 (MQ=255)
                                                                          gtttAAATTTTTAATCAGGTTGGTCATGTAGTGGGGGTATTGGTGACCCTCGTGCAGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATGTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCg                                                                                                                                                                        <  1:407709/237‑1 (MQ=255)
                                                                          aTTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTTAAAGAGTTTTCCTGGTCATTTCTGGGTTGCTGCTTTTATTAACTGGGTTTTTTACCGGCTTCCTCgg                                                                                                                                                                       >  1:573748/1‑239 (MQ=255)
                                                                           tttCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGGGGGTTTAAAAAAGGTTTTCTGGACATTGTTGGGTTTGTTCTTTTTTTACCGGGGttttttttcaccttcccccgg                                                                                                                                                                       >  2:631343/1‑224 (MQ=255)
                                                                                           ggTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCAGCGGCGGAGAAGGTTTAAAAGGTTTTGCTGGATGTTTCTTGTTTTTTGCTTTTTTTTCCGGGGGTTTTTTCCCGCCTCTCCCGGTTTTCtttttttttt                                                                                                                                                        >  2:356976/1‑238 (MQ=255)
                                                                                              tGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTgccgc                                                                                                                                                   <  1:712138/239‑1 (MQ=255)
                                                                                               gCCCATTGATTGGGGGTTTTGCTGACCCTCGGGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTgccgc                                                                                                                                                   <  1:584382/238‑1 (MQ=255)
                                                                                                     gTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCAACGGCGGAGATGGTTAAAAGGGGTTCGCTGGTCATTGCTGGTTTGCTGGTTTTATTACCGGGCTTTTTTACCGAATTCCTCGGGCTTTTAATTTTATTTCCGCCCGTGCa                                                                                                                                            >  1:950581/1‑239 (MQ=255)
                                                                                                      tATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTtgct                                                                                                                                                                                                               <  1:319303/175‑1 (MQ=255)
                                                                                                                aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTGAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGGTTAAAAGAGGTTCGCTGATCATTGCTGGTTTGCTGCCTTTTTTACCGGGCTTTTTTTCCCACTTACTCCGGCTTTCACTTTTTTTGCCCCCGGGGGaaaaaaaaacaa                                                                                                                                 >  2:415928/1‑231 (MQ=255)
                                                                                                                aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGGTTAAAAATGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTTCCGCCGGTGCAAAAACa                                                                                                                                      >  1:857576/1‑234 (MQ=255)
                                                                                                                aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGGAAACCCAGCGGCGGGGGTGATTTAAAAGGTTTCGCTGGTCTATGCTTGGTTGCTGCCTTTTTTACCCGGGGTTTTTTTCGCCTTTCTCGGGTTTTTTCTTTTTTTTTCCCCCGGGCaaaaaaaactta                                                                                                                                 >  2:374491/1‑231 (MQ=255)
                                                                                                                aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAACGGCGGAGATGATTAAAAATGGTTTCCCTATTACTGCCGGTTTGTTGCTTTTTTTACCGGGCGTTTTTTCCCGCTTTCCCCGTTTTTCTCTTTTTTTTCCCTCCGTGGAAAAACATCTgc                                                                                                                                 >  2:585509/1‑238 (MQ=255)
                                                                                                                 tttcTGTCCCTCGTGCTGGTTTTATTCACGTGAGTTTCCGGTATGTCACTGGTACGGAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGAc                                                                                                                               <  2:225421/236‑1 (MQ=255)
                                                                                                                  tgtttcccctgtggCTGTTTTTATTCCCGTTAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACa                                                                                                                              <  2:573748/226‑1 (MQ=255)
                                                                                                                  tGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAAACCAACGGGGGGAGTGATTAAAAAGGGTTTGCTGGTCATTTGTGGTTTTCTCGTTTTTTTTCCCGGGTTTTTTTACCGCCTCCCCCGGCTTTTTCTTTTTTTGCCCCCCGGGAaaaaaaacccgcaa                                                                                                                               >  2:876656/1‑229 (MQ=255)
                                                                                                                         ccTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGtcttc                                                                                                                                                                  <  1:571302/197‑1 (MQ=255)
                                                                                                                                   gTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAATGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACCTCCTCGGGCTTCTACTTTTATTTCCGCCGGTGCAAAAACATCTGGCAGTGGAGGTGATTCCg                                                                                                              >  1:89972/1‑237 (MQ=255)
                                                                                                                                       ttttcccggttgtttttggTTTGTCACTGGTACTTAACCAGGGGTTTAAGAATTTCGGGCTGATGCAGTAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTg                                                                                                        <  2:496011/222‑1 (MQ=255)
                                                                                                                                         ttCCCGTCATTTATCGGTATGTCACTGGTACGGAACCTGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAAGGGCGGCGGGTGAAAACCCAGGGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTgcg                                                                                                      <  2:1005252/239‑1 (MQ=255)
                                                                                                                                          gacaggcAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATGTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTgcg                                                                                                      <  1:243733/232‑1 (MQ=255)
                                                                                                                                           cACTTCAGTTATCGGTATGTCACTGGTACGTAACCGGGGCTTTAAGAATTTCGTGCTGATGCAGCCAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTGCCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTgcg                                                                                                      <  1:158367/237‑1 (MQ=255)
                                                                                                                                             cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGGTTAAAAGGGTTTCGCTGATCATTTCTGGGTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTTTACTTTTTTTTCCGCCGGTGCAAAAGACTCTTACCGTGAAGttttgtcccccttttggcttt                                                                                                  >  1:115720/1‑220 (MQ=255)
                                                                                                                                             cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCt                                                                                                                                                                          <  2:834789/169‑1 (MQ=255)
                                                                                                                                             cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCt                                                                                                                                                                          >  1:834789/1‑169 (MQ=255)
                                                                                                                                             cGTCAGTTAGCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCttt                                                                                                  <  1:374491/239‑1 (MQ=255)
                                                                                                                                              gTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCCCATTTTCGCtttt                                                                                                 >  1:397494/1‑239 (MQ=255)
                                                                                                                                              gTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGAGGCGGAGATGATTAAAAAGGTTTGGCCGGTCATTGCTGGTTTGCTGCCTTTATTACCCGGCTTTTTTACCTGCCTCCCCCGGTTTCTACATTTTTTTACGCCCGTGCAAAAACCTCTGTACTGGAAGTTGATGCCCGATTTTCGCtttt                                                                                                 >  2:124030/1‑239 (MQ=255)
                                                                                                                                                 aGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACCGTGAAGTTGATGCCGCCTTTTCGCtttttt                                                                                               >  1:360739/1‑236 (MQ=255)
                                                                                                                                                               cACTGTTACGTAACCGGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTggcgg                                                                                 <  1:653421/238‑1 (MQ=255)
                                                                                                                                                               cACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTggcggc                                                                                <  1:867587/239‑1 (MQ=255)
                                                                                                                                                                aCTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCCTTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCCACTTCCTCCGGCTTCTACTTTTTTTTCCGCCGGGGCAAAAACATCTGACCGTGCAAGTTATCCCCCATTTTTGCCTTTCTCTCCAGCCTGggggcg                                                                              >  2:867587/1‑239 (MQ=255)
                                                                                                                                                                    gTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTTCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTTATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGttt                                                                          >  1:875630/1‑239 (MQ=255)
                                                                                                                                                                          aaCCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCg                                                                                                                                                                                                                                     >  1:916290/1‑87 (MQ=255)
                                                                                                                                                                          aaCCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCg                                                                                                                                                                                                                                     <  2:916290/87‑1 (MQ=255)
                                                                                                                                                                            ccAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGTCAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCCTGCCTGGCGGCGGGTTTAGCGCCg                                                                  >  1:460613/1‑239 (MQ=255)
                                                                                                                                                                              aGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTTCTTTTATTTCCGCCGGTGCAAAAACATCTGACAGGGAAGTTTATTCCCCATTTGCGCTTTTTTTCCCTGCCTGGCCGCGGTTTTTGTGCCggg                                                                >  1:986935/1‑239 (MQ=255)
                                                                                                                                                                                ggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGGTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGGGCCAAAGCCTCTGACAGTGAAGTTGATGCCGCATTTTCGCTTTTCTCCCATGCCCGGCGGCGGTTTTAGAGCCcggg                                                               >  1:214990/1‑237 (MQ=255)
                                                                                                                                                                                    ttAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGGCGCCGGTGCCAAAGCATCTGACAGGGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGGTTTTAGGCCCgggccggg                                                           >  1:961343/1‑234 (MQ=255)
                                                                                                                                                                                    ttAAGAATTTCGGTCTTATGCAGCAAAAAATGGCGGCGGGTTAAACCCCCGCGGGGGAGATGATTAAAAGTTTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGt                                                          <  2:857576/239‑1 (MQ=255)
                                                                                                                                                                                    ttAAAAATTTCGGGGTGGTGGAGAAAAAAATGGGGGGGGGTGAAAACCCAGGGGCGGAGATGGTTAAAAGTGTTTGGCTGATCATTGCTGGTTTGCTTCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGt                                                          <  2:697582/239‑1 (MQ=255)
                                                                                                                                                                                      aaGAATTTCGTGCTGATGCAACAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTgg                                                        <  1:884267/239‑1 (MQ=255)
                                                                                                                                                                                           tttCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGGGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTTCCGACTTCCCCGGGCTTCTACTTTTATTGCCGCCGGGGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTTCGCTTTTCTCGCATGCCTGGCGGCGGGTTTATGGCCGGGGCCggggggggga                                                   >  1:12919/1‑237 (MQ=255)
                                                                                                                                                                                               gTGCTGATGCAACAAAAAATGGCGGCGGGTGAAAACCCAGGGGGGGGGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATAc                                               <  2:1038391/239‑1 (MQ=255)
                                                                                                                                                                                                    tgtggagcaAAAAATGGGGGGGGGTGAAAACCCAGGGGCGGAAATGATTAAAAGTTTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTg                                          <  2:950581/234‑1 (MQ=255)
                                                                                                                                                                                                    ggggcagcaAAAAGTGGCGGGGGTTGAAAACCCGGCGGCGGAGATGATTAAAAGTGTTTCGTTGTTATTTGCTGGTTTGCTTCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTg                                          <  2:946790/236‑1 (MQ=255)
                                                                                                                                                                                                     gaggagaaaaaaaTGGGGGGGGGTGAAAACCCGGCGGCGGAGATTATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGTTGGCGGTAATACTTTTg                                          <  2:175878/231‑1 (MQ=255)
                                                                                                                                                                                                     aTGCCGGAAAAAATGGCGGCGGGTGAAATCCCGGCGGGGGGGATGATTAAAAGTGTTTCGCTTATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTCCTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGa                                         <  2:177712/239‑1 (MQ=255)
                                                                                                                                                                                                      agcagcaAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTCCCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGa                                         <  1:585509/237‑1 (MQ=255)
                                                                                                                                                                                                            aaaaaaTGGGGGGGGGTGAAAACCCAGCGGGGGAGATGATTAAAAGTGTTTCGCTTATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAg                                  <  2:89972/239‑1 (MQ=255)
                                                                                                                                                                                                            aaaaaaTGGCGGCGGGTGAAAACCCAGTGGCGGATATGATTAAAAGTTTTTCGCTGATCATTTCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGa                                   <  1:358057/238‑1 (MQ=255)
                                                                                                                                                                                                               aaaTGGCGGCGGGTGAAAACCCAGCGGGGGATATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTAc                               <  1:474854/239‑1 (MQ=255)
                                                                                                                                                                                                                    aggaagggTGAAAACCCAGCGGCGGAAATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                          <  1:58712/234‑1 (MQ=255)
                                                                                                                                                                                                                     gggcggGTGAAAACCCAGGGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCGCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                          <  1:570542/237‑1 (MQ=255)
                                                                                                                                                                                                                     cggcggGTGAAACCCCAGCGGGGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGTTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                          <  1:481358/238‑1 (MQ=255)
                                                                                                                                                                                                                     cggcggGTGAAAACCCGGCGGGGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTCCCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                          <  2:139897/238‑1 (MQ=255)
                                                                                                                                                                                                                     cggcggGTGAAAACCCAGGGGCGGTGATGATTAAAAGTTTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                          <  1:356976/238‑1 (MQ=255)
                                                                                                                                                                                                                       gcggGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGaa                        <  1:55601/238‑1 (MQ=255)
                                                                                                                                                                                                                                      ggtgcggAAATGATTAAAAGTGTTTCGCTGATCATTTCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGATTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGa        <  1:736605/236‑1 (MQ=255)
                                                                                                                                                                                                                                      gcggggGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGa        <  1:631343/239‑1 (MQ=255)
                                                                                                                                                                                                                                            gagaTGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGGGGGCGGAAAACTTTTTATGGTGAGTACCAGCGAAAGGATGGTGAggcgggccgcct  >  1:647410/1‑228 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                        
ACCGGCTGCTATTACAGGAGAAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCT  >  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: