Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,563,546 T→G 7.6% intergenic (+28/+241) glpD → / ← yzgL sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding/pseudogene, periplasmic solute binding protein homology

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,563,5460TG 7.6% 512.9 / 26.5 147intergenic (+28/+241)glpD/yzgLsn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding/pseudogene, periplasmic solute binding protein homology
Reads supporting (aligned to +/- strand):  ref base T (57/78);  new base G (6/6);  total (63/84)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

Only 100 of 148 total aligned reads displayed.

ACGGCAGCAACAGCGAGCTGCTGCTCGGCAATGCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAG  >  NC_000913/3563274‑3563807
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                     
              gAGCTGCTGCTCGGCAATGCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACANAACAAGGCATGTGGCTAAATCCGGATCAACAATCTCGTGTGAGTCAGTGGCGGGGTGAGTATACGCAGCATAGGNGATAGATGGAGTCGGAAAGTAAGGTAGTGGNNAT‑AGGTAAGATAGCAATCTGGGCTGAGAc                                                                                                                                                                                                                                                       >  3:58355/1‑275 (MQ=255)
                                gCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACaaaa                                                                                                                                                                                                                                      <  2:223577/274‑1 (MQ=255)
                                gCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACaaaa                                                                                                                                                                                                                                      <  4:342739/274‑1 (MQ=255)
                                           aaGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgcg                                                                                                                                                                                                                          <  2:467941/275‑1 (MQ=255)
                                           aaGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgcg                                                                                                                                                                                                                          <  2:367630/275‑1 (MQ=255)
                                           aaGCGATCTCGGAGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgcg                                                                                                                                                                                                                          <  2:265906/275‑1 (MQ=255)
                                                      gggAAGATTTCGGGAATGAGTTGTTCGAGGCGGGGTTGAAATGCCTGGTGGCTACCGAAGGGGTTCGGTGCGCCGGCGACGCCCTGGGGGGTCGCAAAAAACAAGGCATGTGGCTAAATGCGGATCAACAAGCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                                <  4:287652/274‑1 (MQ=255)
                                                       ggAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                                <  1:44450/273‑1 (MQ=255)
                                                               tCGGTCATTAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATAACGCATGGTTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                                                                                                                                       <  1:132446/274‑1 (MQ=255)
                                                               tCGGTACTGAGTTCGCCAAAGCGGAGCTGAAATACCTGGGGGATAACGAATGGGTCCGCCGAGCCGACGACGCCCTGGGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCGCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                                                                                                                                       <  4:517020/274‑1 (MQ=255)
                                                                   tcgtggggtacaCGGAGCGGAGGTGAAATAGCGGGTGGGTAACGAAGGGGACCGCCGCGCCGGCGACGCCCTGTGGCGTCGCAACAAACAAGGCATGTGGCTAAATGCGGACAAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTAGACGCAGCAGAGGTTATCGCTGGCGTCGTAAAGTAACGTA‑AGGTGATCAGGTCAGATTGCAAGCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGt                                                                                                                                                                                                  <  4:129371/264‑1 (MQ=255)
                                                                    caggAGTTTTACGAAGCGGGGGCGAAAGGCATGGTGGTTCACGAATGGGTCCGGCGGGCCGAGGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGGGTCAGTGGCTGGTGGAGTAGACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGt                                                                                                                                                                                                  <  4:420216/271‑1 (MQ=255)
                                                                       gAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTTAGTGGCTGGTGGAGTATACGTATCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTTATTAGGTCAGATTTCACTCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTAt                                                                                                                                                                                              >  1:11131/1‑275 (MQ=255)
                                                                           tCTAGGAAGCGGACCTGAAATACCTGGTGGGCGACGGATGGGTCCGCCGCGCCGACGACGCCCTGTGGGGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGGGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCa                                                                                                                                                                                           <  4:370800/274‑1 (MQ=255)
                                                                           tCTACGAAGCGGGGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCTGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCa                                                                                                                                                                                           <  2:2415/274‑1 (MQ=255)
                                                                               cGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAatata                                                                                                                                                                                      <  2:100880/275‑1 (MQ=255)
                                                                                  aGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATAtt                                                                                                                                                                                    <  2:175490/274‑1 (MQ=255)
                                                                                          tGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCa                                                                                                                                                                           <  2:83219/275‑1 (MQ=255)
                                                                                                              gAATGGGTCCGCCGTCCCGACGCCGCCCTGTGGCCTCGCACAACACAAGGCATGTGGCTAAATGCGTATAGACAATCTCGTGTGAGTCGGTGGCTGGTGGAGTATACGCAGCAGAGGTTTTCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATaa                                                                                                                                                       <  2:510914/275‑1 (MQ=38)
                                                                                                              agaggggTCCGCCGGGCCCTAGACGCCCTGGGGGGTCGCGAAAAACAGGGCATGTGGGTAAATGCGGATTAACAATTNGGGGTGAGTCAGTGGCTGGTGGAGTATACGCAGCGGGGGTTATCGCTGGAGTCGGAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATaa                                                                                                                                                       <  3:124568/271‑1 (MQ=37)
                                                                                                                     cccnnggcgccgACGAGGCCTGGTGGCGTCGCACAAAACAAGGGATGTGGCTANATGCGGATCAACAGTCTCGTGTGGGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGAGGACAAACACAAAACTGCCTGATGCGCTACg                                                                                                                                                                                                                     <  4:542854/200‑1 (MQ=255)
                                                                                                                                        ccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCt                                                                                                                                                                                                                                                            >  7:1110824/1‑148 (MQ=255)
                                                                                                                                                  cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatg                                                                                                                                                                                                                                                  <  8:1110824/148‑1 (MQ=255)
                                                                                                                                                    cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGAc                                                                                                                                                                                                                                                >  5:34874/1‑148 (MQ=255)
                                                                                                                                                      caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATc                                                                                                                                                                                                                                                                  >  7:128038/1‑128 (MQ=255)
                                                                                                                                                      caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATc                                                                                                                                                                                                                                                                  <  8:128038/128‑1 (MQ=255)
                                                                                                                                                      caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTg                                                                                                                                                                                                                                                                >  7:119256/1‑130 (MQ=255)
                                                                                                                                                      caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTg                                                                                                                                                                                                                                                                <  8:119256/130‑1 (MQ=255)
                                                                                                                                                      caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtg                                                                                                                                                                                                                                                     <  7:906686/141‑1 (MQ=255)
                                                                                                                                                          aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtga                                                                                                                                                                                                                                                    <  6:15833/138‑1 (MQ=255)
                                                                                                                                                                           aTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgcgc                                                                                                                                                                                                                         >  6:536833/1‑148 (MQ=255)
                                                                                                                                                                                gATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAAtt                                                                                     >  2:72816/1‑275 (MQ=255)
                                                                                                                                                                                  tCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCtt                                                                                                                                                                                                                  >  6:33956/1‑148 (MQ=255)
                                                                                                                                                                                     acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAAt                                                                                                                                                                                                                                                                   >  7:528077/1‑96 (MQ=255)
                                                                                                                                                                                     acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAAt                                                                                                                                                                                                                                                                   <  8:528077/96‑1 (MQ=255)
                                                                                                                                                                                     acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTgg                                                                                                                                                                                                                                                               <  5:351168/100‑1 (MQ=255)
                                                                                                                                                                                     acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTgg                                                                                                                                                                                                                                                               >  6:351168/1‑100 (MQ=255)
                                                                                                                                                                                     acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGa                                                                                                                                                                                                                              >  7:249933/1‑133 (MQ=255)
                                                                                                                                                                                     acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGa                                                                                                                                                                                                                              <  8:249933/133‑1 (MQ=255)
                                                                                                                                                                                        aTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTc                                                                                                                                                                                                                                                                      <  6:34165/90‑1 (MQ=255)
                                                                                                                                                                                              tgtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGNTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACNTTCACGTCGAATCCGGCATGAACAAAGAGCACTTTGTTAAAAATCNANACTACTTTTTAATTTAAGGAAATaacac                                                                       >  3:33493/1‑271 (MQ=255)
                                                                                                                                                                                              tgtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCg                                                                        >  2:43472/1‑274 (MQ=255)
                                                                                                                                                                                                 gAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTgg                                                                                                                                                                                                   <  5:33956/148‑1 (MQ=255)
                                                                                                                                                                                                  aGTCAGTGGCTGCGGGAGTATAAGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATAAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAg                                                                   <  2:13590/275‑1 (MQ=255)
                                                                                                                                                                                                   gTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCt                                                                                                                                                                                                         >  6:143794/1‑140 (MQ=255)
                                                                                                                                                                                                   gTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCt                                                                                                                                                                                                         <  5:143794/140‑1 (MQ=255)
                                                                                                                                                                                                    tCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAac                                                                                                                                                                                                                                            >  5:560566/1‑104 (MQ=255)
                                                                                                                                                                                                    tCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCCCATCCGGCCTTCACAAAGCTCACTTTGTCAAAAATCTAACATACTTTTTAATTCCGGGAAATACCGCAAAGc                                                                  >  4:366545/1‑274 (MQ=255)
                                                                                                                                                                                                     cAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                            <  8:190559/135‑1 (MQ=255)
                                                                                                                                                                                                     cAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                            >  7:190559/1‑135 (MQ=255)
                                                                                                                                                                                                     cAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                                                                                                                                               >  7:142344/1‑148 (MQ=255)
                                                                                                                                                                                                       gTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGc                                                                                                                                                                                                           >  4:467860/1‑134 (MQ=255)
                                                                                                                                                                                                        tGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTgag                                                                                                                                                                                                                                                         <  7:197163/87‑1 (MQ=255)
                                                                                                                                                                                                        tGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTgag                                                                                                                                                                                                                                                         >  8:197163/1‑87 (MQ=255)
                                                                                                                                                                                                          gCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                                                                                                                                               >  5:611114/1‑143 (MQ=255)
                                                                                                                                                                                                          gCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                                                                                                                                               <  6:611114/143‑1 (MQ=255)
                                                                                                                                                                                                             ggtggAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTCTCAGGCCTACGTGGTTNATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATTAACAAATCGNACTTTGTCCAACAACCAACCTACTTTTTAANNTATGGAATTTCCCCACACCCCaaccac                                                          >  3:69779/1‑268 (MQ=255)
                                                                                                                                                                                                               tggAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGcc                                                                                                                                                                                                                                 >  8:455430/1‑104 (MQ=255)
                                                                                                                                                                                                               tggAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGcc                                                                                                                                                                                                                                 <  7:455430/104‑1 (MQ=255)
                                                                                                                                                                                                                   gTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATAcc                                                   <  1:58355/274‑1 (MQ=255)
                                                                                                                                                                                                                    tataCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCaa                                                                                                                                                                                                                                                                    <  6:929044/64‑1 (MQ=255)
                                                                                                                                                                                                                    tataCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCaa                                                                                                                                                                                                                                                                    >  5:929044/1‑64 (MQ=255)
                                                                                                                                                                                                                        cGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCa                                                                                                                                                                                           >  7:399656/1‑133 (MQ=255)
                                                                                                                                                                                                                        cGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCa                                                                                                                                                                                           <  8:399656/133‑1 (MQ=255)
                                                                                                                                                                                                                            gcagAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGCCAGATTGAAATCTga                                                                                                                                                                                                                                                               >  8:897701/1‑60 (MQ=37)
                                                                                                                                                                                                                            gcagAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGCCAGATTGAAATCTga                                                                                                                                                                                                                                                               <  7:897701/61‑2 (MQ=37)
                                                                                                                                                                                                                               gagGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAAtt                                                                                                                                                                               >  6:609220/1‑138 (MQ=255)
                                                                                                                                                                                                                               gagGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAAtt                                                                                                                                                                               <  5:609220/138‑1 (MQ=255)
                                                                                                                                                                                                                                agGTTATCGCTGGCGTCGGAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCtt                                                                                                                                                                    <  8:142344/148‑1 (MQ=255)
                                                                                                                                                                                                                                  gTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACCAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACt                                    >  2:21420/1‑274 (MQ=255)
                                                                                                                                                                                                                                    tATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                                                                                                                                               <  5:292666/117‑1 (MQ=255)
                                                                                                                                                                                                                                    tATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAg                                 <  4:43471/275‑1 (MQ=255)
                                                                                                                                                                                                                                     aTCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAACCACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAg                                 >  2:105122/1‑274 (MQ=255)
                                                                                                                                                                                                                                     aTCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTc                                                                                                                                           >  2:470055/1‑168 (MQ=255)
                                                                                                                                                                                                                                      tCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGc                                                                                                                                                              >  5:517258/1‑148 (MQ=255)
                                                                                                                                                                                                                                         cTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCAt                                                                                                                                       <  4:470054/168‑1 (MQ=255)
                                                                                                                                                                                                                                            gCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAAtt                                                                                                                                                                               >  1:18295/1‑125 (MQ=255)
                                                                                                                                                                                                                                            gCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCATTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGcc                                                                 >  2:332909/1‑235 (MQ=255)
                                                                                                                                                                                                                                              gtcgtAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAg                                                                                                                                                      <  6:517258/148‑1 (MQ=255)
                                                                                                                                                                                                                                                cgtAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCa                                                                                                                                                                           <  3:18295/125‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                              <  5:1075819/84‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttAACGTA‑AGGTGATCAGGCCAGATTGAAATCTga                                                                                                                                                                                                                                                               >  8:890305/1‑34 (MQ=2)
                                                                                                                                                                                                                                                          cGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGcc                                                                                                                                                             <  8:943220/129‑1 (MQ=255)
                                                                                                                                                                                                                                                               aGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACaaaa                                                                                                                                                                                                                                      >  4:392424/1‑52 (MQ=255)
                                                                                                                                                                                                                                                                 gTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATAtt                                                                                                                                                                                    <  5:51148/100‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                                                                                                                                               >  7:684809/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                  tGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                                                                                                                                               <  8:684809/88‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ggtacAGGTCCGATTTCAACCTTGCCTGGGACGGATGACAAACACAAAATGTCCTGATGGCTTACTTTTAGCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAATGCGCACTTTGTCAAAAATCTACCCTACTTTTTAATTTAGGGAATTACCGCAAAGCCAACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGaaaaa     <  4:234362/270‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGc                                                                                                                                                                                                                                  <  2:392425/52‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aTCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGcc                                                                 >  1:135080/1‑212 (MQ=255)
                                                                                                                                                                                                                                                                    aTCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAg  >  1:66563/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                     tCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                            <  6:76229/72‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                            >  5:76229/1‑72 (MQ=255)
                                                                                                                                                                                                                                                                       aGGCCAGATTGAAATCTGACCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCa                                                                                                                                                                                           <  7:532720/87‑1 (MQ=18)
                                                                                                                                                                                                                                                                        ggTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCa                                                                                                                                                                                           <  7:977423/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ggTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCa                                                                                                                                                                                           >  8:977423/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                         gTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCa                                                                                                                                                                           <  8:634165/101‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACg                                                                                                                                                                                                                     <  7:999771/58‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGa                                                                                                                                                                                                                              <  7:806783/48‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGa                                                                                                                                                                                                                              >  6:454081/1‑48 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                     
ACGGCAGCAACAGCGAGCTGCTGCTCGGCAATGCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTA‑AGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAG  >  NC_000913/3563274‑3563807

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑