Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,627,069 C→T R386R (CGC→CGT radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,627,0690CT95.3% 262.8 / NA 85R386R (CGC→CGTradADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (39/42);  minor base G (1/3);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGG  >  NC_000913/4626822‑4627309
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
gggAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGAc                                                                                                                                                                                                                                               <  1:541860/251‑1 (MQ=255)
    aggaaCGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGAc                                                                                                                                                                                                                                           <  1:295177/251‑1 (MQ=255)
     ggaaCGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGAcc                                                                                                                                                                                                                                          <  2:67028/251‑1 (MQ=255)
        aCGCGTCCGTTGCTGGTGGGGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCt                                                                                                                                                                                                                                       <  2:64639/251‑1 (MQ=255)
            gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGcc                                                                                                                                                                                                                                    >  1:12899/7‑250 (MQ=255)
                gttgttgggggAGAGTCAGGCGCTGGTAGATCACGCGAGGATGGCGAATCCGCTCCGCGTGGCAGTGTGGCTGCAACAAAACCGTCTGGCAATCCTGCGGGCTGTCTTGCACCGTCACGGTGGTCTGCAAAGGGCCGATCAGGATGGGTTTGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACGTAGCGTTACTGCTGGCGATGGTTTCCAGCCGGCGTGACAGACCGCTGCCACAgg                                                                                                                                                                                                                               <  2:236166/242‑1 (MQ=255)
                  tgctgGTGGATATGCAGGCGCTGGTCGATCTCTCGAGGATGGGGAATTCGCGCGGCGTGGCTGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGCCTTGGACCGTAACGGTGGTCTGCAAATGGCCGATCAGGATGTTTTTGTGAACGGTGTCGGCGTCGTGAAGGTAACCGAAACCAGCGCCGGCTTAGCGTTACTGCTGGCGATGGTTTCCAGACTGGGTGACGGACCGCTGCCTCAGGa                                                                                                                                                                                                                              <  2:296880/250‑1 (MQ=255)
                   gctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATc                                                                                                                                                                                                                            <  1:358132/251‑1 (MQ=255)
                       gtggAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACTAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGATACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCtggt                                                                                                                                                                                                                        >  1:396074/1‑251 (MQ=255)
                          gagaTTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTggtggt                                                                                                                                                                                                                     >  2:110846/1‑251 (MQ=255)
                            gaTTCAGGCGCTGGTCGATCCCTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGgtgt                                                                                                                                                                                                                   <  2:260382/251‑1 (MQ=255)
                             aTTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTGGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGtt                                                                                                                                                                                                                  <  1:481648/251‑1 (MQ=255)
                                 gggCGCTGGTCGATCACTCGATGAGGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCGGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGt                                                                                                                                                                                                              <  2:298548/250‑1 (MQ=255)
                                       tGGTCGATCACTCGCTGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTc                                                                                                                                                                                                        >  1:487516/1‑251 (MQ=255)
                                        ggTCGTTCACTCGATGAGGTCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAATCCGTCTGGCAAGCCTGCTGGCTGTCTTGCACCGTCACGGGGGTCTGCAAATGGCCGATCAGGATGTGTTTGGGAACGTGGTCGGCGGCGTGAAGGTAGCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGt                                                                                                                                                                                                         <  2:192886/249‑1 (MQ=255)
                                             aTCACTCGATGATGGCGAATCCGCGCCGCGAGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTg                                                                                                                                                                                                  >  1:473823/1‑251 (MQ=255)
                                                       gatggcgaatccggggcgcgTGGCAGGGGGGCTGGAACAAAACCGTCGGGCAATCCTGCTGGCTGTCTTGCTCCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCTTTACTGCTGGCGATGGTGTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGgag                                                                                                                                                                                         <  2:418419/234‑1 (MQ=255)
                                                       gatgGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGACCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAgggg                                                                                                                                                                                           >  1:450746/1‑248 (MQ=255)
                                                        atgGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGGGAAGGTAACCGAAACCAGCGCAGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGGGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGGGAAGTCGGGCGGGCAGGGgag                                                                                                                                                                                         >  2:541907/1‑249 (MQ=255)
                                                          ggCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATcc                                                                                                                                                                                     <  1:350882/251‑1 (MQ=255)
                                                            cgaatccgctgcgcgTGGAGGTGGTGCTGGAACAAAACCTTGTGGCAATCCTGCTGGCTGTCTTGCACCGTCCCGGTGGTCTGCAAATGGCCGATCAGGATGTGTTGGTGAACGTGGTCGGCGGCGTGCAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCGTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATACGc                                                                                                                                                                                   <  2:354240/240‑1 (MQ=255)
                                                                  cgcgccggGTTGCAGTGGGGCTGGAACAAAACCGTCTGGCAAGCCGGCTGGCGGTCTTGCATCGTCACGGTGGTCTGCAAAGGGCCGATCAGGATGTGTTTGTGAATGGGGTCGGCTGCGTGAAGGTATCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCGGCCACAGGATCGGGGGGTGTGTGGTTAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTg                                                                                                                                                                             <  2:585259/249‑1 (MQ=255)
                                                                            ggCAGTGGGGCTGGACCCAAACCGTCTGGCAATCCTACTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACCGACCGCTGCCACAGGTTCTGGTTGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGCTCCGCCCGGTGCCCAGCgggc                                                                                                                                                                   >  2:313269/1‑249 (MQ=255)
                                                                                          aaCAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTGGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttca                                                                                                                                                     <  2:384892/251‑5 (MQ=255)
                                                                                            cAAAACCGTCTGGCAATCCTGCGGGTTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttcaga                                                                                                                                                   <  2:362233/251‑7 (MQ=255)
                                                                                              aaaCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAg                                                                                                                                                 <  2:5724/251‑1 (MQ=255)
                                                                                               aaCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgc                                                                                                                                                >  1:55617/1‑251 (MQ=255)
                                                                                               aaCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgc                                                                                                                                                >  1:487810/1‑251 (MQ=255)
                                                                                                aCCGTCTGGCACTCCTGCTGGCTGTCGTGCCCCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcg                                                                                                                                               <  1:337828/251‑1 (MQ=255)
                                                                                                aCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGACACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcg                                                                                                                                               >  1:53969/1‑251 (MQ=255)
                                                                                                aCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACGGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcg                                                                                                                                               <  1:321767/251‑1 (MQ=255)
                                                                                                   gTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGTTCTGCAAATAGCCGAGCAGGATGTGTTTGTTAACGTGTGCGGCGGCGTGAAGGTAACCGAAACCTGCGCCTTCTTAGCGTTACTGCTGGCGATGGTTTCCAGCTGGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcggcg                                                                                                                                            <  1:610333/251‑1 (MQ=255)
                                                                                                         cAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAacac                                                                                                                                      <  2:284588/251‑1 (MQ=255)
                                                                                                              ctgctgGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGttt                                                                                                                                 >  1:161414/1‑251 (MQ=255)
                                                                                                                gctggctgTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTTCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGttt                                                                                                                                 <  1:541907/249‑1 (MQ=255)
                                                                                                                     ctgTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCgg                                                                                                                          >  1:586620/1‑251 (MQ=255)
                                                                                                                       tgctggCACCGTCACGGTGGGCTGCAAAGGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGGGTGAAGGTAACCGAAACCAGCGCCGACTCAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTTACAGACCGCTGCCACAGGATCTGGTGGTGTGTGGTGAAGTCGGGCTTGCAGTGGAGATCCGCCCGGTGCCCAGCTTTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCggg                                                                                                                         <  1:302464/245‑1 (MQ=255)
                                                                                                                            gcgccGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCTGCGCCGACTGAGCGTTACTGCTGGCGATGGTTGCCAGCCTGCGTTACAGTCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTTCCCAGCTGTCATGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTg                                                                                                                   <  1:555816/248‑1 (MQ=255)
                                                                                                                            gCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGAtt                                                                                                                    <  1:459261/250‑1 (MQ=255)
                                                                                                                                 gTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGTGATTGTTCcg                                                                                                              >  2:364362/1‑251 (MQ=255)
                                                                                                                                    aCGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCcggcg                                                                                                           <  2:487516/251‑1 (MQ=255)
                                                                                                                                     cGGTGGTCGGCGAATGGTAGAACAGAATATGTTTCTGAACGTGGAAGGAGGCGAGAAGGTAACCGAAACAAGAGCCAAATTAGCGTTACTGCTCGCGATGGTTTCCAGCCTGAGTGACAGACCGCTGCAACAGGATATGGTGGTGTTTGGCGCAGTAGGGCTGGCAGGGGAGCTCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCtggcgg                                                                                                          >  2:584517/1‑251 (MQ=255)
                                                                                                                                              gCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGACACACGGCTTTCGCCGGGCGATTGTTCCGGCGGCTAACGtac                                                                                                 >  2:48942/1‑249 (MQ=255)
                                                                                                                                              gCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTGTGGTGAAGTCGGGCTGGCAGGGGAGATTGGCCCGGTTCCCAGCGGTCAGGAACGAATTTCAGATGCGGCTAAACACGGTTTTCGCCGGGCGATGGTTCCGGCGGCTAACGta                                                                                                  <  1:242988/250‑2 (MQ=255)
                                                                                                                                                 aaTGGCCGATCCGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCg                                                                                               >  1:251392/1‑250 (MQ=255)
                                                                                                                                                  aTGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACggaccga                                                                                              >  2:280320/1‑244 (MQ=255)
                                                                                                                                                    ggCCGATAAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGGGTTTGGTGAAGTCGGGCTGGCAGGGTAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGAGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTCACGTACCgcaa                                                                                            >  2:534463/1‑247 (MQ=255)
                                                                                                                                                        gATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAgc                                                                                       <  1:461587/251‑1 (MQ=255)
                                                                                                                                                         aTCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGAACCGAAAAAAgcg                                                                                      >  2:218230/1‑251 (MQ=255)
                                                                                                                                                                   gtTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGa                                                                            <  2:165790/251‑1 (MQ=255)
                                                                                                                                                                    tttGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTATCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGAGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAACAAGCGCCGGAAGGGAt                                                                           >  1:51891/1‑251 (MQ=255)
                                                                                                                                                                      tgtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGc                                                                         >  2:39741/1‑251 (MQ=255)
                                                                                                                                                                           aCGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGAtt                                                                    <  1:295407/251‑1 (MQ=255)
                                                                                                                                                                              tGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGAttttt                                                                 >  2:162447/1‑251 (MQ=255)
                                                                                                                                                                                  cggcggcggGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCg                                                             <  2:141703/251‑1 (MQ=255)
                                                                                                                                                                                  cggcggCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGACGGAAGGGATGCAGATTTTTGGCg                                                             >  2:202727/1‑251 (MQ=255)
                                                                                                                                                                                  cggcggCGTGAAGGTAACCGAAACCAGAGCCGACTTAGTGTTACGGCTGGGGATGGTTTACAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGAGCCGCCCGGTGCCCAGCGGTCAGGAACGAATTGCAGAAGCGGCGAAACACGGTTTTCGCCGGGCTATTGTTCCGGCGGCTAACGTTCCGAAAAAAGCGCCGGAAGGGAGGCAGATTTGTGGCg                                                             <  2:579741/251‑1 (MQ=255)
                                                                                                                                                                                       gcgTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTaaa                                                        <  2:546089/251‑1 (MQ=255)
                                                                                                                                                                                         gTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGAGGATTTTTGGCGTTaaaa                                                       >  2:616319/1‑250 (MQ=255)
                                                                                                                                                                                              ggTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAActct                                                 >  1:607675/1‑251 (MQ=255)
                                                                                                                                                                                                tAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTcc                                               >  1:331568/1‑251 (MQ=255)
                                                                                                                                                                                                   ccGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAc                                            <  1:39741/251‑1 (MQ=255)
                                                                                                                                                                                                      aaacaagcgCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGAGTGACAGACAGCTGCCACATGATCTGGGGGTGCTTGGTGAAGTCGGGCTGGCAGGTGCGATACGCCCGGAGCCCAGAGGTCAGGAACGAATTTCAGAAGCGGCGAAACAAGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAATGTACCGAAAAAAACGCCGGAAGGGAAGGCGATTTGTGGCGTTAACAAACTCTCCGAcgcg                                         >  1:299648/9‑251 (MQ=255)
                                                                                                                                                                                                        accagccccGACTTAGCTTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCt                                       <  2:338919/244‑1 (MQ=255)
                                                                                                                                                                                                             cgCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCg                                  >  1:231932/2‑251 (MQ=255)
                                                                                                                                                                                                               ccGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAACCTCTCCGACGCGCTTAGCgtg                                >  2:110918/1‑251 (MQ=255)
                                                                                                                                                                                                                cgtctgagcgtggatgctgGTGATGGTTCCCGGCCGGCGGGACAGACCTCTGGCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGTAGGGGAGATCCTCCCGGGGCCCAGCGGTCAGGAACGAATGTCAGAAGCGTCGAAACACGGTTTTCGCCGGGCGATTGGTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGAGTTTTGGCGTTAATAAACTCTGCGACGCGCTTAGCgtgt                               <  2:558703/237‑1 (MQ=255)
                                                                                                                                                                                                                cGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAGCGAATTTCAGAAGCGGCGAAACACGGTTTTCGGCGGGCGATTGTTCCGGCGGCTAGCGTACCGAAAAAAGCGCCGGAAGGGAGGCCGCTTTTTGGCGTTCAAAAACTCTCCGACGCACTTAGCgggt                               >  1:520155/1‑248 (MQ=255)
                                                                                                                                                                                                                 gACTTAGCGTTACTGCTGGCGATGGTTTCCATCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGtt                              <  1:406290/251‑1 (MQ=255)
                                                                                                                                                                                                                 gACTTAGCGTTACTGCTGGCGATGGTTTCCAGTCTGCGTGACAGACAGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGAGGCAGATTTTTGGCGTTAAAAAAATCTCCGACGCGCTTAGCGTGtt                              >  2:92442/1‑251 (MQ=255)
                                                                                                                                                                                                                     tcgcgctACTGCTGGCGCTGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGTCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGGGGCTAAAGGACCGGAAAAAGGGCAGGACGGGGTGCAGATTTTTGGCGTTAAACAAATCTCCGACGAGCGTAGCGTGGTcgcc                          >  2:450645/7‑249 (MQ=255)
                                                                                                                                                                                                                      aGCGTTACTGCTGGCGATGGTTTCCAGCCGGCGTGGCAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCTTTAAAAAACTCTCCGACGCGCTTAGCGTGTTcgac                          <  2:364232/250‑1 (MQ=255)
                                                                                                                                                                                                                       gCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTcgacg                         >  1:51972/1‑250 (MQ=255)
                                                                                                                                                                                                                        cGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGGGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCgacgac                       >  2:462900/1‑251 (MQ=255)
                                                                                                                                                                                                                         gtaaCTGCTGGCGATGGTTTCCATCCTGCGTGACAGGCCGCTGCCACAGGATCTGGTGGTGTTGGGTGAAGTCGGGCTGGGAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACTAATTTCAGAAGCGGCGAAACACGGTTTGCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCGTAGCGTGTTCGACGACt                      <  1:255521/248‑1 (MQ=255)
                                                                                                                                                                                                                          ttACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACtt                     >  2:323974/1‑251 (MQ=255)
                                                                                                                                                                                                                            aCTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTta                    <  1:323810/250‑1 (MQ=255)
                                                                                                                                                                                                                                  gtcGATGGTTTCCAGCCTGCGGGACAGACCGCTGCCCAAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAATCTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATgag             <  2:512365/249‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCGATGGTTTCCAGCCTGCGTGACAGACCGCGGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCTGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGTTTAGCGTGGTCGACGACTTATAATg               <  2:586620/249‑1 (MQ=255)
                                                                                                                                                                                                                                     gagggTTTCCAGCCTGGGTGACTGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATgaggaa          <  2:98189/248‑4 (MQ=255)
                                                                                                                                                                                                                                      aTGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATgaggaat         >  2:189968/1‑247 (MQ=255)
                                                                                                                                                                                                                                         gTTTCCAGCCGGCGTGACAGACCGCTGACACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGCGATACGCCCGGTGACCAGCGGTCCGGAACGAATTTCAGAAGCGGCGCAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGAACCGAAAACAGCGCCGGAAGGGAGGCAGAGTTTTGGCGTTAAAAAACTCTCCGACGCGCTTCGAGTGTTCGACGACTTATAATgaggaatacg      >  2:570988/1‑244 (MQ=255)
                                                                                                                                                                                                                                          tttCCAGACTGCGTGACAGACCGCTGCCACCGGACCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGCACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTAACGAAAAACGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACgg     >  2:60270/1‑251 (MQ=255)
                                                                                                                                                                                                                                            tCCAGCCTGCGGGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACgga    <  1:356892/250‑1 (MQ=255)
                                                                                                                                                                                                                                             ccAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACggagg  >  2:173334/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
GGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGG  >  NC_000913/4626822‑4627309

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