Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,145,868 A→G 14.8% E868G (GAG→GGG)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,145,8680AG14.8% 229.3 / 7.0 95E868G (GAG→GGG) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (33/47);  new base G (12/3);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 9.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACTATCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCC  >  NC_000913/2145632‑2146103
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
tGACTTTCGTTTTAAATGGGAAGGCCGTGTGCAACGGGGAGGTGCCAGTGCCGGGATTACCTTGATTGATGCCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAgaga                                                                                                                                                                                                                                              <  2:721255/239‑1 (MQ=255)
 gACTATCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGaa                                                                                                                                                                                                                                             <  1:982746/239‑1 (MQ=255)
   ggtataagaggcaggggAAGGCCGTGTGCACCGGGTAGGTGCCAGTTGCGGGATTCCCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATc                                                                                                                                                                                                                                          <  1:1005372/225‑1 (MQ=255)
              tatGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCTTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAAGGAATCCGTTgaggat                                                                                                                                                                                                                                >  1:626983/1‑238 (MQ=255)
              tatGGGAAGGCCGTGGGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGatgatg                                                                                                                                                                                                                              <  1:916059/239‑1 (MQ=255)
                    aaGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGGGCCGGTTTACGAACCGCAGCAACCTGCCGCACCTAGCGAACGGGCGGCGATGTCGCTTGAGTAAAGTGGCGGGGTGATTCAAAAGGACAAGTtgatgtggacgggg                                                                                                                                                                                                                         >  2:957028/1‑230 (MQ=255)
                      ggCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGCTGAACAGTGGCGGATGGTTAAAAGAAATCAGTTGTGGGAGCTCGccc                                                                                                                                                                                                                       >  1:1031789/1‑238 (MQ=255)
                        ccGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGCACCGCAGCAAGCTGCCGCACATAGCGAACGGCGGGCGGGTGCGCTTGATGAAACGGGCCGGGGGATTACAGAGAAACAgaggttggtgccaccccacg                                                                                                                                                                                                                     >  2:199057/1‑220 (MQ=255)
                           tgtgCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGGTGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAAACGCAGCAAGCTGCCGCACCTAGCGAACGGGCGGCCGACTCCCTTGGTGAACCGGGGCGGGCGACTCAAAAGAATCacattatggaacgcgcccccaggg                                                                                                                                                                                                                  >  2:430830/1‑216 (MQ=255)
                           tgtgCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGGTTAAAGAGAATACGTTGGTGATGTTCGCCCACgggg                                                                                                                                                                                                                 >  1:693686/1‑237 (MQ=255)
                            gtgCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACATTGGCGGATGATTAAAGGGAATCAGTTGGTGATGCTGGCCCACgggg                                                                                                                                                                                                                 >  1:1075549/1‑236 (MQ=255)
                                       aGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCa                                                                                                                                                                                                     <  1:810942/238‑1 (MQ=255)
                                        tggtccagggcgcggATTACCTTGATTGGTGACAACAATCATCAAGCGGCAGAAGGGATGCCGCCGGGTGATATCCGCTGTTATTCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCggctcg                                                                                                                                                                                                    <  2:216196/226‑6 (MQ=255)
                                         gggCCAGTGCCCGGATTTCCTTGATTGATGACAACAATCAGCAGGCGGCAGAAGTGATGGCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACg                                                                                                                                                                                                   <  1:674036/236‑1 (MQ=255)
                                         gTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGTTGATTAAAGGGAATCAGTTGATGATGCTCCCCCACGGTGGCCCTTCGCCCCgg                                                                                                                                                                                                  >  1:688703/1‑239 (MQ=255)
                                            ccAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCCCATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGCTGATTTAAAGGAATCGAGTGATGATGCTGGCCCCCGGTTGCGCTTCGCCCCgggt                                                                                                                                                                                                >  1:1082035/1‑236 (MQ=255)
                                              aGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGCTGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGGGTCGCTTTGCCACGGAcccc                                                                                                                                                                                             >  1:125787/1‑239 (MQ=255)
                                                 ggtgggTTTCCCTTGATTGATGACAACAATCATCAGTGGGCAGAAGGGAGGTCGCAGGCTGATATCGCTTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGa                                                                                                                                                                                          <  2:626983/236‑1 (MQ=255)
                                                        tttccTTGATGTATGAAAAAAATCTTCAGGCGGCAGAAGGGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGGGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGcgcg                                                                                                                                                                                    <  1:1016010/235‑1 (MQ=255)
                                                                attgatGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAACTGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCCCCTAGCGAACGGGCGGCGATGTCGCTTGATGAACGGGGGCGGGTGATTAAAGAGGAACAGTTGGTGATCGCTGCCCAACGGGTGGCGTAGCCCGGGACCCCGGAAgcaggaaatctgag                                                                                                                                                                           >  2:416132/1‑227 (MQ=255)
                                                                attgatGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGGGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTGGATGAACAGTGGCGGATGATTTAAGGGAATCCGTTGGTGATGCCCGCCCACGGGGGCGCTTCCCCCCCGGACCCGGAAAcgcggatattttg                                                                                                                                                                           >  2:676471/1‑230 (MQ=255)
                                                                    atgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCCAGCTGCCGCACATAGCGAACGGGGGGCGATGTCGCTTGATGAACCGTGGCGGCAGATTAAAGGAGATCAGGTGGTGATTCACCCCCCCGGGGTCGCTCTGCCACGGGTCCCGGAAGGGCGAAAATTGTGGGTg                                                                                                                                                                       >  2:82601/1‑239 (MQ=255)
                                                                     tgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAAGGAATCAGTTGATGATGCTCCGCCACGGTGTCGCTTCGCCACGGATCCCGGAAGGGGGTAATTTTTGGCTGa                                                                                                                                                                      >  1:681949/1‑239 (MQ=255)
                                                                     tgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGGCGCTTGATGAACAGTGGCGGACGATTAAAGAGAACTCGTTGATGATGCTCGCCCACGGGTGCGCTTCGCCACCGGGCCCgggaaccggagattttgggccgg                                                                                                                                                                      >  2:958277/1‑218 (MQ=255)
                                                                      gaCAACCATCATCAGGCGGCAGAAGTGATGGCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGAt                                                                                                                                                                     <  1:540025/239‑1 (MQ=255)
                                                                      gaCAACAATCATCGGGGGGGAGAAGTGGTGGCGCTGGCTGATATCCGCCGTTATGCCTCCAAAAATGGTGGCCGGGGGCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCAAATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCTCGGATCCCGGAAGCGCGTAATTTGTGGCTGAt                                                                                                                                                                     <  2:894777/239‑1 (MQ=255)
                                                                      gaCAACAATCATCAGGCGGCAGAAGGGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGAt                                                                                                                                                                     <  1:512085/239‑1 (MQ=255)
                                                                           caaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCCCATAGCGAACGGGCGGGGATGTACCTTGCGGTAAAGTGGCGGATGGTTACAAGGAATACAGTTGATGTGCTCCCCCACGGGGGTCCTTCGCCCAGGGACCCCGGAACGCGAGATTTTTGGcgagttttcc                                                                                                                                                                >  2:1093676/1‑230 (MQ=255)
                                                                             atcatcAGGCGGCAGAAGTGATGTCGCAGGCTGATTTCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACtt                                                                                                                                                              <  1:1010859/239‑1 (MQ=255)
                                                                                   cggcggcAGAAGTGGTGTCGCAGGCTGGTTTCGCCTGTTATGCCTCCAAAAATGGGGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGGAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGctc                                                                                                                                                        <  2:499532/238‑1 (MQ=255)
                                                                                             aGTGATGGCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTGATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCg                                                                                                                                              <  1:543967/239‑1 (MQ=255)
                                                                                               tGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGAGTAAAGAGAATCAGTAGAAGAGTCTCGCCCACGGTGTGGCCTCGCCACCGGTGCCGGAAACGCGAGCTTTGGTGGTGTTTTCACCTaaagcctggcggtgggga                                                                                                                                            >  2:949956/1‑223 (MQ=255)
                                                                                                 gtGTCGCGGGGTGATATCGCCTGGTCTGCCTCCAAAAATGGGGGACGGGGCCGGGGGCAGGTTTACGAACCGCAGCAAGCTGCCGGACATAGCGAACGGGCGGCGATGGCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAgg                                                                                                                                          <  2:1027661/238‑1 (MQ=255)
                                                                                                  tgtcggagcgggttaccGCTTGTTAGCCCTCAAAAAATGGTGGGCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGACGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACGTAAGCTCTGGAGTTGCGAAGGc                                                                                                                                         <  2:1016533/223‑1 (MQ=255)
                                                                                                  tggcgcctgCTTATTTCGCCGGTTCTGCCGCAAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGAACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGc                                                                                                                                         <  2:688703/231‑1 (MQ=255)
                                                                                                  tGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAAGGGATACGGTGATGGTGGCCGCCACCCGGGTCCGTTCGCCACGAGTCCCGGAAGCGCGTAAATTTTGGGTGATATCCATTaaaaccaggggggtggaaagg                                                                                                                                         >  2:526773/1‑220 (MQ=255)
                                                                                                  gatttgtataagggcccGCCTGTTATGCCTCCAAAAAGGGTGGCCGTGGCCGGGTGGCGGGTTACGAACCGCAGCAAGCTGCCGCATATAGGGAACTGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGc                                                                                                                                         <  2:337400/223‑1 (MQ=255)
                                                                                                        aGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGat                                                                                                                                    <  1:380713/238‑1 (MQ=255)
                                                                                                         ggCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGGGTCGCTTCGCCACGGATCCCGGAAGCGCGGAATTTGTGGCTGATTTCACTTAAGCTCTGGGGTTGCGAAGGCGGGatta                                                                                                                                  >  1:60469/1‑239 (MQ=255)
                                                                                                            tgaatacGCCTTTTAGGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtg                                                                                                                               <  1:430830/233‑1 (MQ=255)
                                                                                                              atatCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGTGAAGGCGGGATTATtgat                                                                                                                             >  1:227406/1‑239 (MQ=255)
                                                                                                                atCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTTACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACtt                                                                                                                                                              <  1:626775/204‑1 (MQ=255)
                                                                                                                 tCGCCTGTTATGCCTACAAAAATGGGGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGAGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaac                                                                                                                         <  1:812685/239‑1 (MQ=255)
                                                                                                                   gCCTGTTATGCCTCCAAAAAGGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaaca                                                                                                                        <  1:82601/239‑1 (MQ=255)
                                                                                                                    ccTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGTATCCCGGAAGCGCGGAATTTGTGGCTGATTTCACTTAAGCCCTGGAGTTGCCAAGGCGAGATTATTGATGaaaaa                                                                                                                       >  1:226648/1‑239 (MQ=255)
                                                                                                                        ttATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTTAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTTGCTTCGCCCCGGATCCCGGAAGCCCGGAATTTTTGGCTGATTTCCCTTTAGCCCTGGGGTTTGGCAGGGGAGGTTTTtgagaaaaaaaaaa                                                                                                                   >  2:839367/1‑228 (MQ=255)
                                                                                                                        ttATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTTACTTTAGCTCTGGGGTTGCGAAGGCGAGATTATTGATGAACAAACAt                                                                                                                   >  1:483333/1‑239 (MQ=255)
                                                                                                                        ttATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAt                                                                                                                   >  1:245847/1‑239 (MQ=255)
                                                                                                                        ttATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGAGTCGCTTCGCCACGGCTCCCGGAAGCGCGTAATTTTGGGCTGATTTTCTTTAAGCCCTGGAGTTGCGAAGGGCAGATTATTGATGGACAAACAt                                                                                                                   >  1:592521/1‑239 (MQ=255)
                                                                                                                         tATGCCTCCAAAAATGGGGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAtt                                                                                                                  <  1:358459/239‑1 (MQ=255)
                                                                                                                               tCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTTTtgatggttccttttcaacatt                                                                                                            >  1:337400/1‑223 (MQ=255)
                                                                                                                                  aaaaaTGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGGGGTACCGTTGAGGATCTTCGCCCACGGTGTCGCTTCGCCCCGGCTCCCCGCAGCCCGTAATTTGTGGGCGATTTTCCTTAAGCTCTGGAGGTGGCGAGCGGGGAGTATTGGTGAAAAAACATTtctgccggg                                                                                                         >  2:848283/1‑232 (MQ=255)
                                                                                                                                  aaaaaTGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGGGTCGCTTCGCCACGGAACCCGGAAGCGCGTAATTTGTGGCTGATTTTACTTAAGCTCTGGGGTGGGCAAGGGGGGATTCTTGAGGaacaaaaaattttggggag                                                                                                         >  2:521304/1‑226 (MQ=255)
                                                                                                                                           ggccggggccggGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTTGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTCGCTTCAGc                                                                                                  >  1:831083/1‑237 (MQ=255)
                                                                                                                                               ggggccgggGGACGGTTTACGAACCGCAGCAAGCTGCCGCGCATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCg                                                                                            <  1:1036290/239‑1 (MQ=255)
                                                                                                                                                 ggccgggTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGGGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGCAGGCGGGATTATTGATGAACAAACATTTCTTCTTAGCTTCAGCGTTCCgc                                                                                           >  1:574734/1‑237 (MQ=255)
                                                                                                                                                      ggTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGGGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGCTGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGGGCtt                                                                                      >  1:1028535/1‑238 (MQ=255)
                                                                                                                                                          aCGGTTTGCGAACCGCAGCAAGCTGCCGGACATAGAGAACGGGCGGGGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGcc                                                                                  <  1:354871/238‑1 (MQ=255)
                                                                                                                                                          aCGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGGTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGcc                                                                                  <  1:1093676/238‑1 (MQ=255)
                                                                                                                                                           cGGTTTACGAACCGCAGCAAGCTGCCGCAAATAGCGGACGGGCGGCGATGTCGCTTGATGAACAGGGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCa                                                                                 <  1:101428/238‑1 (MQ=255)
                                                                                                                                                                 tcggaCCGCAGCAAGCTGCCGCAAATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCtt                                                                          <  1:866734/235‑1 (MQ=255)
                                                                                                                                                                      ccGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGAcc                                                                      >  1:383787/1‑238 (MQ=255)
                                                                                                                                                                       cGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGaa                                                                                                                                                                                         <  2:779425/122‑1 (MQ=255)
                                                                                                                                                                       cGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGaa                                                                                                                                                                                         >  1:779425/1‑122 (MQ=255)
                                                                                                                                                                        gcagcaAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGGGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGGGAGATTATTGATGAACAAACATTTGTTCGTAGCTTCAGCGATCCGGGGCTTGGCCATACTCTTGAccgcc                                                                   >  2:504615/1‑239 (MQ=255)
                                                                                                                                                                         gagcaAGTCGCCGTACATAGCGAACCGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGAccgccg                                                                  <  1:506539/238‑1 (MQ=255)
                                                                                                                                                                                 tGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTTATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATGTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtg                                                                                                                               <  2:54534/170‑1 (MQ=255)
                                                                                                                                                                                 tGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtg                                                                                                                               >  1:54534/1‑170 (MQ=255)
                                                                                                                                                                                   ccGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGATCTTGACCGCCGGGGATTCCAc                                                        >  1:971651/1‑239 (MQ=255)
                                                                                                                                                                                    cGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGGGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAACGCGGAATTTGTGGCTGATTTAACTTAAGCTCTGGAGTTGCGAAGGGGAGACTATTGTTGAACCAACATTTCCTCGTTACTTTCGCGTCACCGCGCTTAGCCATGTTCTTGACCGCCGGGGGTTCCAcc                                                       >  2:1068945/1‑238 (MQ=255)
                                                                                                                                                                                      cccgttcgctatgtgcggcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGaa                                                     >  1:4302/15‑239 (MQ=255)
                                                                                                                                                                                           aGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCg                                                                                                                                              <  1:526160/145‑1 (MQ=255)
                                                                                                                                                                                           aGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCt                                                                                       <  2:857595/200‑1 (MQ=255)
                                                                                                                                                                                           aGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCt                                                                                       >  1:857595/1‑200 (MQ=255)
                                                                                                                                                                                             cGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGGTTTTTGATGAACAAACCTTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATTCTCTTGACCGCCGGGCATTCCACGAATTTTTcc                                              >  1:655822/1‑239 (MQ=255)
                                                                                                                                                                                                 cgggcgggGATTTTGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                          <  2:245847/239‑1 (MQ=255)
                                                                                                                                                                                                 cGGGCGGCGGTGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                          <  1:199057/239‑1 (MQ=255)
                                                                                                                                                                                                 cGGGCGGCGATGTCGGTTGATGAACAGTGGCGGATGATTAAAGGGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGGGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                          <  1:785819/239‑1 (MQ=255)
                                                                                                                                                                                                     cggcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTATCCATGCACTTGACCGCCGGGTATTTCCCGCAATTTTCCAGCAGGc                                       >  1:1104268/1‑238 (MQ=255)
                                                                                                                                                                                                                 ttGGTGAAAAGTGTCGGAGAATTAAAGAGAATCAGGTGAGGAGGCGCGCCCACGGGGTCGTTTAGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGtt                          <  1:242691/239‑1 (MQ=255)
                                                                                                                                                                                                                 ttGATGAACGTTGGCGGTTGATTAAAGAGAATCAGTTGATTATGCTCGCCCACGGTGTCGCGTCGCCCCGGATCCCGGAAGCGCGTAATTTGGGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGtt                          <  2:1028535/239‑1 (MQ=255)
                                                                                                                                                                                                                  gtgtgaAGAGTGGGGGATGAGTAAAGAGAATCAGGTGGTGATGCGCGCCCGCGGGGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGGGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGGGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTg                         <  2:592521/236‑1 (MQ=255)
                                                                                                                                                                                                                   ttggaACAGTTGGGGGTGATTAAAGAGAATCAGTTGATGATGCTCCCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGGGGGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGc                        <  2:1031789/236‑1 (MQ=255)
                                                                                                                                                                                                                       aaCAGTGGGGGGTGATTAAAGGGAATCAGTTGATTATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGc                        <  2:254934/235‑1 (MQ=255)
                                                                                                                                                                                                                        aCAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTa                   <  1:416132/239‑1 (MQ=255)
                                                                                                                                                                                                                        aCAGGGGGGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTa                   <  1:1068945/239‑1 (MQ=255)
                                                                                                                                                                                                                              gCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAACCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCCCGGATTTTTCCCACAGGCCCCAAAAGCGGTTGCCAGTAAAGGCa             >  1:439542/1‑239 (MQ=255)
                                                                                                                                                                                                                               cGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCa             <  1:958277/238‑1 (MQ=255)
                                                                                                                                                                                                                                 ggctgtttAAAGAGATCCAGTTGATGATGCGGGTCCGCGGTGTCGCTTCGCCACGGATCCAGGAAGGGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGGTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCAt            <  2:1104268/232‑1 (MQ=255)
                                                                                                                                                                                                                                 ggATGATTAAAGGGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCAt            <  1:938700/238‑1 (MQ=255)
                                                                                                                                                                                                                                  gatgatTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGc              >  1:254934/1‑235 (MQ=255)
                                                                                                                                                                                                                                  gatgatTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCACCAGGCCGCAAAACAGGTTGCCAGTAAAGGCATaa          >  1:72498/1‑239 (MQ=255)
                                                                                                                                                                                                                                    tgatTAAAGAGAATCAGTTTATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAg         <  1:949956/238‑1 (MQ=255)
                                                                                                                                                                                                                                      atTAAAGAGAATCAGTTGATGATGCTCGCCCTCGGTTTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAt      <  1:264319/239‑1 (MQ=255)
                                                                                                                                                                                                                                         aaaGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCCGGAGTTGCGAAGGCGAGATTATTGATGCACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGCATTCCTCGAATTTTTCCAGTAGGCCGCAAAATCTGTTGCCAGGTAAGGCATAAGCAtcc    >  1:75227/1‑237 (MQ=255)
                                                                                                                                                                                                                                          aaGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCTTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATctcc  <  1:782678/239‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
TGACTATCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCC  >  NC_000913/2145632‑2146103

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