Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,745,089 A→C 28.3% V84G (GTG→GGG)  trmD ← tRNA m(1)G37 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,745,0890AC28.3% 174.5 / 12.3 89V84G (GTG→GGG) trmDtRNA m(1)G37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (49/11);  new base C (3/25);  total (52/37)
Fisher's exact test for biased strand distribution p-value = 1.86e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCA  >  NC_000913/2744868‑2745310
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
ccGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCTACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCCCTTGTTCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCctt                                                                                                                                                                                                                            >  1:13443/1‑227 (MQ=255)
ccGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCttt                                                                                                                                                                                                                            >  1:1424077/1‑227 (MQ=255)
 cGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTg                                                                                                                                                                                                                           <  2:1396849/227‑1 (MQ=255)
     aTAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCCGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACGGCTAAATCACCTTTGCGCCt                                                                                                                                                                                                                     >  1:1614867/1‑228 (MQ=255)
     aTAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCt                                                                                                                                                                                                                     >  2:1730034/1‑228 (MQ=255)
         aCCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTTAc                                                                                                                                                                                                                 >  1:1721867/1‑228 (MQ=255)
                          aTCAGCGTCATTGCTGTTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTtgc                                                                                                                                                                                                 <  1:1639973/227‑1 (MQ=255)
                            cgccttCATTGCTGTAAACCCCCCACCACTGAAACCGTAATCGCCGTTTGCCCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGCCCGCACCCCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTtgctgc                                                                                                                                                                                              <  2:1612011/223‑1 (MQ=255)
                               cGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCTGCTTTTGCTGCAt                                                                                                                                                                                            >  1:709914/1‑227 (MQ=255)
                                 tCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTTTGGTGACAGATAAATCACCTTTCCGCCTTCACCCGCCCGGGCCTTTGCTGCATGaa                                                                                                                                                                                         >  1:471755/1‑228 (MQ=255)
                                      gCTGGTAACTCCCCCCCACTGAAACCGTACCCGCCGTTTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCg                                                                                                                                                                                    <  2:1142726/228‑1 (MQ=255)
                                      cttgttAACTCCCCACCAGTGAAAACGTAATCCCCAATTGACCATTTTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATCCCTTCGTAGCGACCGCCCACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCg                                                                                                                                                                                    <  2:409375/223‑1 (MQ=255)
                                       ctggaaccccaccaccacGGAGAACGTAATCGCCAATTGCCCTTTCTTCGTCATTTTCGGTTTGGATCCCGCGTTCATCTATACCTTCGTAGCGACCCCACACCGGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGt                                                                                                                                                                                   <  2:1175700/219‑1 (MQ=255)
                                             accacccaccacGGAAACCGTAACGCCCGATGAGCCGTGCTTCGACAATGCGGGTTGGGACCCGGCGCTCATCTATCCCTTCGAGGCGACCGCACACCAGAATCAATGTTTGATTCGTTGCCAGGTCGCTGACGCCCGCTTGACAAAGCTGGCGTCCCTGGGGTGACAGATAAACCCCCTTGGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCa                                                                                                                                                                            <  2:285510/223‑1 (MQ=255)
                                              cTCACCTCCACGTCGATCGTACTCCCATATTGCCCTTTCATCATCATTTTCGGTTTGGATCACGCGCCCATCTTACCCTTCGTAGCGCCCGCACCCCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCaa                                                                                                                                                                           <  2:1421527/228‑1 (MQ=255)
                                                 accaccacTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGGTGTTACAGATAAATCCCCTTTGCGCCTTTGCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAg                                                                                                                                                                          >  2:1286707/1‑226 (MQ=255)
                                                     ccacTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTg                                                                                                                                                                     >  1:1194286/1‑227 (MQ=255)
                                                      cacTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTTTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGTTGCATGAAGGGCGTCCCGCCAAGGGtgcc                                                                                                                                                                    >  2:1244980/1‑225 (MQ=255)
                                                       ggactccggctAAGCGCCGTATGCCCCTTCTTCCTCAATTTGGGTTTGGAACACCGCCTTCTTTATTCCTTCGGAGCCACCGCACCCCAAAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGAACAAGCTTGCCGCCCTGTGGTGAAAGAAAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCa                                                                                                                                                                   <  2:1564914/217‑1 (MQ=255)
                                                       ctggaaaccGTAACCGCCGTTTGCCCATTCTTCGTCATTTCCGGTTTGGATCCCGCCCTCATCTATCCCTTCGTAGCGCCCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTCCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCa                                                                                                                                                                   <  2:1189086/219‑1 (MQ=255)
                                                       acTGAACACGTAATCGCTGATTGACCATTCTTCGTCAATTTCGGTTTGGATCCCGCGCTCATTTATACCTTCGTAGCGACCGCACCCCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCa                                                                                                                                                                   <  2:13416/227‑1 (MQ=255)
                                                        ttGAGAACGTAACCGCCGTTTGCCCTTTCTTCGTCAATTTCGGTTTGGATCACGCGCTTATCTATCCCTTCGTAGCGCCCGCACCCCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCa                                                                                                                                                                   <  2:376726/225‑1 (MQ=255)
                                                           agaACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGc                                                                                                                                                                    >  1:377348/1‑222 (MQ=255)
                                                                       ccGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCCAGGGTTGCACCATCATTATCATctcc                                                                                                                                                  >  2:1131545/1‑228 (MQ=255)
                                                                          aTTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTCCGCCTTCCCCCGCCCCGGCTTTTGCTGCAGGAATGGCTTCCCGCAAGGGTTGCACCACCCTTatccggctcct                                                                                                                                                >  1:1567696/1‑217 (MQ=255)
                                                                              ccctttcctcGACATTTTCGGTTGGACACACGCGTCTATCTTTCCTTTCGTCGCGCCCGCACCCCAAAATCAATTTTTGATGCTTTGCCAGTTCTCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCCCCATCATTAACATCCCCGGTccgc                                                                                                                                           <  2:631998/224‑1 (MQ=255)
                                                                                      gctaaagtgcgggTGTGAATCGCGCGCACACCTAGCCGTACGGAGCGACCGCACCACAAAATCAATTTGTGATTGTTTCCCAGTCCGTTGCCGCCCGTTGGAACAAGCTTGCGTCCCTGGGGTACCAGAAAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTa                                                                                                                                    <  2:237871/216‑1 (MQ=255)
                                                                                       cGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGc                                                                                                                                                                    >  1:574177/1‑194 (MQ=255)
                                                                                        gTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCa                                                                                                                                                                   <  2:573241/194‑1 (MQ=255)
                                                                                             tttCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGa                                                                                                                             >  1:1175105/1‑227 (MQ=255)
                                                                                              ttCGGTTTGGATCTCGCGCTCACTTATCCTTCTGTGCCACCCCCACCCAAGATTCATTTTTTGATTCTTTGCCAGTTCGCTAACGCCCGTTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAACTCCCCTTTGCGCCTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGAt                                                                                                                            <  2:1173000/227‑1 (MQ=255)
                                                                                                   tttGGCTCACGTTCTTATCTATCCCTGCATAGCCACCGCACCCCAGTATAAATTTTTGATTCGTTGCCGTTTCGCTACCGCCGTTTTGATCAAGCTGGCGCCCCGTGGGGTAAAGATAAACCCCCTTTGCCCTTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGCCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCa                                                                                                                      <  2:1875227/228‑1 (MQ=255)
                                                                                                   tttGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGCCGATCGGcc                                                                                                                       >  1:1397668/1‑227 (MQ=255)
                                                                                                            cgcgCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGGCCacggtacggt                                                                                                             >  1:519729/1‑228 (MQ=255)
                                                                                                                  cATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATTGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCGGCCGTAAGGACGATCGTCCACGGCACGGTGCCgg                                                                                                        >  1:507567/1‑227 (MQ=255)
                                                                                                                       ataCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACTGTGCCGGTCATGc                                                                                                  >  1:34191/1‑228 (MQ=255)
                                                                                                                             tCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTTCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCACCATTAACATCCCCGGTccgacgc                                                                                                                                        >  1:691603/1‑184 (MQ=255)
                                                                                                                             tCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCTGCCGCGGCTTTTGCTGCATGCATGGCGTCCCCCAAGGGTTGCACCCTCCTTCACATCCCCGGTCCGCCGCCGCAAGGACGATCGTCCCCGGTACGGCGCCGGCCATGCGTCaa                                                                                             >  1:545294/1‑227 (MQ=255)
                                                                                                                             tCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGaa                                                                                             >  2:1183121/1‑227 (MQ=255)
                                                                                                                             tCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGaa                                                                                             >  2:1967017/1‑227 (MQ=255)
                                                                                                                             tCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGGGaa                                                                                             >  1:584859/1‑227 (MQ=255)
                                                                                                                                gTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggt                                                                                                                  >  1:1356767/1‑204 (MQ=255)
                                                                                                                                 tagcacccGCCCACCAGAATCATTTTTTGATTGTTTGCCAGTTCGCTGACGCCGGCTTGATCAAGTTGGCCTCCCTGTGGTGACAGAAAAACCCCTTTGGCGCTTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggta                                                                                                                 <  2:1354329/198‑1 (MQ=255)
                                                                                                                                     gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTCCGGTGCCGGTCATGCGGTAAGTcgcg                                                                                       >  2:1420369/1‑226 (MQ=255)
                                                                                                                                     gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTcgcg                                                                                       >  2:764065/1‑226 (MQ=255)
                                                                                                                                      aCCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTcgcg                                                                                       >  1:570627/1‑225 (MQ=255)
                                                                                                                                      aCCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGtt                                                                                                                                                                      >  2:1151820/1‑146 (MQ=255)
                                                                                                                                      aCCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTcgcg                                                                                       >  2:1021580/1‑225 (MQ=255)
                                                                                                                                       ccGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTGGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTg                                                                                                                                                                     <  1:1153887/146‑1 (MQ=255)
                                                                                                                                          caccacAGTATCTATTTTGGATGTCGTGCCCATTCCTCTGTCGCCCGTTTGAACAAGCTGGCTTCCTGGTGGTGACAGAAAAACCCCTTTTGCCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACt                                                                                 <  2:1395163/223‑1 (MQ=255)
                                                                                                                                          cacaCCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGCGGACTc                                                                                >  1:181539/1‑228 (MQ=255)
                                                                                                                                                    tCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTgg                                                                       >  1:1240688/1‑227 (MQ=255)
                                                                                                                                                      gaTTTTGGATCCTGTGCCAGTCCGCTGCCGCCCGGTTGCCCCAGCTCTCCCCCGTGTGGTGACACAAAAACCCCCTTGGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATg                                                                    <  2:506737/227‑1 (MQ=255)
                                                                                                                                                              aTTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGc                                                                                                                                                                    >  2:162721/1‑124 (MQ=255)
                                                                                                                                                               ttCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCa                                                                                                                                                                   <  1:163005/124‑1 (MQ=255)
                                                                                                                                                                                tGACGCCCGTTTGACAAGCTTGCGGTCCTTGTGGTGAAAGATAACTCCCTTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGcc                                            <  1:1367525/227‑1 (MQ=255)
                                                                                                                                                                                   cGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGACGTCCCGAGGACTCCAGCCATGGATGCTCAGCAGGCCATTTTTAACTGc                                             >  1:1660017/1‑223 (MQ=255)
                                                                                                                                                                                   cGCCCGCTTGATCAAGCTTGCGTCCCTGGGGTGACAGATAAACCACCTTGGCGCCTGCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGc                                        <  1:1209646/228‑1 (MQ=255)
                                                                                                                                                                                     cccGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCTTGCGTGAAGTCGCGAGGACTCCAGCTCTGGATCCTCAGCAGGCCATTTTTAACTGCCCGGCCa                                      >  2:889913/1‑228 (MQ=255)
                                                                                                                                                                                       cGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCCCAGCAGGCCCTTTTTAACTGCCCGGcccg                                     >  1:1599018/1‑225 (MQ=255)
                                                                                                                                                                                       cGCTTCATCCAGCTTGCTCCCTGGGTGTAACGATAAACTCCCTTTGGCCTCTTCCCCCGCCGCGGCTTTTGCTGCATGACTGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGt                                    <  2:89186/228‑1 (MQ=255)
                                                                                                                                                                                              tCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggta                                                                                                                 >  1:1449849/1‑144 (MQ=255)
                                                                                                                                                                                              acaggcttgcttCCGTTGGGTAACAGAAAACTCCCCTTCGCGCCTTCCCCCGCCGTGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCg                             <  2:181225/224‑1 (MQ=255)
                                                                                                                                                                                               cAAGGTTGGCTCCGTGTGGTGACAGAAAAACCCCCTTTCCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCCCCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggtac                                                                                                                <  2:1447257/144‑1 (MQ=255)
                                                                                                                                                                                                agccTTGCTCCCTTGTGGTGAAAGAAAAACCCCTTTGGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTa                           <  1:906587/225‑1 (MQ=255)
                                                                                                                                                                                                   cTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGCTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGCGAAGTCGCGAGGACTCCAGct                                                                           >  1:295569/1‑177 (MQ=255)
                                                                                                                                                                                                       cGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCCTTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGt                     >  1:1908847/1‑227 (MQ=255)
                                                                                                                                                                                                       cGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTa                    >  2:416874/1‑228 (MQ=255)
                                                                                                                                                                                                       cGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTCCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGGCCGGCCAGTTACCCCGTCATCGGTa                    >  1:1869843/1‑228 (MQ=255)
                                                                                                                                                                                                         gTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGt                     >  1:1022609/1‑226 (MQ=255)
                                                                                                                                                                                                             cTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCg                                                                                                                                                                               >  1:1281728/1‑69 (MQ=255)
                                                                                                                                                                                                             cTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTg                >  1:186765/1‑227 (MQ=255)
                                                                                                                                                                                                             cTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGCTACTTg                >  2:529419/1‑227 (MQ=255)
                                                                                                                                                                                                              tgtgGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGc                                                                                                                                                                             <  2:1279433/69‑1 (MQ=255)
                                                                                                                                                                                                              tgtgGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTgcg              >  2:1797750/1‑228 (MQ=255)
                                                                                                                                                                                                               gtggtaACAGATAAATCACCTTGGCGCCTCCCCCCGCCGCGGCTTTGGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACGTCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTgcg              <  1:1422917/227‑1 (MQ=255)
                                                                                                                                                                                                                gtgtgtgaaaaTAAATCCCCTTTGCCACCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTgcg              <  2:820045/218‑1 (MQ=255)
                                                                                                                                                                                                                 ggtgACAGATAAATCCCTTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCgg           <  1:446376/228‑1 (MQ=255)
                                                                                                                                                                                                                    gACATATAGATCACCTTGGCGCTTTCGCCCGCCGCGGCTTTTGCTGCATGAAGGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGGAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGaa         <  1:1133579/227‑1 (MQ=255)
                                                                                                                                                                                                                     aaCAGATAATTCCCTTTGCCGCTTTCCCCCCCCGCGGCTTTTGCTGCAGGACGGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGaa         <  2:1905462/226‑1 (MQ=255)
                                                                                                                                                                                                                      aCAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAgg                                                                                                                                                                         >  2:1430389/1‑66 (MQ=255)
                                                                                                                                                                                                                      aCAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGa                                                                                      >  2:281563/1‑149 (MQ=255)
                                                                                                                                                                                                                      aCAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACa       >  2:731845/1‑228 (MQ=255)
                                                                                                                                                                                                                       caaaaaaaTCCCCTTGGCGCTTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACCTCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACa       <  1:103041/222‑1 (MQ=255)
                                                                                                                                                                                                                       cAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAggg                                                                                                                                                                        <  1:1432954/66‑1 (MQ=255)
                                                                                                                                                                                                                       cAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAg                                                                                     <  1:282044/149‑1 (MQ=255)
                                                                                                                                                                                                                        aGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACCtt     >  1:1169652/1‑228 (MQ=255)
                                                                                                                                                                                                                        aGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACAGCCCGGCCAGTTACCCCGTAATCGGGAATTTCGCGGAGCAtt     >  2:1446474/1‑228 (MQ=255)
                                                                                                                                                                                                                           tAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGCAATTGCGCGGAACATTTCa  >  1:1558546/1‑228 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
CCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCA  >  NC_000913/2744868‑2745310

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: