Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 417,574 C→G 7.2% R145G (CGA→GGA)  phoB → response regulator in two‑component regulatory system with PhoR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913417,5740CG 7.2% 228.8 / 7.6 69R145G (CGA→GGA) phoBresponse regulator in two‑component regulatory system with PhoR
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base G (2/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

CCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGATATCGTTTTTCAACCCGCTTTTAACGCCTTGCTCATCGGACGCAGAGCAGGG  >  NC_000913/417275‑417859
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                             
ccTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCg                                                                                                                                                                                                                                                                                              <  2:237719/301‑1 (MQ=255)
   ggCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGCAGGAGCTGGTGGCGCGCATCAAAGCGGTAATGCGCCGTCTTTCGCCAATGGCGGTGGACGAGGTGATTGAGATGCAGGGATTAAGTCTCCACCCGACATCTCACCGCGt                                                                                                                                                                                                                                                                                           >  1:543050/1‑301 (MQ=255)
            tAATTTTCCTGGACTGGATGTTCGCTTGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATGCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATggcggg                                                                                                                                                                                                                                                                                  <  1:123164/301‑1 (MQ=255)
            tAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGCCATCACTCGGTGAGATgtcggg                                                                                                                                                                                                                                                                                  >  2:64053/1‑301 (MQ=255)
            tAATTCTCCGCGACGGGAGGTTACCTGGCGGCTCCGGGTTCCAGTGCATCAAACACCGCAAGCGCGAGGCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGCCATCACTCGGTGAGATgtcggg                                                                                                                                                                                                                                                                                  <  1:12083/301‑1 (MQ=255)
             ctcgacgctcttcggttctttcccTGGCGGCTCCGGTATCCAGTGCATCAAACACCGCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTGTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGCCATCACTCGGTGAGATgtcgggt                                                                                                                                                                                                                                                                                 <  1:589649/279‑2 (MQ=255)
                cctcctgGACGTGATGTTACCGGGCGGCTCCGGTATCCAGTTCATCAAACACCGCAAGCGCGTGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGAGCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGCCATCACTCGGTGAGATgtcgggtcga                                                                                                                                                                                                                                                                              <  1:63985/300‑5 (MQ=255)
                 cgccttgccTGGATGTGTCCTGGTGGTTCCTGTATCCAGTTCATCAAACTCCCCAAGCGCGAGTTGATTGCCCGGGATATTCCAGTGGTGTTTTGGACCGCCAGGGGGGAAGAAGAAGATCGCGTTCGGGGCCTTGCACCCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAag                                                                                                                                                                                                                                                                             <  1:535289/293‑1 (MQ=255)
                       gACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGc                                                                                                                                                                                                                                                                       <  2:543723/301‑1 (MQ=255)
                       gACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGAGAAGAGCCGc                                                                                                                                                                                                                                                                       >  2:148857/1‑301 (MQ=255)
                         cTGGTTGTTACCTGGCTGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGTCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTg                                                                                                                                                                                                                                                                     <  2:44180/301‑1 (MQ=255)
                          ggggggttttcTGGTGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTgg                                                                                                                                                                                                                                                                    <  1:172494/291‑1 (MQ=255)
                                ttACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCACTCGGTGAGATgtcgggtagatcggaagagcacacgg                                                                                                                                                                                                                                                              >  2:590410/1‑281 (MQ=255)
                                                    ccAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCCGGAGATGGGGCCGACAGAATTTAAActg                                                                                                                                                                                                                                          >  1:562943/1‑301 (MQ=255)
                                                        ttCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACtgctgc                                                                                                                                                                                                                                      <  2:563657/301‑1 (MQ=255)
                                                            tcaAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACtt                                                                                                                                                                                                                                  >  2:407023/1‑301 (MQ=255)
                                                               aaaacCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACttttt                                                                                                                                                                                                                               <  2:370931/298‑1 (MQ=255)
                                                                        aGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCATATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCa                                                                                                                                                                                                                      >  2:175164/1‑301 (MQ=255)
                                                                                       cccGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTa                                                                                                                                                                                                       >  2:413588/1‑301 (MQ=255)
                                                                                             cgtgtgctcttccgATCTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCcg                                                                                                                                                                                                 <  2:596159/288‑1 (MQ=255)
                                                                                                 tCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCgagaa                                                                                                                                                                                             >  2:554674/1‑300 (MQ=255)
                                                                                                   cAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATGGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCGGAGCGCGTGTACAGCCGCGagca                                                                                                                                                                                           <  2:455048/301‑1 (MQ=255)
                                                                                                         ggATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCATATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGtt                                                                                                                                                                                     <  1:174967/300‑1 (MQ=255)
                                                                                                          gATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTa                                                                                                                                                                                    >  2:342547/1‑301 (MQ=255)
                                                                                                            tGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAgatcggaagagagtcgggtagggaaagagg                                                                                                                                                                                  >  1:163872/1‑272 (MQ=255)
                                                                                                            tGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTaaa                                                                                                                                                                                  >  2:424068/1‑301 (MQ=255)
                                                                                                            tGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGCGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCCGGAGCTGGGGCCGACAGAAGTTAAAATGCTGACCTTTTTTATGACGCATCCTGAACGCGTTTAAAGCCACAAGCAACCGTTaaa                                                                                                                                                                                  >  1:292538/1‑301 (MQ=255)
                                                                                                               tGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGCAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGAAGCATCCTGAGCGCGTGTACAGCCgagatcggaagcgcggcgg                                                                                                                                                                               >  1:595398/1‑286 (MQ=255)
                                                                                                                       ggagGGGAAGAAGAAGAGCGCGGGCGCGTCCGTGAAAGCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTggg                                                                                                                                                                       <  1:464182/300‑1 (MQ=255)
                                                                                                                        cccgacgcttttccgttttcgTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTgggg                                                                                                                                                                      <  1:241949/282‑1 (MQ=255)
                                                                                                                               agaagaagaTCGCGTGCGCGGCCTTGGAACCGGCGGGGATGACTATATCACCAATCCGTTTTCGCCGAAGGAGCGGGTGGCGCGAATCAAAGCGGGAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAAc                                                                                                                                                               <  1:342084/301‑1 (MQ=255)
                                                                                                                                    acgttcggtggctgggcttggAAACCGGTGCGGATGACTGTATCACCAAGCCGTTTTCGCCGAATGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTGCGCCAATGGCGGTGGAATAGGTGATTGAGATGCAGGGATTAAGGCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTa                                                                                                                                                          <  1:554001/281‑1 (MQ=255)
                                                                                                                                     agaTCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTACACCACGACTGGGGAACTAACGTGTaa                                                                                                                                                         >  2:113853/1‑300 (MQ=255)
                                                                                                                                       tagtcatccgcgccggtttcaaggccgcgGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTAtgt                                                                                                                                                       >  2:605996/26‑301 (MQ=255)
                                                                                                                                        agtcatccgcgccggtttcaaggccgcgGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGAGTaaggt                                                                                                                                                      >  1:91766/25‑297 (MQ=255)
                                                                                                                                          gggggtgcgGTCTTTAAACCGGGGCGGATGACTATATTACCAAGCCGGTGTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGa                                                                                                                                                    <  1:320248/295‑1 (MQ=255)
                                                                                                                                           cgTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGCCAGAAGTTAAACTGCTGCACTTTTTTATGAAGCACCCTGAGCGCGTGTCAAGCCAACCGCAGCGGGTAAAAAACGTCTGGGGAACTAACGTGTATTTgggg                                                                                                                                                   >  1:2646/1‑299 (MQ=255)
                                                                                                                                           cgTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGagatcggaagagcacacgt                                                                                                                                                   >  2:242216/1‑283 (MQ=255)
                                                                                                                                            ttccgcgccggtttcaaggccgcgGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAg                                                                                                                                                  <  2:91890/281‑1 (MQ=255)
                                                                                                                                                     cTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACgg                                                                                                                                         >  2:77789/1‑301 (MQ=255)
                                                                                                                                                     cTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAAAGTGTATGTGGAAAACCGCAAgg                                                                                                                                         >  2:148031/1‑301 (MQ=255)
                                                                                                                                                       tGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTc                                                                                                                                       >  2:105564/1‑301 (MQ=255)
                                                                                                                                                          aaCCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGaa                                                                                                                                    >  2:84714/1‑300 (MQ=255)
                                                                                                                                                                 gcgGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTTGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCcaaa                                                                                                                             >  2:197321/1‑299 (MQ=255)
                                                                                                                                                                   ggATGACTATTTCACTAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACAtt                                                                                                                           <  2:409932/301‑1 (MQ=255)
                                                                                                                                                                     aTGACTATTTCACCTAGCCTTGTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGGATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTGATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATtcg                                                                                                                         <  1:426297/301‑1 (MQ=255)
                                                                                                                                                                        aCTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATtcgtcg                                                                                                                      >  2:287802/1‑301 (MQ=255)
                                                                                                                                                                          atgtgaaCCAGGCCGTTTTCGGCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGcc                                                                                                                    <  1:605235/295‑1 (MQ=255)
                                                                                                                                                                              tCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCg                                                                                                                <  2:498147/301‑1 (MQ=255)
                                                                                                                                                                              tCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCg                                                                                                                <  2:292942/301‑1 (MQ=255)
                                                                                                                                                                                   aaGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAg                                                                                                           <  2:2646/301‑1 (MQ=255)
                                                                                                                                                                                        gTTTTCGCCGAAGGAGCTTGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGGGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGGGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTg                                                                                                      <  1:105440/301‑1 (MQ=255)
                                                                                                                                                                                            tCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGc                                                                                                  >  2:371368/1‑301 (MQ=255)
                                                                                                                                                                                               ccGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCACAGTAGTTCGCCAATGGCGGTGGCAGCGGTGACTGAGATGCTGGGAATAAGTCACGACCCGACATCTCACCGAGCGATGGCGGGCGAAGAGCCGCCGGAGATGGGGCCGACAGAATTTAAACTGCCGCTCTATTTCATGCCGCATCCTGAGCGCGTGTACCGCCCCGCGCACCTGTTAAACTAAGTCTGAGTCACTAACGTGTTAGTGTAAGAACGCACGGTCCAAGACCAAATTAGACGCCAACGTAAAacaccgaagaacc                                                                                               >  1:39413/1‑292 (MQ=255)
                                                                                                                                                                                                  aaGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCcggcg                                                                                            <  2:334191/301‑1 (MQ=255)
                                                                                                                                                                                                   gggggCTGGTGGGGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGAGGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCcggcgg                                                                                           <  1:84590/297‑1 (MQ=255)
                                                                                                                                                                                                      gggtggtggGGGGAATCAAAGCGGTAATGGGCCGTATTTCGCCAATGGCGGTGGAAGGGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCGCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCa                                                                                        <  1:413025/298‑1 (MQ=255)
                                                                                                                                                                                                        gggtggtgtggAATCAAAGCGGTATGGCTCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATg                                                                                      <  1:287400/291‑1 (MQ=255)
                                                                                                                                                                                                                      cccctacacgacgctcttccgatctccAATGGCGGTGGAAGGGGTGATTGAGATGCAGGGAGTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAg                                                                        <  1:631237/276‑1 (MQ=255)
                                                                                                                                                                                                                                         aTTTCGCCGAGTGGGGTGGAGGAGGTGGTTGAGATTCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTGATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGat                                                     <  1:423505/301‑1 (MQ=255)
                                                                                                                                                                                                                                               ccAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGagatcggaagagcacacgtctgaaa                                               >  2:632114/1‑277 (MQ=255)
                                                                                                                                                                                                                                                   tGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGATATCGTTTTTc                                           <  2:584018/301‑1 (MQ=255)
                                                                                                                                                                                                                                                       ggTGGAAGAGGTGAGTGAGATGCAGGGAGTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCGGGGGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGATATCGTTTTTCAAcc                                       <  1:77699/301‑1 (MQ=255)
                                                                                                                                                                                                                                                          ggAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGATATCGTTTTTCAACACGc                                    >  2:601645/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                   gATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGATATCGTTTTTCAACCCGCTTTTAACGc                           >  2:544152/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                       gagaTGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGATATCGTTTTTCAACCCGCTTTTAACGCCttt                       >  1:117677/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                         gaTGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCCTGGTGCAGACCGTTCGCGGTACAGGATATCGTTTTTCAACCCGCTTTTAACGCCTTGcc                     >  1:134013/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                               gtgATTAAGTCTTGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTGTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGATATCGTTTTTCAACCCGCTTTTAACGCCTTGCTCATCgg               <  1:600884/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            gACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGATATCGTTTTTCAACCCGCTTTTAACGCCTTGCTCATCGGACGCAGAGCAggg  >  2:19131/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                             
CCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGAGCCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAAGCACTGGAGCCCGGCGGGCATGACCGCATGGTGCAGACCGTGCGCGGTACAGGATATCGTTTTTCAACCCGCTTTTAACGCCTTGCTCATCGGACGCAGAGCAGGG  >  NC_000913/417275‑417859

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