Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,145,868 A→G 16.2% 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,8680AG16.2% 200.4 / 7.1 82E868G (GAG→GGG) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (25/43);  new base G (11/3);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 6.56e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACTATCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCC  >  NC_000913/2145632‑2146103
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             
tGACTTTCGTTTTAAATGGGAAGGCCGTGTGCAACGGGGAGGTGCCAGTGCCGGGATTACCTTGATTGATGCCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAgaga                                                                                                                                                                                                                                             <  2:721255/239‑1 (MQ=255)
 gACTATCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGaa                                                                                                                                                                                                                                            <  1:982746/239‑1 (MQ=255)
              tatGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCTTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAAGGAAT‑CCGTTgaggat                                                                                                                                                                                                                               >  1:626983/1‑238 (MQ=255)
              tatGGGAAGGCCGTGGGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGatgatg                                                                                                                                                                                                                              <  1:916059/239‑1 (MQ=255)
                    aaGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGGGCCGGTTTACGAACCGCAGCAACCTGCCGCACCTAGCGAACGGGCGGCGATGTCGCTTGA‑GTAAAGTGGCGGGGTGATTCAAAAGGAC‑AAGTtgatgtggacgggg                                                                                                                                                                                                                        >  2:957028/1‑230 (MQ=255)
                      ggCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGCTGAACAGTGGC‑GGATGGTTAAAAGAAAT‑CAGTTGTGGGAGCTCGccc                                                                                                                                                                                                                       >  1:1031789/1‑238 (MQ=255)
                           tgtgCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGGTTAAAGAGAAT‑ACGTTGGTGATGTTCGCCCACgggg                                                                                                                                                                                                                 >  1:693686/1‑237 (MQ=255)
                            gtgCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACATTGGC‑GGATGATTAAAGGGAAT‑CAGTTGGTGATGCTGGCCCACgggg                                                                                                                                                                                                                 >  1:1075549/1‑236 (MQ=255)
                                       aGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCa                                                                                                                                                                                                     <  1:810942/238‑1 (MQ=255)
                                         gggCCAGTGCCCGGATTTCCTTGATTGATGACAACAATCAGCAGGCGGCAGAAGTGATGGCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACg                                                                                                                                                                                                   <  1:674036/236‑1 (MQ=255)
                                         gTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGTTGATTAAAGGGAAT‑CAGTTGATGATGCTCCCCCACGGTGGCCCTTCG‑CCCCgg                                                                                                                                                                                                  >  1:688703/1‑239 (MQ=255)
                                            ccAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCCCATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGCTGATTTAAAGGAAT‑CGAGTGATGATGCTGGCCCCCGGTTGCGCTTCG‑CCCCgggt                                                                                                                                                                                                >  1:1082035/1‑236 (MQ=255)
                                              aGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGCTGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGGGTCGCTTTG‑CCACGGAcccc                                                                                                                                                                                             >  1:125787/1‑239 (MQ=255)
                                                 ggtgggTTTCCCTTGATTGATGACAACAATCATCAGTGGGCAGAAGGGAGGTCGCAGGCTGATATCGCTTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGa                                                                                                                                                                                          <  2:626983/236‑1 (MQ=255)
                                                        tttccTTGATGTATGAAAAAAATCTTCAGGCGGCAGAAGGGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGGGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGcgcg                                                                                                                                                                                    <  1:1016010/235‑1 (MQ=255)
                                                                attgatGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGGGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTGGATGAACAGTGGC‑GGATGATTTAAGGGAAT‑CCGTTGGTGATGCCCGCCCACGGGGGCGCTTCC‑CCCCCGGACCCGGAAAcgcggatattttg                                                                                                                                                                           >  2:676471/1‑230 (MQ=255)
                                                                    atgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCCAGCTGCCGCACATAGCGAACGGGGGGCGATGTCGCTTGATGAACCGTGGC‑GGCAGATTAAAGGAGAT‑CAGGTGGTGATTCACCCCCCCGGGGTCGCTCTG‑CCACGGGTCCCGGAAGGGCGAAAATTGTGGGTg                                                                                                                                                                       >  2:82601/1‑239 (MQ=255)
                                                                     tgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAAGGAAT‑CAGTTGATGATGCTCCGCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGGGGGTAATTTTTGGCTGa                                                                                                                                                                      >  1:681949/1‑239 (MQ=255)
                                                                      gaCAACCATCATCAGGCGGCAGAAGTGATGGCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGAt                                                                                                                                                                     <  1:540025/239‑1 (MQ=255)
                                                                      gaCAACAATCATCGGGGGGGAGAAGTGGTGGCGCTGGCTGATATCCGCCGTTATGCCTCCAAAAATGGTGGCCGGGGGCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCAAATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCTCGGATCCCGGAAGCGCGTAATTTGTGGCTGAt                                                                                                                                                                     <  2:894777/239‑1 (MQ=255)
                                                                      gaCAACAATCATCAGGCGGCAGAAGGGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGAt                                                                                                                                                                     <  1:512085/239‑1 (MQ=255)
                                                                           caaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCCCATAGCGAACGGGCGGGGATGTACCTTGCGGTAAAGTGGC‑GGATGGTTACAAGGAATACAGTTGATG‑TGCTCCCCCACGGGGGTCCTTCGCCCAGGGACCCCGGAA‑CGCGAGATTTTTGGcgagttttcc                                                                                                                                                                >  2:1093676/1‑230 (MQ=255)
                                                                             atcatcAGGCGGCAGAAGTGATGTCGCAGGCTGATTTCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACtt                                                                                                                                                              <  1:1010859/239‑1 (MQ=255)
                                                                                   cggcggcAGAAGTGGTGTCGCAGGCTGGTTTCGCCTGTTATGCCTCCAAAAATGGGGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGGAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGctc                                                                                                                                                        <  2:499532/238‑1 (MQ=255)
                                                                                             aGTGATGGCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTGATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCg                                                                                                                                              <  1:543967/239‑1 (MQ=255)
                                                                                                 gtGTCGCGGGGTGATATCGCCTGGTCTGCCTCCAAAAATGGGGGACGGGGCCGGGGGCAGGTTTACGAACCGCAGCAAGCTGCCGGACATAGCGAACGGGCGGCGATGGCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAgg                                                                                                                                          <  2:1027661/238‑1 (MQ=255)
                                                                                                  tggcgcctgCTTATTTCGCCGGTTCTGCCGCAAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGAACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGc                                                                                                                                         <  2:688703/231‑1 (MQ=255)
                                                                                                        aGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGat                                                                                                                                    <  1:380713/238‑1 (MQ=255)
                                                                                                         ggCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGGGTCGCTTCG‑CCACGGATCCCGGAAGCGCGGAATTTGTGGCTGATTTCACTTAAGCTCTGGGGTTGCGAAGGCGGGatta                                                                                                                                  >  1:60469/1‑239 (MQ=255)
                                                                                                            tgaatacGCCTTTTAGGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtg                                                                                                                               <  1:430830/233‑1 (MQ=255)
                                                                                                              atatCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGTGAAGGCGGGATTATtgat                                                                                                                             >  1:227406/1‑239 (MQ=255)
                                                                                                                atCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTTACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACtt                                                                                                                                                              <  1:626775/204‑1 (MQ=255)
                                                                                                                 tCGCCTGTTATGCCTACAAAAATGGGGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGA‑GATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaac                                                                                                                         <  1:812685/239‑1 (MQ=255)
                                                                                                                   gCCTGTTATGCCTCCAAAAAGGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaaca                                                                                                                        <  1:82601/239‑1 (MQ=255)
                                                                                                                    ccTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGTATCCCGGAAGCGCGGAATTTGTGGCTGATTTCACTTAAGCCCTGGAGTTGCCAAGGCGAGATTATTGATGaaaaa                                                                                                                       >  1:226648/1‑239 (MQ=255)
                                                                                                                        ttATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTTAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTTGCTTCG‑CCCCGGATCCCGGAAGCCCGGAATTTTTGGCTGATTTCCCTTTAGCCCTGGGGTTTGGCAGGGGAGGTTTTtgagaaaaaaaaaa                                                                                                                   >  2:839367/1‑228 (MQ=255)
                                                                                                                        ttATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTTACTTTAGCTCTGGGGTTGCGAAGGCGAGATTATTGATGAACAAACAt                                                                                                                   >  1:483333/1‑239 (MQ=255)
                                                                                                                        ttATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAt                                                                                                                   >  1:245847/1‑239 (MQ=255)
                                                                                                                        ttATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGAGTCGCTTCG‑CCACGGCTCCCGGAAGCGCGTAATTTTGGGCTGATTTTCTTTAAGCCCTGGAGTTGCGAAGGGCAGATTATTGATGGACAAACAt                                                                                                                   >  1:592521/1‑239 (MQ=255)
                                                                                                                         tATGCCTCCAAAAATGGGGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAtt                                                                                                                  <  1:358459/239‑1 (MQ=255)
                                                                                                                                  aaaaaTGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGGGGTA‑CCGTTGAGGATCTTCGCCCACGGTGTCGCTTCG‑CCCCGGCTCCCCGCAGCCCGTAATTTGTGGGCGATTTTCCTTAAGCTCTGGAGGTGGCGAGCGGGGAGTATTGGTGAAAAAACATTtctgccggg                                                                                                         >  2:848283/1‑232 (MQ=255)
                                                                                                                                           ggccggggccggGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTTGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTCGCTTCAGc                                                                                                  >  1:831083/1‑237 (MQ=255)
                                                                                                                                               ggggccgggGGACGGTTTACGAACCGCAGCAAGCTGCCGCGCATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCg                                                                                            <  1:1036290/239‑1 (MQ=255)
                                                                                                                                                 ggccgggTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGGGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGCAGGCGGGATTATTGATGAACAAACATTTCTTCTTAGCTTCAGCGTTCCgc                                                                                           >  1:574734/1‑237 (MQ=255)
                                                                                                                                                      ggTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGGGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGCTGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGGGCtt                                                                                      >  1:1028535/1‑238 (MQ=255)
                                                                                                                                                          aCGGTTTGCGAACCGCAGCAAGCTGCCGGACATAGAGAACGGGCGGGGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGcc                                                                                  <  1:354871/238‑1 (MQ=255)
                                                                                                                                                          aCGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGGTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGcc                                                                                  <  1:1093676/238‑1 (MQ=255)
                                                                                                                                                           cGGTTTACGAACCGCAGCAAGCTGCCGCAAATAGCGGACGGGCGGCGATGTCGCTTGATGAACAGGGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCa                                                                                 <  1:101428/238‑1 (MQ=255)
                                                                                                                                                                 tcggaCCGCAGCAAGCTGCCGCAAATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCtt                                                                          <  1:866734/235‑1 (MQ=255)
                                                                                                                                                                      ccGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGAcc                                                                      >  1:383787/1‑238 (MQ=255)
                                                                                                                                                                       cGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGaa                                                                                                                                                                                         >  1:779425/1‑122 (MQ=255)
                                                                                                                                                                       cGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGaa                                                                                                                                                                                         <  2:779425/122‑1 (MQ=255)
                                                                                                                                                                        gcagcaAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGGGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGGGAGATTATTGATGAACAAACATTTGTTCGTAGCTTCAGCGATCCGGGGCTTGGCCATACTCTTGAccgcc                                                                   >  2:504615/1‑239 (MQ=255)
                                                                                                                                                                         gagcaAGTCGCCGTACATAGCGAACCGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGAccgccg                                                                  <  1:506539/238‑1 (MQ=255)
                                                                                                                                                                                 tGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTTATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATGTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtg                                                                                                                               <  2:54534/170‑1 (MQ=255)
                                                                                                                                                                                 tGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtg                                                                                                                               >  1:54534/1‑170 (MQ=255)
                                                                                                                                                                                   ccGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGATCTTGACCGCCGGGGATTCCAc                                                        >  1:971651/1‑239 (MQ=255)
                                                                                                                                                                                    cGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGGGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAACGCGGAATTTGTGGCTGATTTAACTTAAGCTCTGGAGTTGCGAAGGGGAGACTATTGTTGAACCAACATTTCCTCGTTACTTTCGCGTCACCGCGCTTAGCCATGTTCTTGACCGCCGGGGGTTCCAcc                                                       >  2:1068945/1‑238 (MQ=255)
                                                                                                                                                                                           aGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCg                                                                                                                                              <  1:526160/145‑1 (MQ=255)
                                                                                                                                                                                           aGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCt                                                                                       >  1:857595/1‑200 (MQ=255)
                                                                                                                                                                                           aGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCt                                                                                       <  2:857595/200‑1 (MQ=255)
                                                                                                                                                                                             cGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGGTTTTTGATGAACAAACCTTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATTCTCTTGACCGCCGGGCATTCCACGAATTTTTcc                                              >  1:655822/1‑239 (MQ=255)
                                                                                                                                                                                                 cgggcgggGATTTTGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                          <  2:245847/239‑1 (MQ=255)
                                                                                                                                                                                                 cGGGCGGCGGTGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                          <  1:199057/239‑1 (MQ=255)
                                                                                                                                                                                                 cGGGCGGCGATGTCGGTTGATGAACAGTGGC‑GGATGATTAAAGGGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGGGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                          <  1:785819/239‑1 (MQ=255)
                                                                                                                                                                                                     cggcgATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTATCCATGCACTTGACCGCCGGGTATTTCCCGCAATTTTCCAGCAGGc                                       >  1:1104268/1‑238 (MQ=255)
                                                                                                                                                                                                                 ttGGTGAAAAGTGTC‑GGAGAATTAAAGAGAAT‑CAGGTGAGGAGGCGCGCCCACGGGGTCGTTTAG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGtt                          <  1:242691/239‑1 (MQ=255)
                                                                                                                                                                                                                 ttGATGAACGTTGGC‑GGTTGATTAAAGAGAAT‑CAGTTGATTATGCTCGCCCACGGTGTCGCGTCG‑CCCCGGATCCCGGAAGCGCGTAATTTGGGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGtt                          <  2:1028535/239‑1 (MQ=255)
                                                                                                                                                                                                                  gtgtgaAGAGTGGG‑GGATGAGTAAAGAGAAT‑CAGGTGGTGATGCGCGCCCGCGGGGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGGGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGGGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTg                         <  2:592521/236‑1 (MQ=255)
                                                                                                                                                                                                                   ttggaACAGTTGG‑GGGTGATTAAAGAGAAT‑CAGTTGATGATGCTCCCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGGGGGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGc                        <  2:1031789/236‑1 (MQ=255)
                                                                                                                                                                                                                       aaCAGTGGG‑GGGTGATTAAAGGGAAT‑CAGTTGATTATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGc                        <  2:254934/235‑1 (MQ=255)
                                                                                                                                                                                                                        aCAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTa                   <  1:416132/239‑1 (MQ=255)
                                                                                                                                                                                                                        aCAGGGGG‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTa                   <  1:1068945/239‑1 (MQ=255)
                                                                                                                                                                                                                              gC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAACCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCCCGGATTTTTCCCACAGGCCCCAAAAGCGGTTGCCAGTAAAGGCa             >  1:439542/1‑239 (MQ=255)
                                                                                                                                                                                                                               c‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCa             <  1:958277/238‑1 (MQ=255)
                                                                                                                                                                                                                                 ggctgtttAAAGAGATC‑CAGTTGATGATGCGGGTCCGCGGTGTCGCTTCG‑CCACGGATCCAGGAAGGGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGGTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCAt            <  2:1104268/232‑1 (MQ=255)
                                                                                                                                                                                                                                 ggATGATTAAAGGGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCAt            <  1:938700/238‑1 (MQ=255)
                                                                                                                                                                                                                                  gatgatTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGc              >  1:254934/1‑235 (MQ=255)
                                                                                                                                                                                                                                  gatgatTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCACCAGGCCGCAAAACAGGTTGCCAGTAAAGGCATaa          >  1:72498/1‑239 (MQ=255)
                                                                                                                                                                                                                                    tgatTAAAGAGAAT‑CAGTTTATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAg         <  1:949956/238‑1 (MQ=255)
                                                                                                                                                                                                                                      atTAAAGAGAAT‑CAGTTGATGATGCTCGCCCTCGGTTTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAt      <  1:264319/239‑1 (MQ=255)
                                                                                                                                                                                                                                         aaaGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCCGGAGTTGCGAAGGCGAGATTATTGATGCACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGCATTCCTCGAATTTTTCCAGTAGGCCGCAAAATCTGTTGCCAGGTAAGGCATAAGCAtcc    >  1:75227/1‑237 (MQ=255)
                                                                                                                                                                                                                                          aaGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCTTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATctcc  <  1:782678/239‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             
TGACTATCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGATGCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCC  >  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: ‑