Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,820 A→G V303V (GTA→GTG radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,8200AG98.8% 271.6 / ‑4.7 84V303V (GTA→GTGradADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (36/47);  total (36/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

TAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCA  >  NC_000913/4626573‑4627061
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
tAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATACAAACCTCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGGCACCCTTCGGCCATTTTCTTAAGTCGGGGAGAGGAAGTGGCCTTCGGAAGCTCCGTGATGGTGGTGTgg                                                                                                                                                                                                                                                  >  2:519503/1‑251 (MQ=255)
       aGATGGCTCGTTGGCGGGCCCTAAAGTGCGGGAACACGGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAATCCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTGCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAATCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGggaaggaa                                                                                                                                                                                                                                           <  2:492594/251‑1 (MQ=255)
        gatgtctCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAAc                                                                                                                                                                                                                                          <  2:453332/245‑1 (MQ=255)
             ctCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGGCTCCCGTTTTCGCACTTTACGCAGCAATAACAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGTTGGTGGTGTGGGAAGGAACGCGTc                                                                                                                                                                                                                                     >  1:288140/2‑251 (MQ=255)
               cGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTcc                                                                                                                                                                                                                                    <  1:560386/250‑1 (MQ=255)
                gttggcgggcCCTAAAGTGCTGGAACCCTGTATCGACGGTTCGTGGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTGCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTccgt                                                                                                                                                                                                                                  <  2:286672/249‑3 (MQ=255)
                   ggcgggcCCTAAAGAGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTccgttgc                                                                                                                                                                                                                               >  1:539755/5‑246 (MQ=255)
                     cgggcCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTccgttgctg                                                                                                                                                                                                                             >  1:386556/3‑244 (MQ=255)
                         ccctaaaGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCtggtgg                                                                                                                                                                                                                         <  1:463039/247‑1 (MQ=255)
                         ccctaaaGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCtggtgg                                                                                                                                                                                                                         >  1:269328/5‑251 (MQ=255)
                          cctaaaGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCtggtgg                                                                                                                                                                                                                         >  1:348859/4‑250 (MQ=255)
                            taaaGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGACGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCtggtgg                                                                                                                                                                                                                         >  1:71762/2‑248 (MQ=255)
                            taaaGTGATGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCtggtgg                                                                                                                                                                                                                         >  2:35280/2‑248 (MQ=255)
                               aGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTc                                                                                                                                                                                                                   <  1:87298/251‑1 (MQ=255)
                                 tGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGCGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTc                                                                                                                                                                                                                   >  1:614445/1‑249 (MQ=255)
                                        acacTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTgg                                                                                                                                                                                                          <  1:339890/251‑1 (MQ=255)
                                         cacTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCa                                                                                                                                                                                                                  <  2:224577/242‑1 (MQ=255)
                                         cacTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCa                                                                                                                                                                                                                  >  1:224577/1‑242 (MQ=255)
                                         cacTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGCGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGCCCTCCGGCAGCTCCGTGATGGCGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGCTTCAGGCGATGGt                                                                                                                                                                                                         >  2:594237/1‑251 (MQ=255)
                                          acTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTc                                                                                                                                                                                                        >  1:259083/1‑251 (MQ=255)
                                             gTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGAt                                                                                                                                                                                                     >  1:665018/1‑251 (MQ=255)
                                                 cGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACt                                                                                                                                                                                                 >  1:458909/1‑251 (MQ=255)
                                                         cGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCgatgat                                                                                                                                                                                          <  1:360637/250‑1 (MQ=255)
                                                          gtggcgtgtGGATGGCGATGCCGACTACCGTTGTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTGCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCTGTGGCTGGTGGAGATTCAGGCGCTGGTCGATCACGCGATGATgg                                                                                                                                                                                        <  2:221082/243‑1 (MQ=255)
                                                           gTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGc                                                                                                                                                                                       <  1:35280/251‑1 (MQ=255)
                                                           gTGCTTTTGGATGGCGATGCCGACGCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCGTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGc                                                                                                                                                                                       <  1:416438/251‑1 (MQ=255)
                                                                tttGGAGGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGaat                                                                                                                                                                                   <  2:259083/250‑2 (MQ=255)
                                                                 ttGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGaatc                                                                                                                                                                                  <  1:580576/250‑3 (MQ=255)
                                                                      gggCGATGCCGACTCCCGTTTGCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCAcgcc                                                                                                                                                                            <  1:607631/250‑1 (MQ=255)
                                                                       ggCGATGCCGATTCCCGTTTTCGCCCTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGcgccg                                                                                                                                                                           <  2:508403/251‑1 (MQ=255)
                                                                           aTGCCGACTCCCGTTTGCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTg                                                                                                                                                                       <  1:540012/251‑1 (MQ=255)
                                                                              ccGACTCCCGTGTTCGCACTTTACGCAGCAATAAAACCCGCGTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGGCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCa                                                                                                                                                                    <  2:162822/251‑1 (MQ=255)
                                                                               cGACTCTCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAg                                                                                                                                                                   <  1:246410/251‑1 (MQ=255)
                                                                                 agtCCCGTTTTCGCCATTTACGCAGACATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGGGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTACGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGGCGATCACTCGATGATGGAGAATCCGCGCCGCGTGGCAGTg                                                                                                                                                                 <  2:418580/249‑1 (MQ=255)
                                                                                 aCTCACGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGCGAATGAGCTTGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCCGCAACCCTTCGGCAATTATCTTAAGTCGCGGCGGTGAAGTGACCTCCGGCAGCTCCGGGATGGTGGTGTGGGAAGGACAGCGTCCGTTGCTTATGGAGCTTCAGGAGCTTGTCGATCACTCGATGCTGGCGAAACCACGCCGCGTGGCGGTg                                                                                                                                                                 >  2:499498/1‑251 (MQ=255)
                                                                                    cccGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTgggg                                                                                                                                                              <  1:377265/251‑1 (MQ=255)
                                                                                      cGTTTTCGCACTATACGCAGCCATAAAACCCGCTTCGGCTCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCGGCGTGAAGTCAGCAACCCTTCGGCAAGTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTTGGAAGGAACGCGTCCGTTTCTGGTGGAGATTCAGGCGCTGGTCGATCACGCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCt                                                                                                                                                            <  2:407205/251‑1 (MQ=255)
                                                                                           tcgaactttacgaggCCAGAAAAACCGCTGCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAAGTTTCTGAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTTCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGGTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGaa                                                                                                                                                        <  2:198512/236‑1 (MQ=255)
                                                                                               actttacgcAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTGCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACaaaa                                                                                                                                                   <  1:92529/245‑1 (MQ=255)
                                                                                                   tacgcAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGTGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGTTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGa                                                                                                                                                         <  1:501565/239‑1 (MQ=255)
                                                                                                    acgcAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGCGGCTGGCACAAAACCGTc                                                                                                                                              >  2:33711/2‑251 (MQ=255)
                                                                                                       cAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGACAGCAACCCTTCGGCAATTTTCTTAAGACGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGATGGAACGCGTCCGTTGATGGTGGAGGTTCCGGCGCTGGTCGATCACTCGATGATGGCGAATACGCGCCGCAGGGCAGTGGGGCTGGGACAAAACAGTCt                                                                                                                                             >  2:606953/1‑249 (MQ=255)
                                                                                                            aTAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTGCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCGCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATc                                                                                                                                      <  2:578348/251‑1 (MQ=255)
                                                                                                              aaaaaCCGCTTCGGCGCTGTGAATGAGCTGGGCGGCTTCGCGATGACCGCACAGGGGCTGCGTGAAGTCATCGACCCTTCGGCAATTTTCTGAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGAGCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATcc                                                                                                                                     <  2:458909/250‑1 (MQ=255)
                                                                                                              aaaaaCCGCTTCGGCGCGGGAAAGGAGCGGGGCGTCTGCGCGATGACCGAACAGGGGCTGCGGTAAGTCAGCAACCCTGCGGCAATTGTCTTAAGCCGCGGAGAGGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGATAGTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCct                                                                                                                                    <  2:351375/251‑1 (MQ=255)
                                                                                                               aaaaCCGCTGCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCctg                                                                                                                                   <  1:317887/251‑1 (MQ=255)
                                                                                                                  aCCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAAAAAAACGGTCTGGCCATCctgcgg                                                                                                                                >  2:285628/1‑249 (MQ=255)
                                                                                                                          ggCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCtgtctt                                                                                                                        <  1:322922/251‑4 (MQ=255)
                                                                                                                            cgcgGTGAATGAGCTGGTCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGc                                                                                                                      <  1:519503/251‑1 (MQ=255)
                                                                                                                            cgcgGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGCACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGc                                                                                                                      >  1:480991/1‑251 (MQ=255)
                                                                                                                              cgGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTCATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCAc                                                                                                                    <  1:348741/251‑1 (MQ=255)
                                                                                                                                 ggAATGAGCTGGGCTTCGTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGt                                                                                                                 <  2:558923/250‑1 (MQ=255)
                                                                                                                                    aggAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGGGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGGGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACGCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCAc                                                                                                              <  2:310367/249‑1 (MQ=255)
                                                                                                                                           gggCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTc                                                                                                       <  1:355585/251‑1 (MQ=255)
                                                                                                                                            ggCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTGCGGCAATTTTCTGAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGGCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAAGCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCt                                                                                                      <  2:338213/251‑1 (MQ=255)
                                                                                                                                               gTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGc                                                                                                    >  1:18979/1‑250 (MQ=255)
                                                                                                                                                 cTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCaa                                                                                                  >  1:550322/1‑250 (MQ=255)
                                                                                                                                                        aTGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCCCTCGATGATGGCGAATCCGCGCCGCGTGGCCGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGGCTTGCACCCTCACGGGGGTCTGCAACTGGCCg                                                                                           >  2:674603/1‑250 (MQ=255)
                                                                                                                                                         tGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGAt                                                                                         >  1:505363/1‑251 (MQ=255)
                                                                                                                                                            ccGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAg                                                                                      <  1:674603/251‑1 (MQ=255)
                                                                                                                                                            ccGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAg                                                                                      <  2:18979/251‑1 (MQ=255)
                                                                                                                                                               aaCAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGGTGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAAAGCGTCCGTTGCTGGTGGGGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGAGCCGCGGGGCAGTGGGGCTGGAACAAAACCGTCAGGCAATCCTGCTGGCTGTCTTGCACCGTAACGGTGGTCTGCAAATGGCCGATCAGGAt                                                                                   >  1:97828/1‑251 (MQ=255)
                                                                                                                                                                      gCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGGCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTtg                                                                            <  2:269328/251‑1 (MQ=255)
                                                                                                                                                                                  cAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGTCCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGATATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTc                                                                <  1:220549/250‑1 (MQ=255)
                                                                                                                                                                                  cAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTcg                                                               >  1:81565/1‑251 (MQ=255)
                                                                                                                                                                                        cccTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAg                                                                                                                                                                   >  1:343071/1‑146 (MQ=255)
                                                                                                                                                                                        cccTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGgcggcg                                                         <  1:130138/251‑1 (MQ=255)
                                                                                                                                                                                        cccTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATGCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCTCCGCGTGGCAg                                                                                                                                                                   <  2:343071/146‑1 (MQ=255)
                                                                                                                                                                                             cGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAg                                                    <  1:458050/251‑1 (MQ=255)
                                                                                                                                                                                                      tCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGCCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGGGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGa                                            >  1:34326/1‑250 (MQ=255)
                                                                                                                                                                                                         tAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAAcc                                        >  1:272108/1‑251 (MQ=255)
                                                                                                                                                                                                              cgcgGAGATGAAGTGACCTCCGGCCGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGCTTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTTGCAGTGGGGCTGGAACAAATCCGTCTGGCAATACTGCTGGCTGTCTTGCACCGACACGGTGGTCTGCCAATGGCCGTTCAGGCTGTGTTTGTGAACGTGGTCGGCGGCGGGCAGGTAACCGaacacagagc                                   >  2:123332/1‑243 (MQ=255)
                                                                                                                                                                                                                 ggagATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCg                                 >  1:438021/4‑250 (MQ=255)
                                                                                                                                                                                                                 ggagATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGAGGAGATTCAGGCGCTGGACGATCACTCGATGATGGCGAATCCGCGCCGCGAGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGGAACGGTGGTCTGCACCTGGCCGATCAGTGTGTGTTTGTGAACGTGGTCGGCGGCGTGCAGGTAACCGCAAcccgcgccgg                                >  2:3872/4‑243 (MQ=255)
                                                                                                                                                                                                                     aTGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACtt                             >  1:290367/1‑250 (MQ=255)
                                                                                                                                                                                                                       gcaGTGACCTCCGCCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGTTGGAGATTCAGGCGTTGGTCGATCACTCGATGATGGCGACTCCGCGCCTCGTGGCAGTGGTGCTGGACCAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATTGCCGATCAGGATGTGTTTGTTAACGTGTTCGGCGGCTTGTAGGTAACCGAAACCAGCGCCGACTTAGc                          <  2:348859/249‑1 (MQ=255)
                                                                                                                                                                                                                        aaGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCg                         >  1:79731/1‑251 (MQ=255)
                                                                                                                                                                                                                             cccTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTAct                   <  2:550322/250‑1 (MQ=255)
                                                                                                                                                                                                                              ccTCCGGCAGCTCCGTGATGGTGGTGTGGGCAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGGCTGCAAATGGCCGTTCAGGATGTGTTTGTGACCGTGGTCGGCGGCGTGAAGGTAACCGAAACCCGCGCCGACTTAGCGTTActg                  >  2:313528/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cggcaGCTCCGTGATGGTGGTGTGGGAAGTAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAATCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAATCCAGCGCCGACTTATCGGTActgcttg              <  1:131991/247‑3 (MQ=255)
                                                                                                                                                                                                                                     caGCTCAGTGATGGTGGGGTGGGAAGGAAAGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGAGGGGCTGGGACAAAACAGGCTGGCAATACTGCTGGCTGTCTTGCACCGGCACGGTGGTCTGCAACTGGCCGATCAGGATGTGTTTGTGAACGTAGTCGTCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTACGTTTACTGCGGGCg            >  2:516594/2‑250 (MQ=255)
                                                                                                                                                                                                                                      aGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTtg                                                                            <  1:614453/187‑1 (MQ=255)
                                                                                                                                                                                                                                        ctccgTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCAGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCGGGGTGTGTTTGTGAACGGTGGCGGCGGCGGGAAGGTAACCGAAAACCGCTACGACTTAGCGTTACTGCTGGGGATgg        >  2:454887/5‑251 (MQ=255)
                                                                                                                                                                                                                                              gATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCa  <  1:641681/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
TAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCA  >  NC_000913/4626573‑4627061

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: