Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,956,003 G→A D225N (GAT→AAT)  ilvA → threonine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,956,0030GA97.8% 277.8 / NA 91D225N (GAT→AAT) ilvAthreonine deaminase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (34/55);  total (36/55)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCG  >  NC_000913/3955775‑3956249
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
cagGGGTTCACCTGTGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTGATGCCGCAAATCAAAAGGATAGCCGTAGAAGCGGGAGAATCCGCATGCCTGAAAGCCGCGCTGGAGGCGGGGCATCCGGGTGATCt                                                                                                                                                                                                                                   >  2:395687/1‑251 (MQ=255)
   gggTTCCCCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGcc                                                                                                                                                                                                                                <  1:110707/251‑1 (MQ=255)
      ttCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAAGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTTAAGCCGCAAATCACAGTGATCGCCGTAGAAGCGGGAAACTCCGCCAGCCTGAAAGCGGCGCCGAATGCGGGTCATCCGGTTGCTCTTCcgcg                                                                                                                                                                                                                             >  1:151834/1‑251 (MQ=255)
       tCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCCAAGTGGTCGCCGTAGAAACGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCCTCCGGGTGATCTGCcgcgc                                                                                                                                                                                                                            >  1:386516/1‑251 (MQ=255)
         cccggggTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGt                                                                                                                                                                                                                          <  2:382732/247‑1 (MQ=255)
          ccTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTGATGCCGCAAATCAAAGGGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCcgcgcgg                                                                                                                                                                                                                          >  1:264017/1‑249 (MQ=255)
             gggTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCGGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGATGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCAGCCTGTTGATCTGCCGCGCGTCggg                                                                                                                                                                                                                      <  2:261457/251‑1 (MQ=255)
                  gcgccgTGCGACCAGCCGATTGTTATTGCCGGTCTAGGCACTCTGGCTCTGTGAATGCTCCAGTAGGACGCTCATCTCGATCGCGTATTTGGGCAAGGCGTCTGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACAGATGCCGCCAATCAAAGAGATCGCCGTAGAAGTGGAAGCCTCCGCCTGCCTGAAAGCAGCGCTTAATGCGGGTCATCCGGTTGATCTGCAGCGCGTAGGGCTAtt                                                                                                                                                                                                                 <  2:349644/250‑1 (MQ=255)
                      cgTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGGTGCCGCAAATATAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGGAAGCAGCGCTGAATGCGGGGCATCCGGGTGATCTGCCGCGCGTAGGGCTATTTGCt                                                                                                                                                                                                             >  1:292749/1‑251 (MQ=255)
                       gTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTg                                                                                                                                                                                                           <  1:27653/251‑1 (MQ=255)
                          cGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGAGATCGCCGGAGAAGCGGAAGACTCCGCCTGCCCGAAAGCAGCGCTAAATGAGGGTCATCCGGTTGATCTGGCGCGCGTAGGGCTATTTTCTGa                                                                                                                                                                                                          >  1:388085/1‑249 (MQ=255)
                            aCCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCAGAACAAAAAAATGAAGGCGCAAATCCAAGGGATCGCCGGAGAAGCGGAAAACACCGCCTGCCTGAAAGAAGCGCTGAAAGCGGGGAATCCGGTTGATCTGCCGGGCGGCGGGGTATTTGGTGAAggg                                                                                                                                                                                                      >  1:307629/1‑250 (MQ=255)
                              cATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGt                                                                                                                                                                                                    <  1:314913/251‑1 (MQ=255)
                                 ccGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAg                                                                                                                                                                                                  <  1:399074/250‑1 (MQ=255)
                                 ccGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGc                                                                                                                                                                                                 <  1:153957/251‑1 (MQ=255)
                                    agggTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCGCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGGTGTCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGGGAGCGCCGTAGAAGCGGAATACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGGTCTTCCGCGCGTGGGGCTATTTGCTGAAGGCGTAGCGGt                                                                                                                                                                                              <  2:59459/249‑1 (MQ=255)
                                      ggtggtTGCCGGGCAAGGCACGCTGTCGCTGGAACTGCTCCAGCAGGACGCCTATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCGGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCTCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCGGAAAGCAGCGCTGAATGCGGTTCAGCCGGTTGATCTGCCGCGCGGAGAGCTATTTGTTGAAGGCGTGGCGGTaa                                                                                                                                                                                            <  2:182639/251‑1 (MQ=255)
                                      ggTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGGTGATCAAACAAATGATGCCGCAAATAAAAGGGATCGCCGGAGAAGCGGAAGACTCCGCCTGTCTGTAAGCAGCGCTGAATGCGGGTCATCCGGGTGATCTGCCGCGCGTAGGGCTATTTTCTGAAGGAGTAGCGGTa                                                                                                                                                                                             >  2:373205/1‑250 (MQ=255)
                                      ggTGATTGCCGGGAAAGGCACGCTGGCGCTGGATCTGCTACAGCAGGACGCCCATCTCGACCGCGTATTTGGGCCAGTCGGCGGTGGCGGTCTGTGTGCTGGCGTGGCGGTGCTGATCAAACAACCGATGCCGCAAATCAAAGTGATCGCAGTAGAAGCGGAAGGCTCCGCGTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTaa                                                                                                                                                                                            <  1:385888/251‑1 (MQ=255)
                                        tGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGGGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTGGGAGACGGAGACTCCGCCTGCCCGGAAACAGCGCTGAATGCGGGTCATCCGGTTGGTATGCCGCGCGTAGGGCTATTTGCTTGAAGCGGAGCGGTaaa                                                                                                                                                                                           >  1:56908/1‑250 (MQ=255)
                                         gATTGCCGGGCAAGGCACGCTGGCGCGGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTAGTTGCTGAAGGCGTAGCGGTAAAAc                                                                                                                                                                                         <  2:343319/251‑1 (MQ=255)
                                              ccGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATc                                                                                                                                                                                    <  1:43353/251‑1 (MQ=255)
                                                  gCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACCTAAGGCGCAAATCAAAGGGATCGCCGGAGGAGCGGTAGACTCCCCCTGACTGAACACAGCGCGGAATTCGGGGCATCCGGGTGGTCTGCCGCGCGTCTGGGTATTTGCTGAAGGGGGAGCGGGTAAACGCCTCgggg                                                                                                                                                                                >  2:216852/1‑249 (MQ=255)
                                                    aaGGCACGCTTGCGATGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAAGCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGATGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTTAAGGCGTAGCGTTAAAACGCATCGGTGAc                                                                                                                                                                              <  1:317442/251‑1 (MQ=255)
                                                    aaGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCCAATCCAAGTGATCGCCCTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCCGCGCTGAATGCCGGGCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGGAACGGGAAAACGCATCgggggc                                                                                                                                                                              >  1:258231/1‑247 (MQ=255)
                                                         aCGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGTGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAAc                                                                                                                                                                         <  1:337901/251‑1 (MQ=255)
                                                             tggcgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGAGTGGCGGTGCTGATCAAACACCTGCTGGCGCACATCAAAGTGATCGCCGGAGAAGCGGAAGACTCCCCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATTTGCCGCGCGAAGGGCTCTTTGCTGCAGGCGTAGCGGGAAAACGGATCGGTGACGGAACCTTc                                                                                                                                                                     >  2:425503/1‑251 (MQ=255)
                                                               gcgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTTATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTAGTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAATCCTTCCg                                                                                                                                                                   <  2:270401/251‑1 (MQ=255)
                                                                       aaCTGCTCCAGCAGGCCGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGAGGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGTCTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACACATCGGTGACGAAACCTTCCGTTTATGc                                                                                                                                                            <  2:99868/251‑1 (MQ=255)
                                                                        aCTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTGATGCCGCCAATCAAAGGGATCGCCGTAGAAGCGGCAGACTCCGCCTGCCTGAAAGGAGGGCTGAATGAGGGGCATCCGGGTGATCTGCCGCGCGAAGGGGTATTTGGTGAAAGCGGAGCGGGAAAACGAAACGGTGACGAAACCTTCCGGTTATggc                                                                                                                                                           >  1:280415/1‑249 (MQ=255)
                                                                        aCTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGcc                                                                                                                                                           <  2:58709/251‑1 (MQ=255)
                                                                          tGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGCTCAAACAACCGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGGAGCGGGAAAACGCATCGGTGACGAAACCTTCCGTTTATGcc                                                                                                                                                           >  1:410946/1‑249 (MQ=255)
                                                                              ccAGCAGGCCGCCCATCTCGACCGCGTATGTGGGCCAGTCGGCGGCGGCGGTCTTGCTGCTGGCGTTGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTTATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTTCCGCGCGTAGGGCTATTTGCTGAAGTCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGt                                                                                                                                                     <  2:307629/251‑1 (MQ=255)
                                                                                cgcagGACGCCCATCTCGACCGCGTATTGGGGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTAt                                                                                                                                                   <  2:172833/250‑1 (MQ=255)
                                                                                agcagGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCCAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGGTGATCTGCCGCGCGTAGGGCGATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAggat                                                                                                                                                   >  2:215583/1‑248 (MQ=255)
                                                                                   agGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGGGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGCATAtctc                                                                                                                                                >  1:353445/1‑251 (MQ=255)
                                                                                    ggACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTTCTGATCCAACAACTGATGCCCCAAATCCAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCCCGCGTAGGGCTATTTGGTGAAGGCGTTGCGGTAAAACGCAACGGGGACGAAACCTTCCGTTTATGCCAGGAggataccg                                                                                                                                               >  1:371311/1‑244 (MQ=255)
                                                                                        gCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCCAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAAAGCGGGGCATCCGGGTTATCTGCCGCGCGGAGGGCTCTTTGCTGGAGGCGGACCGGTAAAACGCACCAGGGACGAAACCTTCCCGTTATGCCAGGAGGATCTcgtcgc                                                                                                                                           >  2:328855/1‑247 (MQ=255)
                                                                                           cATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGTCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGCcat                                                                                                                                        <  2:326982/251‑1 (MQ=255)
                                                                                                cGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGTAGTATCTCGACGAcatca                                                                                                                                      <  2:392755/248‑1 (MQ=255)
                                                                                                cGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAAAAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAAACTCCGCCTGCATGAAAGAAGCGCTGAATGCGGGTCATACGGGTGAACTGCCGGGGGTAGGGCTATTTGCTGAAGGCGGAGCGGTAAAACGCATCGGGGAAGAAACCTTCCGGTTTTGCCAGGAGTATTTTGACGAcataata                                                                                                                                    >  1:84899/1‑246 (MQ=255)
                                                                                                  aCCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCAcc                                                                                                                                 <  1:106332/251‑1 (MQ=255)
                                                                                                   ccGCGTAGTTGTGCCAGGCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCCCCg                                                                                                                                <  2:334560/251‑1 (MQ=255)
                                                                                                            tgtgCCAGGCGGCGGCGGCGGTATGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgc                                                                                                                        <  2:132850/250‑1 (MQ=255)
                                                                                                            ggggcaaGTCGGCGGCGGCGGGCTGTCTGCTGGCGTGGCGGTGCTGATCAAACAACTTATGCCGCAAATCAAAGTGAGCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGTTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgc                                                                                                                        <  1:214444/244‑1 (MQ=255)
                                                                                                               gCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgat                                                                                                                     <  1:431472/250‑1 (MQ=255)
                                                                                                                 cAGTCGGCGGCGTCGGGCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgatgcg                                                                                                                  <  1:55007/251‑1 (MQ=255)
                                                                                                                      ggcggcggcggTCTGGCTGCTGGCGTGGCGGTGCCGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCCCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTTAAGGCGTAGGGGTAAAACGCCTAGGTGACGAAACCTTCCGTTTATGCCAGGCGAATCTCGACGACCTCATCACCGTCGATAGCGATGCGATCtg                                                                                                             >  1:15725/1‑251 (MQ=255)
                                                                                                                         ggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAAAAACTGATGCCGCAAATCAAAGGGATCGCCGTAGAAGCGGGAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGGGATTTGGTGGAGGCGGAGAGGGAAAACGCATCGGGGACGAAAACTTCCGTTTATGCCAGGGGGAACTCGGCGACATCATCACCGCCGATAGCGATGCGCTCtgtgg                                                                                                          >  2:344116/1‑250 (MQ=255)
                                                                                                                               ggTCGGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCGGCAAATCAAAGTGATCGCCGGAGAATCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCAGCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAGCGCATCGTTGACGAAACCTTCCGTGTATGCCAGGAGTATCTCGACGACATCAGCACCGTCGATAGCGATTCGATCTATGCTGCGAt                                                                                                    <  2:409032/251‑1 (MQ=255)
                                                                                                                                         ctgGCGTGGCGGTGCTGATCAAAAAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAGAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACCTCATCACCGCCGATAGCGCTGCGATCTGTGCGGCGATGAAGGAttt                                                                                           >  1:104400/1‑250 (MQ=255)
                                                                                                                                             cgtggcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCg                                                                                      <  2:353445/251‑1 (MQ=255)
                                                                                                                                                   ggTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCTCCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGAtg                                                                                <  2:383722/251‑1 (MQ=255)
                                                                                                                                                   ggTGCTGATCAAACAACGGATGCCGCAAATCAAAGTGCTCGCCGTAGAAGCGGCAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGCTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGAtg                                                                                >  2:420184/1‑251 (MQ=255)
                                                                                                                                                      gCTGATCAAACAACTGATGCCGCAAATAAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCGGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTgc                                                                             <  2:239535/251‑1 (MQ=255)
                                                                                                                                                        tGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGCAGATGTgcgc                                                                           >  2:322240/1‑251 (MQ=255)
                                                                                                                                                                     gATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAAc                                                              >  2:244207/1‑251 (MQ=255)
                                                                                                                                                                      aTGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAACGCAGCGATGAATGCGGGTCATCCGGTTGATCTGCCGCGAGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGAGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAAcc                                                             >  2:225081/1‑251 (MQ=255)
                                                                                                                                                                        gCCGCAAATCAAAGTGGTCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCct                                                           <  2:367489/251‑1 (MQ=255)
                                                                                                                                                                          cGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCctct                                                         <  1:357706/251‑1 (MQ=255)
                                                                                                                                                                           gCAAATCAAAGTGATCGCCGTAGAAGCGGCAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGGGGCGGCACCCTCTg                                                        >  2:372716/1‑251 (MQ=255)
                                                                                                                                                                            cAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTgg                                                       <  2:264017/251‑1 (MQ=255)
                                                                                                                                                                            cAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTgg                                                       >  1:356979/1‑251 (MQ=255)
                                                                                                                                                                             aaaTCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGATCCCTCTGgc                                                      <  2:386516/251‑1 (MQ=255)
                                                                                                                                                                                     gggATCGCCGTAGAAGCGGAAGACGCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCGCTGGCgcgctggc                                              <  2:242856/249‑1 (MQ=255)
                                                                                                                                                                                     gTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggc                                              <  1:373073/251‑1 (MQ=255)
                                                                                                                                                                                      tGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggcg                                             >  1:168404/1‑251 (MQ=255)
                                                                                                                                                                                        aTCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggcgct                                           <  1:91461/251‑1 (MQ=255)
                                                                                                                                                                                         gcGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATGTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggcgctg                                          <  1:412225/250‑1 (MQ=255)
                                                                                                                                                                                          cGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTGCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggcgctg                                          <  1:425503/250‑1 (MQ=255)
                                                                                                                                                                                           gCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCgctggcgctggc                                        <  2:124105/251‑1 (MQ=255)
                                                                                                                                                                                           gCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCgctggcgctggc                                        <  1:325051/251‑1 (MQ=255)
                                                                                                                                                                                              gTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCggg                                     <  1:344116/251‑1 (MQ=255)
                                                                                                                                                                                                aGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGa                                    <  2:366253/250‑1 (MQ=255)
                                                                                                                                                                                                 gAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAAt                                  <  1:215583/251‑1 (MQ=255)
                                                                                                                                                                                                 gAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAAt                                  >  1:130180/1‑251 (MQ=255)
                                                                                                                                                                                                 gAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTGGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTGATGCCAGGAGTATCTCGACGACAGCATCGCCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATTTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAAt                                  <  1:180792/251‑1 (MQ=255)
                                                                                                                                                                                                  aaGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGGCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCTGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATg                                 <  1:194204/251‑1 (MQ=255)
                                                                                                                                                                                                  aaGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATg                                 <  1:167793/251‑1 (MQ=255)
                                                                                                                                                                                                      cGGAAGACCGCGCATGCCTGAAGGCAGCGCAGAATGCGGGTCATCCGGTTGCTCTGGCGCGCGTAGGGCTAGTTGCTGAAGGCGTAGCGGTAAAACGCAGCGGTGACTAAGCCTGCCGTTTATGCCTGGAGTATCTCGACGACATCATCGCCGTCGATAGCGAGGCGTTCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGGGGTGGCGGAACCCCCTGGCGCGCTGGCGCTGGCGGGAATGaaa                              <  1:169834/251‑1 (MQ=255)
                                                                                                                                                                                                      cGGAAGAATCTGCCTGCCTGAATGCAGGGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGAGTAGAGGTAAATCGCATCGGTGACGACACCTTCCGATTATGCCAGGAGTATCTCGACGACATCAACACCGTCGATAGCGATGCGATCTGTGCGGCGATGAGGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGAGCTGGCGGGAATGaa                               >  2:181838/1‑250 (MQ=255)
                                                                                                                                                                                                       ggAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaaaa                             <  1:395687/251‑1 (MQ=255)
                                                                                                                                                                                                        gAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaaaaa                            <  1:134641/251‑1 (MQ=255)
                                                                                                                                                                                                         aaGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAACCCTTCCGTTTATGCCAGGAGAATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGGGCTGCCGCTGGCGGGAATGaaaaa                            >  2:143722/1‑250 (MQ=255)
                                                                                                                                                                                                           gacgccGCCTGCCTGAAAGCAGCGCGGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTTCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCTGGAATGAAAAAata                         <  2:292749/247‑1 (MQ=255)
                                                                                                                                                                                                                      actgAAAGCAGCGCTGAATGCGGGTCATACGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAGGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCCGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGACGATGCAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCTGGAATGAACAAATATATCGCCCTGc              >  2:263862/2‑251 (MQ=255)
                                                                                                                                                                                                                         gAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATATCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGTTGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAACAATATATCGCCCTgaaa            >  1:433174/1‑247 (MQ=255)
                                                                                                                                                                                                                           aaGCAGTGCTGAATGCGGGGCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATAGATCGCCCTGCacaac         <  2:56908/251‑1 (MQ=255)
                                                                                                                                                                                                                           aaGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAACACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTCTCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCacaac         >  1:427501/1‑251 (MQ=255)
                                                                                                                                                                                                                             gcagcGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACAt       <  2:241916/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gcTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGCAAAAATAAATCGCCCTGCACCACATTcgcg  >  2:300067/1‑251 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
CAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCG  >  NC_000913/3955775‑3956249

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