Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 187,969 C→A 7.8% D228Y (GAT→TAT)  glnD ← uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913187,9690CA 7.8% 254.1 / 8.8 78D228Y (GAT→TAT) glnDuridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (33/38);  new base A (2/5);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATAACCCTTCCAGCATGCACTCTTCAAGCGTGCGCACGCTATGACCGACTTCCAGCTTTACATCCC  >  NC_000913/187688‑188256
                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                               
gcgTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCaa                                                                                                                                                                                                                                                                                           >  1:345403/1‑288 (MQ=255)
 cgTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAg                                                                                                                                                                                                                                                                                          >  1:44131/1‑288 (MQ=255)
        cgcgGAGGGAACCCTTAACAATCCGGTGGTCCGGTTCGTTACGTTCACCAATGGAATGCAGACGCGGGGCGACGGTAAGCGGGCGATGGATAAAAAGTGGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATAGGCAGACATTCGTTTAATTCCGCACGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAgg                                                                                                                                                                                                                                                                                  <  2:350767/289‑1 (MQ=255)
           ggAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGgc                                                                                                                                                                                                                                                                                <  1:578145/288‑1 (MQ=255)
                       tcatcaTCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTc                                                                                                                                                                                                                                                                   <  1:260562/289‑1 (MQ=255)
                            atcCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGt                                                                                                                                                                                                                                                               >  1:462621/1‑288 (MQ=255)
                             tcCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGGGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTc                                                                                                                                                                                                                                                             <  1:732533/289‑1 (MQ=255)
                              ccGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAGTGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCa                                                                                                                                                                                                                                                            <  1:579706/289‑1 (MQ=255)
                                                aCCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCa                                                                                                                                                                                                                                          <  2:637507/289‑1 (MQ=255)
                                                    tCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGc                                                                                                                                                                                                                                            >  1:600504/1‑283 (MQ=255)
                                                    tCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGc                                                                                                                                                                                                                                            <  2:600504/283‑1 (MQ=255)
                                                    tCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGt                                                                                                                                                                                                                                      >  1:720418/1‑289 (MQ=255)
                                                    tCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGt                                                                                                                                                                                                                                      >  2:152623/1‑289 (MQ=255)
                                                    tCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGAAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGt                                                                                                                                                                                                                                      >  1:393502/1‑289 (MQ=255)
                                                              aaTTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCGTTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGAt                                                                                                                                                                                                                             <  2:119306/288‑1 (MQ=255)
                                                              aaTTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGAt                                                                                                                                                                                                                             <  2:78385/288‑1 (MQ=255)
                                                               aTTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGCGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTCGCTGGTGCCATGGTAACGCTGACgg                                                                                                                                                                                                                           >  1:317966/1‑289 (MQ=255)
                                                                ttCAGACGCTGGGCTACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCAGGTTCAAGGTTGTAGCTGGTGCCATGGACACGCTTACGgc                                                                                                                                                                                                                          >  2:156851/1‑289 (MQ=255)
                                                                  cAGACGCTGGGGGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGgcgc                                                                                                                                                                                                                        <  1:30414/289‑1 (MQ=255)
                                                                  cAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTGTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGgcgc                                                                                                                                                                                                                        <  2:30234/289‑1 (MQ=255)
                                                                  cAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGAAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGgcgc                                                                                                                                                                                                                        <  2:393502/289‑1 (MQ=255)
                                                                  cAGAAGCGGGGGGACGCTAAGCTGGCGATCGAATAACCGGCGATTATCGTAACGGCTGACGACAAGATGCAGGGCAAAGAGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCAGTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGgcgc                                                                                                                                                                                                                        <  2:629926/289‑1 (MQ=255)
                                                                     aCGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTTCCATGGTAACGCTGATGGCGCTgg                                                                                                                                                                                                                     >  2:32627/1‑289 (MQ=255)
                                                                     aCGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGgcgc                                                                                                                                                                                                                        <  2:626039/286‑1 (MQ=255)
                                                                     aCGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGgcgc                                                                                                                                                                                                                        >  1:626039/1‑286 (MQ=255)
                                                                       gCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTc                                                                                                                                                                                                                                                                   >  2:131282/1‑241 (MQ=255)
                                                                       gCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTc                                                                                                                                                                                                                                                                   <  1:131282/241‑1 (MQ=255)
                                                                         tGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGgttc                                                                                                                                                                                                                  <  1:156851/288‑1 (MQ=255)
                                                                           ggCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGgttctgt                                                                                                                                                                                                               >  1:705729/1‑289 (MQ=255)
                                                                                     aGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATTCAGACATTCGTTTAATTCCGCCCGTTCCGCTGAGGTTAAGAAGCCACACCCGACCATTTCATCCAGCAATGTTGCGCAAACATGACGGCGAGCCACCCATTGCAGAGGGTGGATATCGCGCAAGCCGTCCGGGCTGATTTAGATGTCCGGTTCAAGGTTGTATTCGGTGCCATGGTATCGCAGCTGGCGCTGGTTCTGTTTTTCAACtt                                                                                                                                                                                                     >  2:190046/1‑289 (MQ=255)
                                                                                         tgCGATCGAGTACCAGGCGATCAGCGTAACGGTTGACGACCTGAAGCAGGGCAGGGCGAAGCAGCCACAGTAGATGCGGACATTCGTATAAGTCCGCCGGCTCGGCTGAGGTTAAGACACCAAAACCGACCAGTTCATCCAAAGGCGATGCGGCAAAATGACGGCGGGCCAACCATTGGAGCGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACGTTcg                                                                                                                                                                                                  <  2:215160/287‑1 (MQ=255)
                                                                                          gCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGgtt                                                                                                                                                                                                                   >  1:423419/1‑270 (MQ=255)
                                                                                          gCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGgtt                                                                                                                                                                                                                   <  2:423419/270‑1 (MQ=255)
                                                                                              tCGAATAACAGGCGATTATCGGAACAGCAGACGACCAGATGCAGGGCAAAGCGAATAAGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCcgcg                                                                                                                                                                                             <  1:585440/288‑1 (MQ=255)
                                                                                                                           gacgacCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCg                                                                                                                                                               <  1:32627/289‑1 (MQ=255)
                                                                                                                                          ggCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCACCTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGtttt                                                                                                                                                >  2:197955/1‑289 (MQ=255)
                                                                                                                                          ggCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGtttt                                                                                                                                                >  1:583324/1‑289 (MQ=255)
                                                                                                                                             aaaGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTg                                                                                                                                              <  2:500989/288‑1 (MQ=255)
                                                                                                                                             aaaGCGAAGACGCCACAGAATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTg                                                                                                                                              <  2:720418/288‑1 (MQ=255)
                                                                                                                                              aaGCGAATACGCCCATATAGACGCAGACATACGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAAGGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCa                                                                                                                                            <  2:317966/289‑1 (MQ=255)
                                                                                                                                               aGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAg                                                                                                                                           >  2:116161/1‑289 (MQ=255)
                                                                                                                                                    aTACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGaa                                                                                                                                                                              >  2:90676/1‑249 (MQ=255)
                                                                                                                                                    aTACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGaa                                                                                                                                                                              <  1:90676/249‑1 (MQ=255)
                                                                                                                                                        gCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGaa                                                                                                                                  <  1:712216/289‑1 (MQ=255)
                                                                                                                                                                    gCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATc                                                                                                                      >  1:267944/1‑289 (MQ=255)
                                                                                                                                                                         cATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGcc                                                                                                                   <  1:152623/287‑1 (MQ=255)
                                                                                                                                                                          aTTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTc                                                                                                                                                                                                                                                                   <  1:9940/142‑1 (MQ=255)
                                                                                                                                                                          aTTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTc                                                                                                                                                                                                                                                                   >  2:9940/1‑142 (MQ=255)
                                                                                                                                                                          aTTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCt                                                                                                                                                                  <  2:50641/239‑1 (MQ=255)
                                                                                                                                                                          aTTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCt                                                                                                                                                                  >  1:50641/1‑239 (MQ=255)
                                                                                                                                                                          aTTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAAt                                                                                                                >  1:667232/1‑289 (MQ=255)
                                                                                                                                                                          aTTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAAt                                                                                                                >  2:566673/1‑289 (MQ=255)
                                                                                                                                                                           ttCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCaa                                                                                                                                                                                                                                                                                           <  2:621971/117‑1 (MQ=255)
                                                                                                                                                                           ttCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCaa                                                                                                                                                                                                                                                                                           >  1:621971/1‑117 (MQ=255)
                                                                                                                                                                           ttCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAAtt                                                                                                               >  1:88934/1‑289 (MQ=255)
                                                                                                                                                                                         cccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTc                                                                                                 >  1:181251/1‑289 (MQ=255)
                                                                                                                                                                                                            gaagCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAAAAGCGCAACATCGCCAAGTAATAAGCGGTATTCGATTAAATTGGTGGCCACg                                                                              >  2:125555/1‑289 (MQ=255)
                                                                                                                                                                                                                 cAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGTAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGAAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTTGCGTCGGTTaa                                                                         >  2:429007/1‑289 (MQ=255)
                                                                                                                                                                                                                  aaaaCCCACCAGTTCATCCAGCGATGTTGCGCCAAAATGACGTCGGGCCACCCATTGCAGGGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTaaa                                                                        <  2:44131/289‑1 (MQ=255)
                                                                                                                                                                                                                  aaaCCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTaaa                                                                        <  1:107672/289‑1 (MQ=255)
                                                                                                                                                                                                                    aCCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAAt                                                                       >  1:313240/1‑288 (MQ=255)
                                                                                                                                                                                                                        gACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGAt                                                                  <  1:390782/289‑1 (MQ=255)
                                                                                                                                                                                                                        gACCATTTCACCCAGCAGTGTTGCGCCAAAGTGCCGGCGGGCACACCATTGCAGAGGGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGAt                                                                  <  2:667232/289‑1 (MQ=255)
                                                                                                                                                                                                                          ccATTTCATCCAGCGATGTTGCGCCAAAAGGACGGCGGGCCACCCATTGCAGAGGGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATaa                                                                <  2:345403/289‑1 (MQ=255)
                                                                                                                                                                                                                                        gagtgtGCGCCAAAATGACGGCGGGCCAGCCATTGCAGAGTTTGGATAGCGCGCAAGCCGGCAGGGCTGCTTTTTAGGTATGGTTCAAGGTTTTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAAGCCTTCGCGGAAAATATGTTTGTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATAACCCTTCCAGCATg                                                   <  2:529933/283‑1 (MQ=255)
                                                                                                                                                                                                                                          tGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATAACCCTTCCAGCATGCa                                                 <  1:197955/288‑1 (MQ=255)
                                                                                                                                                                                                                                               cgcCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGc                                                                                                                                                                                                                                            <  1:502898/96‑1 (MQ=255)
                                                                                                                                                                                                                                               cgcCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGc                                                                                                                                                                                                                                            >  2:502898/1‑96 (MQ=255)
                                                                                                                                                                                                                                                        ggcggcggGCCACCCATTGCAGAGTGTGGATATCGAGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATAACCCTTCCAGCATGCACTCTTCAAGCGTgcg                                  <  2:57209/287‑1 (MQ=255)
                                                                                                                                                                                                                                                          cggcggGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATAACCCTTCCAGCATGCACTCTTCAAGCGTGCGCa                                <  2:462621/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aCCCATTGCAGTGTGTGGATTTCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAAGCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATACCCCTTCCAGCATGCACTCTTCAAGCGTGCGCACGCTATGAc                       <  1:190046/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                tgtgGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATAACCCTTCCAGCATGCACTCTTCAAGCGTGCGCACGCTATGACCGACTTCCAGCtt          >  2:107947/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                tgtgGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATAACCCTTCCAGCATg                                                   <  2:12777/248‑1 (MQ=255)
                                                                                                                                                                                                                                                                                tgtgGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATAACCCTTCCAGCATg                                                   >  1:12777/1‑248 (MQ=255)
                                                                                                                                                                                                                                                                                  tgGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTCGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATAAccc                                                             <  2:639320/236‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  tgGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTCGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATAAccc                                                             >  1:639320/1‑236 (MQ=255)
                                                                                                                                                                                                                                                                                       atCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATAACCCTTCCAGCATGCACTCTTCAAGCGTGCGCACGCTATGACCGACTTCCAGCTTTACATcc   <  2:88934/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        tCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATAACCCTTCCAGCATGCACTCTTCAAGCGTGCGCACGCTATGACCGACTTCCAGCTTTACATccc  >  1:577592/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                               
GCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAAAATATGTTTTTGCAGTTCGAGGAACAGCGCAACATCGCCAATTAATAAGCGGGATTCGATTAAATTGGTGGCGACGGTTAAATCCGATAACCCTTCCAGCATGCACTCTTCAAGCGTGCGCACGCTATGACCGACTTCCAGCTTTACATCCC  >  NC_000913/187688‑188256

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