Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,523 A→G A82A (GCT→GCC)  pyrE ← orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,5230AG98.8% 264.1 / NA 82A82A (GCT→GCC) pyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (13/68);  total (14/68)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAG‑CGGTTGTGGTGGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATG  >  NC_000913/3815282‑3815769
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                           
cgccgcGCCCCCGTTCCTGACGAGCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGt                                                                                                                                                                                                                                                       <  2:566388/247‑1 (MQ=255)
   cgcgCCCCCGTTCCTGACGATCAAGCGAAATCAATCCGCCAGCAAGCTGCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAgtggtc                                                                                                                                                                                                                                                <  2:515153/251‑2 (MQ=255)
   cgcgCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAgtggtc                                                                                                                                                                                                                                                <  1:52170/251‑2 (MQ=255)
    gcgcCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAgtggtcg                                                                                                                                                                                                                                               <  2:696029/251‑3 (MQ=255)
     cgctcccgTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTGGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAgtggtcgc                                                                                                                                                                                                                                              <  1:729566/245‑4 (MQ=255)
       cccccgTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGGTTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCaa                                                                                                                                                                                                                                            <  1:351316/247‑1 (MQ=255)
       cccccgTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCCTGATGCCCCGCCAGCGCCAcgg‑cgggagtgggcccc                                                                                                                                                                                                                                             >  1:515487/5‑239 (MQ=255)
        ccccgTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAAt                                                                                                                                                                                                                                           <  2:240874/248‑1 (MQ=255)
              cccTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACAGCATCTACCAGCATTTCGCGTCCTTGTAACGCGCTACCAACCAGATTTCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑At                                                                                                                                                                                                                                    <  2:175822/250‑1 (MQ=255)
               ccTGACGATCAAGCGAAATCAATCCGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTGGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCTCGG‑CGGTAGTGGTCGCAATCGGG‑ATc                                                                                                                                                                                                                                   <  1:671277/251‑1 (MQ=255)
                 ggACGATCAAGCGAAATCAATACGCCAGCAANCGTCGCCCCATTCGCCTGAATAATCGCCATCGATTCGCGAATCGCCGTTCCGGCGGTGCTCACATCATTTTCCAGCATTACGCGACCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATccc                                                                                                                                                                                                                                 <  2:368143/250‑1 (MQ=255)
                    cGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTGGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCttg                                                                                                                                                                                                                              <  2:487258/251‑2 (MQ=255)
                      atcaaGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGt                                                                                                                                                                                                                            <  1:345796/247‑1 (MQ=255)
                            cGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAg                                                                                                                                                                                                                      <  1:221125/251‑1 (MQ=255)
                            cGAAATCAAAACGCCAGCAAGCGGCGCGCCATTAGCCTGAATAATCTCCAGCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTACCGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGAGGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAg                                                                                                                                                                                                                      <  2:475513/251‑1 (MQ=255)
                              aaaTCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAggg                                                                                                                                                                                                                    <  1:267607/251‑1 (MQ=255)
                               aaacaatacGCCAGCAAGCGTCGCGCCATTANCCTGAATAATCTCCATCGATGCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTCCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCGCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGc                                                                                                                                                                                                                   <  2:350629/244‑1 (MQ=255)
                                aTCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGcc                                                                                                                                                                                                                  <  1:145037/251‑1 (MQ=255)
                                         ccAGCAAGCGTCGCGCCAGTAGCCTTAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTCCCAGCATGACGCGTCCTTGTAACGCGCCACCAACCAGATTGCCGCCTTCACCGTTGTCTTTTGCCTCTTTGCGGGTAAAGCAGCACGGCAGTTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAAcagcag                                                                                                                                                                                                         <  2:244346/251‑1 (MQ=255)
                                             cAAGCGTCGCGCCATTAGCCTTAATAATCTCCATCGATTCGCGAAGCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCg                                                                                                                                                                                                     <  2:577760/251‑1 (MQ=255)
                                               aGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGTCATGATTCTCCGCCCGCGCCACGG‑CGGGAGTGGTCGCAAACGGG‑GTCCCCTTGTAAGCAGGGCCaaaaagacgaacga                                                                                                                                                                                                    >  1:84872/1‑239 (MQ=255)
                                                     gcgcCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCGTTTGCCTCTTTGGGGTTAAAGCAGTACGTCAGGTCAAGGGCATGATGCTCCGCCAGCGCGTCGG‑CGGTAGTTGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAAt                                                                                                                                                                                             <  1:435270/251‑1 (MQ=255)
                                                       gcCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCAGTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGc                                                                                                                                                                                           <  1:387123/251‑1 (MQ=255)
                                                        ccATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCTTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGcc                                                                                                                                                                                          <  1:708632/251‑1 (MQ=255)
                                                             agCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCCTGATGCTCCGCCCGCGCC‑CAGGCGGGTGTGGTCGCAAACGGGGATCC‑TTTGTAAGCCGGGGCAAAcacccggtcggaatcaatgcccaaa                                                                                                                                                                                      >  1:50363/2‑227 (MQ=255)
                                                                    ataatCTCCACCGATGCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCCTTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCGGCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGTCCAATCAGCAGATCGAACTCAATGCCAGAATCCTCCaa                                                                                                                                                                              <  2:212675/251‑1 (MQ=255)
                                                                           ccATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCt                                                                                                                                                                          <  1:517035/248‑1 (MQ=255)
                                                                               cgattCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGt                                                                                                                                                                   <  1:516382/247‑1 (MQ=255)
                                                                                           cGCCTTTCCGGCGGGGATCACAGCATCTACCAGCATTACGCGTCCTTTTAACGCGCTACCAACCAGATTGCCGCCTTCTCCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCCGCGCCGCGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCt                                                                                                                                                        <  2:256392/250‑1 (MQ=255)
                                                                                           cGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTa                                                                                                                                                       <  1:172190/251‑1 (MQ=255)
                                                                                            gCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTaa                                                                                                                                                      <  1:569429/251‑1 (MQ=255)
                                                                                                  ccGGCTGTGATCCCATCATCTACCAGCATTACGCGGCCTTGTAACGCGCGACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCTCCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGc                                                                                                                                               <  2:654657/251‑1 (MQ=255)
                                                                                                  ccGGCGGTGATCACATCATCTACCAGCATTANGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGc                                                                                                                                               <  2:355911/251‑1 (MQ=255)
                                                                                                        gTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGGAGTGGTCGCCATCCGG‑ATCCCTTTGTAAGCCGGGCCCAACAGCCGATCGAACCCCATGCCAGAAACCACCAATGCTTCAGCGTAAAAACGGCCCAAC‑AGTGcccaggc                                                                                                                                         >  1:377878/1‑246 (MQ=255)
                                                                                                          gATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGcc                                                                                                                                                                                          <  1:151135/201‑1 (MQ=255)
                                                                                                             cccatcatCTAACAGCATTACGCGTCCTTGTAACGCGCTACTAACCAGAGTGCCGCCTTCACCGTGGTCTTTCGCCGCTTGGCGGTTAAAGAAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCGCCAATGCTTCAGCGTAAAAACGGCCTAAC‑ATTGCC‑AGGTCGCGcc                                                                                                                                   <  2:331598/249‑1 (MQ=255)
                                                                                                             cacaTCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCgcgc                                                                                                                                    <  1:461766/250‑1 (MQ=255)
                                                                                                              ccatcatCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGccc                                                                                                                                  <  2:536933/250‑1 (MQ=255)
                                                                                                                          aGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCTTGATGCTCCGCCAGCGCCACGG‑CGGTAGGGGGCGCCATCGGG‑GTCCCTTTGTTAGCAAGGGCAAACCGCAGATCGCACTCCATGCCAGAATCCCCCAATGCTTCCGCGTAAAAACCGCCTAACCAG‑GCC‑AAGTCGCGCCCCGGATTAAAc                                                                                                                        >  1:706601/1‑249 (MQ=255)
                                                                                                                           gCATTTCGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGc                                                                                                                     <  2:84872/251‑1 (MQ=255)
                                                                                                                               tACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCg                                                                                                                  <  1:126447/250‑1 (MQ=255)
                                                                                                                                   cgTCCTTGTACCGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCAGGATGCGCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGtt                                                                                                             <  2:515487/251‑1 (MQ=255)
                                                                                                                                   cgTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCCGCGCCACGG‑CGGGAGTGGTCGCCATCCGG‑AACCCTTTGTAAGCCGGGGCCAACAGCAGATCGAACTCCATGCCAGAATCCACCAAAGCTTCCGCGGAAAAAACGCCTAAA‑AATGCC‑CGGTAGCGCCCGGCATTAAACAGCCCGGCGtt                                                                                                             >  2:66617/1‑251 (MQ=255)
                                                                                                                                   cgTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGtt                                                                                                             <  2:3697/251‑1 (MQ=255)
                                                                                                                                    ggccTTGTAACGCGCTACCAACCAGATTGCCNCCTTCACCGGGGTCTTTTGCCTCTGTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCAACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGATCTCAATGCCAGAATCCACCAATGCGTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCTGTATTAAACAGCCCGGCGTTg                                                                                                            <  2:342186/249‑1 (MQ=255)
                                                                                                                                      cctcgTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTgaa                                                                                                          <  1:683579/247‑1 (MQ=255)
                                                                                                                                                 gcTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGc                                                                                               <  1:348279/251‑1 (MQ=255)
                                                                                                                                                       aaccaGATGGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTgc                                                                                         <  2:307266/251‑1 (MQ=255)
                                                                                                                                                          caGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGGGGTCGCAATCGGG‑ATCCCTTTTTAAGCCGGGCCAAACCGCAGAACGAACTCCATGCCAGAATCCACCCATGCTTCAGCGTAAAAACGGCCCAAC‑AATGCC‑CGGGCGCGCCCGGTATTAAACAGCCCCGCGGTGAAGAAATAAGGGCTTTTgcggc                                                                                      >  1:654874/1‑249 (MQ=255)
                                                                                                                                                           aGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGccc                                                                                     <  2:763100/251‑1 (MQ=255)
                                                                                                                                                                gccgccTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑tt                                                                               <  1:409945/251‑1 (MQ=255)
                                                                                                                                                                gccgccTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCCAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAACCAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCAGGA‑tt                                                                               <  1:449889/251‑1 (MQ=255)
                                                                                                                                                                  cgccTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGGGGGCGCCATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCAACAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGGCGCGCCCGGGATTAAAAAGCCCGGCGTTGGAGAAAAAAGGGCTTTTGCGCCCGGG‑tttt                                                                             >  1:532983/1‑250 (MQ=255)
                                                                                                                                                                       tCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGc                                                                          <  1:12760/249‑1 (MQ=255)
                                                                                                                                                                         aCCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGAAGTGGTCGCCATACGG‑ATCCCTTTGTAAGAAGGGGCAAACAGCAGAACGAACTCAAAGCCAGAAACCACCCATGCTTCAGCGTAAAAACGGGCTAAC‑AATGCCCAG‑GCGCGCCCGGTATTAAACAGACCGACGTGTAAGAAATAGGGGCTTTTTCGCCCGGAGTTTCacgta                                                                        >  1:555646/1‑246 (MQ=255)
                                                                                                                                                                          ccGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTaaa                                                                     <  1:254614/251‑1 (MQ=255)
                                                                                                                                                                            gTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCATAATCCACCAATTCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACt                                                                   <  1:287802/251‑1 (MQ=255)
                                                                                                                                                                            gTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGTCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACt                                                                   <  1:111336/251‑1 (MQ=255)
                                                                                                                                                                              ggTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTc                                                                  <  1:380240/250‑1 (MQ=255)
                                                                                                                                                                              ggTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTc                                                                  <  1:496505/250‑1 (MQ=255)
                                                                                                                                                                              ggTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCg                                                                 <  1:326281/251‑1 (MQ=255)
                                                                                                                                                                              ggTCTTTTGCCGCTTTGCGGTTAAAGCAGGNCGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTc                                                                  <  2:377878/250‑1 (MQ=255)
                                                                                                                                                                                 cctgtgactttttGCGGTTAACGCTGTACGNGAGGTTGAGGTCATGATGCTCCGCCAGCGCTACGG‑AGGTAGTCGTCGCAATCGGG‑ATCCCTTGGTAAGCTGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCGCCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTATGCAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTAAGCGTAAAATCGcc                                                               <  2:343210/240‑1 (MQ=255)
                                                                                                                                                                                  ttttGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCa                                                              <  2:189182/250‑1 (MQ=255)
                                                                                                                                                                                   tttGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCa                                                              <  1:713375/249‑1 (MQ=255)
                                                                                                                                                                                       cctcttTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAc                                                                                                                                                     <  2:70603/159‑1 (MQ=255)
                                                                                                                                                                                            ttGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTaa                                                      <  1:430913/248‑1 (MQ=255)
                                                                                                                                                                                                ggTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCtgc                                               <  1:188804/251‑1 (MQ=255)
                                                                                                                                                                                                 gTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCCCGG‑CGGTAGTGGTCGCTATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGACCTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCtgct                                              <  1:383456/251‑1 (MQ=255)
                                                                                                                                                                                                     aaGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCtgcttgct                                          <  1:665828/251‑1 (MQ=255)
                                                                                                                                                                                                              ggCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGAACGAACTCAATGCCAGAATCCAACAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCaa                                  >  1:273190/1‑250 (MQ=255)
                                                                                                                                                                                                               gCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAt                                <  1:535397/251‑1 (MQ=255)
                                                                                                                                                                                                               gCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAt                                <  1:59480/251‑1 (MQ=255)
                                                                                                                                                                                                                 aggtcaaggtcaTGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTc                              <  1:368407/251‑1 (MQ=255)
                                                                                                                                                                                                                    tcaaggtcATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTAAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAAATCCCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCaa                            >  1:719131/4‑250 (MQ=255)
                                                                                                                                                                                                                         gtcatGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCGTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACt                      <  2:239696/247‑1 (MQ=255)
                                                                                                                                                                                                                                 gCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTg               <  2:286183/250‑1 (MQ=255)
                                                                                                                                                                                                                                 gCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTg               <  2:715506/250‑1 (MQ=255)
                                                                                                                                                                                                                                     cGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATg          >  1:602887/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cagcgCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACCGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCCCTgcatttggt       >  2:400353/5‑243 (MQ=255)
                                                                                                                                                                                                                                         agcgCCACGG‑CGGTAGTGGGCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCAACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGCGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGGTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAAATGGAGCTGATATGGttt      >  2:507852/4‑251 (MQ=255)
                                                                                                                                                                                                                                          gcgCCACGG‑CGGGTGTGGGCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGGCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAGTCAATAAACTGGCGCTGATATGGTTTc     <  2:654874/249‑1 (MQ=255)
                                                                                                                                                                                                                                             ccACGG‑CGGTAGGGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCCAACAGCAGATCGGACTCAATGCCAGAATCCACCCATGCTTCAGCGTAAAAACGGCCTAAC‑AGGGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAAATTTAACACCTGCTTGCTAAGCGCCAATTCAATAAAATGGCGCTGATATgggtttatg  >  2:145440/1‑244 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                           
CGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAG‑CGGTTGTGGTGGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATG  >  NC_000913/3815282‑3815769

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