Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,396,772 T→C F236F (TTT→TTC)  amiB → N‑acetylmuramoyl‑L‑alanine amidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,396,7720TC98.8% 259.9 / ‑4.4 82F236F (TTT→TTC) amiBN‑acetylmuramoyl‑L‑alanine amidase B
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (38/41);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

GAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGA  >  NC_000913/4396527‑4397016
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
gAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACTCTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCTATGTTCAAAgg                                                                                                                                                                                                                                                 <  1:593473/251‑1 (MQ=255)
  aCCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCg                                                                                                                                                                                                                                               <  2:465303/251‑1 (MQ=255)
     ggCGCGCAAGCCGGTTAAAACGGAAAGTCACCGCACGACGGGTGTTATCAGCAGTAATCCGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCAGCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACGTTGCTCAATGACGATCCGATGTTCAAAGGCGttt                                                                                                                                                                                                                                            <  2:250263/251‑1 (MQ=255)
            aaTCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCg                                                                                                                                                                                                                                     >  1:72649/1‑251 (MQ=255)
             aTCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGt                                                                                                                                                                                                                                    <  1:291698/251‑1 (MQ=255)
               ccGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTg                                                                                                                                                                                                                                   <  2:145837/250‑1 (MQ=255)
                    tAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACggg                                                                                                                                                                                                                              <  1:601580/250‑1 (MQ=255)
                       aaCGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGGGGTACGCGGGAGAAAAATGTCACCATCGCCCTCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGAtt                                                                                                                                                                                                                          >  1:87296/1‑251 (MQ=255)
                                  aCCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCg                                                                                                                                                                                                               <  1:435094/251‑1 (MQ=255)
                                          ccGGGTGGTACCAGAAGTAATTCGGTAACGCGACCGGCAGGTCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTAGTGAGGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGCGATggg                                                                                                                                                                                                       <  2:355614/250‑1 (MQ=255)
                                                   aTCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTg                                                                                                                                                                                               >  1:189451/1‑250 (MQ=255)
                                                   aTCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTg                                                                                                                                                                                               <  1:304167/250‑1 (MQ=255)
                                                         agTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGc                                                                                                                                                                                        <  1:90153/251‑1 (MQ=255)
                                                           tAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAAAAATGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGGTGTTCCAAGGCGTTTTAACCCGTGACGGGGATTACTTTATTTCGGTGATGGGGAGCAGTGATGTGGCAc                                                                                                                                                                                      >  1:186534/1‑251 (MQ=255)
                                                              aaCGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTa                                                                                                                                                                                   <  2:72649/250‑1 (MQ=255)
                                                               aCGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGGGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGGGATGTGGCACGTa                                                                                                                                                                                   >  1:166250/1‑250 (MQ=255)
                                                                    aaCGCGTCCGGCAGCGCGCGCGACGGATAACACTTGGGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTAAGGCCCCTTGTGCTATCGGCACCGGTGGGACGCTGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGAGCCGATGTTGAAAGGCGTTTTATCGCGTGACGGGGATTACGTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCaa                                                                                                                                                                              <  2:361025/250‑1 (MQ=255)
                                                                           ccGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGGGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGGCCGTAAGCAAAACGCCa                                                                                                                                                                       >  1:580168/1‑250 (MQ=255)
                                                                            cGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCCGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCCTCGCCATCGCGCGTAAATTGCGTACTTTGCTCCATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCCGGGATGTGGCACGTAAGCAAAACGCCAAt                                                                                                                                                                     >  1:238376/1‑251 (MQ=255)
                                                                              gCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAAttt                                                                                                                                                                   <  1:42771/251‑1 (MQ=255)
                                                                              gCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAAttt                                                                                                                                                                   >  1:74972/1‑251 (MQ=255)
                                                                                  cgcgcgcgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTc                                                                                                                                                               <  1:697798/251‑1 (MQ=255)
                                                                                           cGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGAtt                                                                                                                                                      <  1:520816/251‑1 (MQ=255)
                                                                                            ggCTAACACTGGCGCTAAAATTATCATCGCTATTGATTCCGGGCACGGCGGTCATGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAAGGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACggaa                                                                                                                                                                                  <  2:467686/222‑4 (MQ=255)
                                                                                            ggCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCCGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGGAAAATGTCACCCTCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGGTCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGGGCCGGGATGTGGCACGTaa                                                                                                                                                                                  >  1:467686/1‑222 (MQ=255)
                                                                                                    cTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTCAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGGGGAGGGGTGAGGCAAGTAAGCCAAAAGCACATTTCCTCGGGGCGATTCACGCTgcg                                                                                                                                             >  2:474014/1‑249 (MQ=255)
                                                                                                     tGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGGGGGACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGGAAATTGCGTACTTTGCTCAATGACGATCCGCTGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGGGGTGGGGCGCCGTGGTGTGGCACGGAAGCAAAACGCCAATTTCCTCGTGTCGGATCACGcggg                                                                                                                                              >  2:254954/1‑246 (MQ=255)
                                                                                                      ggCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGAACCTTGCGCTATCGGCCCCCGGGGTACGCGGGAGGAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTAATTTGCTCAATGAAGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGCGTCGATTCCAGCTGATGc                                                                                                                                           >  1:146468/1‑251 (MQ=255)
                                                                                                        cGATAAAATTATCATCGCTATTGATGCCGGCCACGGCGGTCGAGACCCTGGAGCTTTCGGACCCGGTGGTCCGCGGGGGAAAAAAGTCACCATCGCCATCCCGCGAAAATTGCGGACTGTGCTCCATACCGATCCGATGTTCAAAGGCGTTTTCACCCGGTACGGGGTGTACTTTTTATTGGGGGTGGGGCGTGGGGGTGTGGACCGGCAGCAAAACGCCAAATTCCTCGTGTCGATTcccggttgtgccc                                                                                                                                         >  2:338663/1‑239 (MQ=255)
                                                                                                        cGATAAAATTATAATCGCTATTGATGCCGGACACGGAGGTCAGGAACCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGGAAAATGTCACCATCGCCATCGAGCGTAAATTGCGTACTTTGCTCAATGACGATACGATGTTCAAAGGCGTTTTAAACCGTGACGGGGAATACTTTATCTCGGGGATGGGGCGCAGTGATGGGGCACGTAAGAAAAAAGCCCATTTCCTCGTGTCGCTTCCCGATGGTGCCg                                                                                                                                         >  2:471509/1‑251 (MQ=255)
                                                                                                         gATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCCATCGGCCCCGGTGGGACGCGGGAGAAAAATGTCACCATCGCCATCTCGCGTAAATTGCGTACTTTTCTCAATGACGCTCCGATGTTCAAAGGCATTTTAACCCGTGACGGGGATTACTTTATCTCGGTAATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGc                                                                                                                                        >  1:473314/1‑251 (MQ=255)
                                                                                                          aTAAAATTAGCAGCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTGTTAACCCGTGACGGGGATTTCTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGAGGCCGCa                                                                                                                                       <  2:460851/251‑1 (MQ=255)
                                                                                                                  atcatcGCTATTGATGCCGGACACGGCGGTCAGGCACCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAAAGTCACCATCGCCCTCGCGCGGAAATTGCGTACTTTGCTCAATGACGATCCGATTTTCAAAGGCGTTTTAACCCGTGGCGGGGATTACTTTATCTCGGTGCTGGGGCGCAGGGGTGTGGCACGTAAGCAAAAAGCCAATTTCCTCGTGTCGATTCACGCGGATGCCGCACCGAAccc                                                                                                                               >  2:297713/1‑250 (MQ=255)
                                                                                                                  atcatcGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGTGCTCTCGGTCCCGGGGGTACGCGGGGGAAAAATGTCACCCTCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGGTCCGAGGTTCAAAGGCGGTTTAAACCGGGACGGGGATTACTTTATCTCGGTGCTGGGGGGCAGTGAGGGTGCACGTAAGCAAAAAGCCCATTTTCTCGTGTGGATTCACGCTGTTGCCGCCCCgccacg                                                                                                                               >  2:16423/1‑246 (MQ=255)
                                                                                                                   tcatcGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGGTGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGc                                                                                                                              >  1:206790/1‑251 (MQ=255)
                                                                                                                         ctcttggggtcGGACACGGCGGTCAGGTACCTGGCGCCTTCGGCCCCGGTGGTACGCGGGATAAAAATGTCACTAGCGCCATCGCGCTTAAAGTTCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGATGATGCCGCAATGAACCGCAGt                                                                                                                           <  2:132267/238‑1 (MQ=255)
                                                                                                                         cTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCg                                                                                                                        >  1:644644/1‑251 (MQ=255)
                                                                                                                          tATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGGGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCgg                                                                                                                       >  1:175103/1‑250 (MQ=255)
                                                                                                                                 gCCGGACACGGCGGTCAGGGCCCCGTCGCCAGCGGCCCCGGTCGTACGCGGGAGAAACAAGTCACCACCGACATCGCGCGTAAATACCGTACTTTGCGCAAGGCCGATCCGATGGTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATTTGGCACGTAAGCAAAACGCCAATTTCCGCGTGTCGATGCACGCTGATGCCGCACCGAACCGCAGTGCGACTGgcg                                                                                                                 <  2:367950/250‑1 (MQ=255)
                                                                                                                                    ggACGCGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGGCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTc                                                                                                             <  2:580168/251‑1 (MQ=255)
                                                                                                                                    ggACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAAGGCGACTGGCGCTTc                                                                                                             >  1:463086/1‑251 (MQ=255)
                                                                                                                                       cacGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGt                                                                                                          <  1:621432/251‑1 (MQ=255)
                                                                                                                                         cggcggTCAGGACACTGGCGCTATCGGACCCGGTGGGACGCGGGGGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCCAAGGCGTTTTAACCCGTGTCGGGGATTACCTTATCTTGGTGATGGGGCGCAGTGATGTGGCAAGTAAGCAAAACGCCCATTTTCTCGTGTCGATTCACGCTGATTCCGCACCGAACCGCAGGGCGGCTGGCGCTTCCGtct                                                                                                        >  2:19234/1‑249 (MQ=255)
                                                                                                                                          ggcggTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATg                                                                                                       >  1:702242/1‑251 (MQ=255)
                                                                                                                                              gTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTTCTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTg                                                                                                   <  1:471822/251‑1 (MQ=255)
                                                                                                                                                     cccTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGCACCGCAGTGCGACTGGCGCTTCCGTATGGGGGctctct                                                                                             >  2:371457/1‑250 (MQ=255)
                                                                                                                                                        tGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGCCTGGCGCTTCCGTATGGGTGCTCTCTAAcc                                                                                         >  1:148546/1‑251 (MQ=255)
                                                                                                                                                         ggCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCg                                                                                        >  1:160774/1‑251 (MQ=255)
                                                                                                                                                            gcTATCGGCCCCGGTGGTACGCGGGAGAAACATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCCGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGGTTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGATTCCGTTTGGGTGCTCTCTAACCGTcg                                                                                     >  2:474697/1‑251 (MQ=255)
                                                                                                                                                              tATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTcgcg                                                                                   >  1:233585/1‑251 (MQ=255)
                                                                                                                                                                 cGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCa                                                                                <  1:149400/251‑1 (MQ=255)
                                                                                                                                                                    ccccGGTGGTACGAGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGCCTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGccccc                                                                              >  2:40411/1‑247 (MQ=255)
                                                                                                                                                                       cGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTg                                                                          <  1:509999/251‑1 (MQ=255)
                                                                                                                                                                            gtACGCGGGAGAAAAAGGTAGCCAGCGCCATCGCGCGAAAATTGCGTCCTTTGCACAATGACGTTACGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGATTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAAt                                                                      <  2:240135/250‑1 (MQ=255)
                                                                                                                                                                                    gataaaaaTGTCCCCATCGCCATCGCGCGTAAATGGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCGCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAgctg                                                             <  2:74972/248‑1 (MQ=255)
                                                                                                                                                                                    gagaAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGATGATCCGATGTTCAAAGGCGGTTTAACCCGTGACGGGGATTACTTTATCGCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAgctg                                                             <  1:231373/251‑1 (MQ=255)
                                                                                                                                                                                      aaaaaaGGTCACCATCTCCATCGGGCGTAAATTGCGTTCTGTGCTCAATGCCGATCCTATGTTCAATGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTGCCTCGTGTCGAGTCACGCTGATGCCGCGCCGAACCTCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAgctggc                                                           <  2:410942/250‑1 (MQ=255)
                                                                                                                                                                                          aatgaccccatcgccaccgCGCGTAAAGTGAGTACCTTGCTCAATGGCGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATGACTTTATCGCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCGCCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCGCTCTAACCGTCGCGCCAACAGTGAAATGGCCAgctggctgta                                                       <  2:287351/244‑3 (MQ=255)
                                                                                                                                                                                           aTGTCACCATCGCCATCGCGCGTAACTTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGctggctgg                                                        >  1:201013/1‑249 (MQ=255)
                                                                                                                                                                                           aTGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGag                                                      <  1:297713/251‑1 (MQ=255)
                                                                                                                                                                                                   atcgccatcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGACCCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGCGCAGCACg                                               >  2:337524/1‑250 (MQ=255)
                                                                                                                                                                                                   atcgccatcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTCACCGTCGCGCCCACAGTGGAATGGCCAGCTGGCTGGa                                                       >  2:281866/1‑242 (MQ=255)
                                                                                                                                                                                                   atcgccatcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGa                                                       <  1:281866/242‑1 (MQ=255)
                                                                                                                                                                                                       ccagcgcgcgTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGaaa                                          <  2:189451/247‑1 (MQ=255)
                                                                                                                                                                                                         atcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCACTTTCCTCGTGTCGATTCACGCTGATGCCGCACCGACCCGCAGTGCGACTGGCGCTTCCGTATGGGGGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGCGAAAcc                                        >  2:615725/1‑250 (MQ=255)
                                                                                                                                                                                                            gcgcgTAACTAGCGTACTTTGCTCAATGACGATCCGATGATCAAGGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGCTGTGGCACGTAAGCAAAACGACAATTTCCTCGTGTCGATTCACGATGATGCCGCACCGAACCGCAGTGCGACAGGCGCTTCCGTATGGGTGCTATCTAACCGTCGCGCCACCAGTGAAATGGCCAGcaggctggggacaccagagcaacaatc                                     >  2:229280/1‑232 (MQ=255)
                                                                                                                                                                                                                     ttGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAgctgct                            >  1:109969/1‑251 (MQ=255)
                                                                                                                                                                                                                      tGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGAGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGCCGCTTAACGCTGCTGCCGCACCGAACCGCCGTGCGACTGGCGCTTCAGTATGGGTGCTCTCTAACCGTCGCGCCAACAGGGAAATGGACAGCAGGCTGGGACAGCACGGAAAAAAGTAGGAgctgt                             >  2:519401/1‑248 (MQ=255)
                                                                                                                                                                                                                            cTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCtgggtggg                     >  1:72075/1‑251 (MQ=255)
                                                                                                                                                                                                                              ttGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTgggg                    <  1:471509/250‑1 (MQ=255)
                                                                                                                                                                                                                                ttttaaTGACTATCCGATGTTATAAGGCTGTTGAACCCGTGGCGGGGATTACGTGATCGCTGTTATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCGCGTGTCGATTCACGCGGATGCCGCCCCGACCCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAAGGGCCAGCTGGCTGGAGCAGAACGAGAAACAGTCGGAGCTGCTGGGTGGGGCg                  <  2:651392/246‑1 (MQ=255)
                                                                                                                                                                                                                                ccTCAATGCCGATCCGATGTTCAAAGGCGTTGTAACCCGGGACGGGGAGTACTTTATCTCGGTGAGTGGGCTCAGTGAGGTGGCACGTAAGCAAAACGCCAATGTCCGCGTGGCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACCGGCGCTTCCGTTGGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCCGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCGGGGTGGGGCg                  <  2:555967/249‑1 (MQ=255)
                                                                                                                                                                                                                                          gATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCt       <  1:474014/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCt       <  1:571898/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aTCCGATGTTCAAAGGGGTTTTAACCCGTGACGGGGATTACTTTAGCTCGGTGCTGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTg      <  2:238376/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aTCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGaa                                           >  2:430322/1‑214 (MQ=255)
                                                                                                                                                                                                                                           aTCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGaa                                           <  1:430322/214‑1 (MQ=255)
                                                                                                                                                                                                                                           aTCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTg      <  1:254954/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aTCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTg      <  1:188838/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ccGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGc    <  1:19234/251‑1 (MQ=255)
                                                                                                                                                                                                                                                aTGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCg   <  1:634969/249‑1 (MQ=255)
                                                                                                                                                                                                                                                aTGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGGGGGGCGGGGGATGGGCTGGCgg  >  1:234123/1‑249 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
GAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGA  >  NC_000913/4396527‑4397016

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