Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,993,145 C→A 31.2% R665S (CGT→AGT)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,993,1450CA31.2% 139.9 / 77.1 89R665S (CGT→AGT) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (35/25);  new base A (12/15);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

ATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTA  >  NC_000913/3992873‑3993418
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                 
atgtgctataGCTGGCGGGCTTCAGTTTTNNNGGGAATCAAAATTGCCTGTTAGGTAGCGTTGACCTTCTGTAACGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCAAGCAAGCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGt                                                                                                                                                                                                                                                                                 <  3:352456/266‑1 (MQ=255)
 tttCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGCTGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGt                                                                                                                                                                                                                                                                                 >  4:404893/1‑274 (MQ=255)
                          tttgggcaaaaTCAAAATTGCCGGGTCGGTGCGGTTGCCCTGCTGTACCCCAACTCGTGGAACGAAGTGCGTACGCTGCCCTTCAACGTCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCgg                                                                                                                                                                                                                                                       <  4:839282/266‑1 (MQ=255)
                                        aaaTTGCCTGGTAGGTAGCGTTGACCGGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGgcg                                                                                                                                                                                                                                          <  1:278765/274‑1 (MQ=255)
                                                   taggtagCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCa                                                                                                                                                                                                                                                                       >  4:276092/1‑234 (MQ=255)
                                                       tagCGTTGACCTGCTGTACCGCAACTCGTGGAACTAAGTGCGTACGCGGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATAAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATATGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTc                                                                                                                                                                                                                          <  2:563557/275‑1 (MQ=255)
                                                       tagCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCg                                                                                                                                                                                                                                                                   <  2:276095/234‑1 (MQ=255)
                                                        agCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCa                                                                                                                                                                                                                         <  2:785608/275‑1 (MQ=255)
                                                             tGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCt                                                                                                                                                                                                                    >  1:142776/1‑275 (MQ=255)
                                                                 ctgctgTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTTGGTCTTTCCAGCCCCCGCCAGGAAACCGGGGGTTCCAAGGCGCTGGGCGTTTCTGGTCAAACCTggg                                                                                                                                                                                                                 >  4:533375/1‑274 (MQ=255)
                                                                              aCTGGTGGAACGAAGGGGGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGa                                                                                                                                                                                                    <  4:759235/274‑1 (MQ=255)
                                                                                 cGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTggg                                                                                                                                                                                                                 >  4:581182/1‑258 (MQ=255)
                                                                                     gAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgtt                                                                                                                                                                                                             <  2:581185/258‑1 (MQ=255)
                                                                                             tGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCACGGAAACCGGGCGTTTCAAGGAGCTGCGCGTTTCTGGTCAACCCTGGGGGTTGTTCTTCGAACGCCTGAATTTATCg                                                                                                                                                                                    >  4:657548/1‑275 (MQ=255)
                                                                                                       gCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCCCCAGATAGCGTGGCAGTCTTCTGTTATCGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTGTCTGAGTGTATTGAATTGCGTCTTTCCTGCACCCGCCAGGGAACCGGGCGTTTCAAGGCGCGTCGCGTTTCCGTTCCAACCTGGGGGTTGTTCTTGCACCTCCTGATAGTCTCGGTACAGAAAc                                                                                                                                                                          >  4:62334/1‑275 (MQ=255)
                                                                                                                cGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACg                                                                                                                                                                 >  2:290337/1‑275 (MQ=255)
                                                                                                                  gCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCaa                                                                                                                                                                                                                        >  4:189962/1‑218 (MQ=255)
                                                                                                                     aGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGACGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTACTTAGTACTCGCGTGCAGCAGCTGGTTTNTGAGTCTATTGAATTGCGTCTGTCCAGCACCCGTCAGGATCCCGGGCGTTTCACGGCTCCGCGCGTTTGCGTCCNAACTTGTGGGTTGGTCGTCGCCGGCCTGCANNTATCGCTCCCGACATTGGAAAACGCCa                                                                                                                                                              >  3:153347/1‑273 (MQ=255)
                                                                                                                     aGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCt                                                                                                                                                                                                                              >  2:138911/1‑209 (MQ=255)
                                                                                                                      gCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCt                                                                                                                                                                                                                    <  2:189965/218‑1 (MQ=255)
                                                                                                                      gCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATc                                                                                                                                                            <  2:392552/274‑1 (MQ=255)
                                                                                                                         aTCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTc                                                                                                                                                                                                                          <  4:138908/209‑1 (MQ=255)
                                                                                                                              tgatCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGt                                                                                                                                                                                                                                                                                 <  6:102793/149‑1 (MQ=255)
                                                                                                                                    aaGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGttctt                                                                                                                                                                                                       >  4:437219/1‑217 (MQ=255)
                                                                                                                                     aGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCa                                                                                                                                                                                                                                                                          <  7:51364/149‑1 (MQ=255)
                                                                                                                                       cccTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTCTCTGGTCAANCCTGGGGTTTGTTCTTCGAACGCCTGCAGGTATCTGTACAGCAAGTGGAAAACGCAATCGGGTTTTCTGgcggg                                                                                                                                             >  3:292346/1‑270 (MQ=255)
                                                                                                                                        ccTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGaa                                                                                                                                                                                                   <  2:437222/217‑1 (MQ=255)
                                                                                                                                              aaCTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGcc                                                                                                                                                                                                                                                                 >  8:89719/1‑149 (MQ=255)
                                                                                                                                               aCTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCa                                                                                                                                                                                                                                                                <  8:28449/149‑1 (MQ=255)
                                                                                                                                                   ttCTCGGCAAAATGCATCAGGAGTCCGACCCGCCGGTTGGCGTGGAAGTTTTATGTTATAGCCAGAATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACaa                                                                                                                               <  4:590384/274‑1 (MQ=255)
                                                                                                                                                        ggCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGttctt                                                                                                                                                                                                       >  2:661007/1‑197 (MQ=255)
                                                                                                                                                            aaaTGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGaa                                                                                                                                                                                                   <  4:661004/197‑1 (MQ=255)
                                                                                                                                                             aaTGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGcc                                                                                                                    >  4:334174/1‑275 (MQ=255)
                                                                                                                                                                    cAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTg                                                                                                            <  2:404896/275‑1 (MQ=255)
                                                                                                                                                                     aGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTg                                                                                                            <  2:533378/275‑1 (MQ=255)
                                                                                                                                                                         cgccgcACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCggg                                                                                                                                                                                                                                                      >  2:924221/1‑133 (MQ=255)
                                                                                                                                                                         cgccgcACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAg                                                                                                         >  2:725401/1‑274 (MQ=255)
                                                                                                                                                                            cgcACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGttctt                                                                                                                                                                                                       >  4:401089/1‑177 (MQ=255)
                                                                                                                                                                            cgcACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCg                                                                                                                                       >  4:760097/1‑241 (MQ=255)
                                                                                                                                                                             gcACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGtt                                                                                                                                                                                                                                                  <  4:924218/133‑1 (MQ=255)
                                                                                                                                                                               aCCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGaa                                                                                                   <  2:334177/274‑1 (MQ=255)
                                                                                                                                                                                ccGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGaa                                                                                                                                                                                                   <  2:401092/177‑1 (MQ=255)
                                                                                                                                                                                ccGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATa                                                                                                                                   <  2:760100/241‑1 (MQ=255)
                                                                                                                                                                                           cGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAAAAAAc                                                                                                                             >  2:137490/1‑236 (MQ=255)
                                                                                                                                                                                             tGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCaat                                                                                    >  4:111891/1‑274 (MQ=255)
                                                                                                                                                                                             tGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCCAACGCCTGAATGTATCGGTACAGAAACTGGAAACAGCCATCGAGTTTTATGGCGAGATTTCGCATACAAAACTGCACGGCCTGTCAGTGCAGCGTGAAAcaatactcggcaaa                                                                                    >  4:11912/1‑262 (MQ=255)
                                                                                                                                                                                             tGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCGTCGAACGCCTGGATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGCTTTCGCATAACAAACTGCACGGCCTTGCATTGCAGGTTGAAAGCAATCACGtacaa                                                                                    >  1:243550/1‑271 (MQ=255)
                                                                                                                                                                                              ggAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTggg                                                                                                                                                                                                                 <  5:14685/149‑1 (MQ=255)
                                                                                                                                                                                               gAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTgggg                                                                                                                                                                                                                >  5:55141/1‑149 (MQ=255)
                                                                                                                                                                                                   tcttctGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgttg                                                                                                                                                                                                            >  5:86708/1‑149 (MQ=255)
                                                                                                                                                                                                     ttctGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTg                                                                                                                                                                                                                             >  1:266823/1‑130 (MQ=255)
                                                                                                                                                                                                     ttctGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGGTCTTCGAACGCCCGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTGCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACccgc                                                                            >  1:110543/1‑275 (MQ=255)
                                                                                                                                                                                                      tctGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAATTTACcgg                                                                             <  2:111894/273‑1 (MQ=255)
                                                                                                                                                                                                         gTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCa                                                                                                                                                                                                                         <  3:266823/130‑1 (MQ=255)
                                                                                                                                                                                                             taGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGttcttc                                                                                                                                                                                                      >  6:66741/1‑145 (MQ=255)
                                                                                                                                                                                                             taGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGttcttc                                                                                                                                                                                                      <  5:66741/145‑1 (MQ=255)
                                                                                                                                                                                                             taGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGaa                                                                                                                                                                     >  2:862427/1‑178 (MQ=255)
                                                                                                                                                                                                                 cAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACg                                                                                                                                                                 <  4:862424/178‑1 (MQ=255)
                                                                                                                                                                                                                  aGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTgg                                                                                                                                                                       >  2:473307/1‑171 (MQ=255)
                                                                                                                                                                                                                    cATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCAt                                                                                                                                    >  1:165573/1‑204 (MQ=255)
                                                                                                                                                                                                                    cATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCtt                                                             >  2:496914/1‑275 (MQ=255)
                                                                                                                                                                                                                    cATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGCCGCGATTTc                                                                                                                                        >  1:43236/1‑200 (MQ=255)
                                                                                                                                                                                                                      tCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGaaaa                                                                                                                                                                   <  4:473304/171‑1 (MQ=255)
                                                                                                                                                                                                                        tGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGaaa                                                                                                                                                                                                                                                           >  4:985986/1‑81 (MQ=255)
                                                                                                                                                                                                                          cgcgGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGttct                                                                                                                                                                                                        >  4:809072/1‑130 (MQ=255)
                                                                                                                                                                                                                            cgGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCgg                                                                                                                                                                                                                                                       <  2:985989/81‑1 (MQ=255)
                                                                                                                                                                                                                              gCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGa                                                                                                                                                                                                    <  2:809075/130‑1 (MQ=255)
                                                                                                                                                                                                                                        tACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAAc                                                                                                                                                                                                                                                          >  4:500123/1‑66 (MQ=255)
                                                                                                                                                                                                                                         aCTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCCGAATGTATCGGTACAGAAACTGGAACACGCCATCGAGTTTTCTGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGGGGTGGACGGCCTTGCCAGCGAAGGGCGCATCCa                                        >  4:780493/1‑275 (MQ=255)
                                                                                                                                                                                                                                         aCTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGc                                                                                                                          >  4:782009/1‑193 (MQ=255)
                                                                                                                                                                                                                                         aCTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGc                                                                                                                     >  2:789399/1‑198 (MQ=255)
                                                                                                                                                                                                                                            cgcgTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCggg                                                                                                                                                                                                                                                      <  2:500126/66‑1 (MQ=255)
                                                                                                                                                                                                                                             gcgTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACgg                                                                                                                      <  2:782012/193‑1 (MQ=255)
                                                                                                                                                                                                                                             gcgTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGt                                                                                                                 <  4:789396/198‑1 (MQ=255)
                                                                                                                                                                                                                                              cgTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGGGGTGGACGGCTTTGCCAGCGag                                                   >  4:89082/1‑258 (MQ=255)
                                                                                                                                                                                                                                                  cagcaACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGa                                               <  2:89085/259‑1 (MQ=255)
                                                                                                                                                                                                                                                      aaCTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAg                                                                                                                                                         >  8:93407/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      aaCTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAg                                                                                                                                                         >  5:64980/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           gTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAAc                                                                                                                             >  1:27708/1‑172 (MQ=255)
                                                                                                                                                                                                                                                             ttCTGAGTTGATTTAATTGAGTCTTTCCAGAACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACg                     <  1:292346/274‑1 (MQ=255)
                                                                                                                                                                                                                                                               cTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCa                                                                                                                         <  3:27708/172‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gtATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACgag            >  2:81581/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                     gtATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACgag            >  4:58130/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                      tATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCa                                                                                                                                                                                                                                               >  1:204727/1‑47 (MQ=255)
                                                                                                                                                                                                                                                                       aTTGAATTGGGTCTTTCCCGAACCCGCCAGGGACCCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGGTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGaa          <  4:578435/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ttGAATTGCGTGTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGTTGCAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGTTTTCg                                                                                                                                       >  6:20993/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          gAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAggt                                                                                                                                                                                                                                           <  3:204727/47‑2 (MQ=255)
                                                                                                                                                                                                                                                                            atggCGGCTTTCCAGCAACCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCt     <  4:81578/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                               gCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTa  <  2:58132/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                                cGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGttctt                                                                                                                                                                                                       >  2:463997/1‑77 (MQ=255)
                                                                                                                                                                                                                                                                                cGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGc                                                                                                                          >  4:441361/1‑154 (MQ=255)
                                                                                                                                                                                                                                                                                cGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGAGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCATGAATGTATCGGTACAGAAACTGGCAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACACACTGCACGGCCTGTCAGTGCAGGTGGAAACCAATCAAGTCAAATTACCGGCCGTGGTGGACGTCTTTGCCAGCGAAGGGAACGACCAGTTCTTTTTCGAAGACACGCAAGACGAGAATGGCTTTa  >  4:423060/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                 
ATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTA  >  NC_000913/3992873‑3993418

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: