Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,344,033 A→G 6.2% L105P (CTG→CCG)  yciT ← global regulator of transcription; DeoR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,344,0330AG 6.2% 370.9 / 10.9 114L105P (CTG→CCG) yciTglobal regulator of transcription; DeoR family
Reads supporting (aligned to +/- strand):  ref base A (53/51);  new base G (3/4);  total (56/58)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-02

Only 100 of 118 total aligned reads displayed.

CAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGACCC  >  NC_000913/1343774‑1344302
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                             
cAGCACGGCATTGACAAAATCGGAACGCATCATATCCCGACCGGAAAATCCATTTTAAGGTTGCCAGCCATCAATACCAATAAATGCTTTTCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCAGACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTa                                                                                                                                                                                                                                                                <  4:2570/275‑1 (MQ=255)
  gCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTCCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTAAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACCCCACCGACCAAAAAAACTTCAAAAGGCGCGTCTGTCAGCAAATGCACAATGTAGCTGCTGAACGTGCTGATAGTGACATTTTNNNTATGTCTGCACATAGGGCGGGCGAGTaa                                                                                                                                                                                                                                                               >  3:142482/1‑273 (MQ=255)
                  aTCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGCCATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGATGCCAt                                                                                                                                                                                                                                              >  1:19690/1‑275 (MQ=255)
                  aTCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTATTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGCGCAAAATAACTTCACACGCCCCCTCTTTCAGCAACTACCCGAAGTAGCTGCAGACCGTGATGATAGATTCATTTTTCTTCGCCTTGGCCAGAATCCGGGCGAGTAGGGCCTTGCTGCTTCCAt                                                                                                                                                                                                                                              >  4:416514/1‑275 (MQ=255)
                       tACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCg                                                                                                                                                                                                                                         >  2:290415/1‑275 (MQ=255)
                                                                     aacAATACCAATAAATGCGTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCgg                                                                                                                                                                                           <  2:495842/273‑1 (MQ=255)
                                                                                  aaTGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCtt                                                                                                                                                                              <  2:136719/275‑1 (MQ=255)
                                                                                          gCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGta                                                                                                                                                                       >  2:33035/1‑274 (MQ=255)
                                                                                           cTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGAGACATTTTTCTTCTGCTCGCCCAGCGTCCGGGCGGGTAGGGCATTGCTGCTNCCATTTTCGATCAAGGTGGTTTCGCCCGGTTGACCCAGCGACGGGGCAAACCTGGAGCGTTCACGGTGCAggggct                                                                                                                                                                      >  3:116934/1‑269 (MQ=255)
                                                                                                        ggtggatacagctgCGTGTCAAAGGGCCAACCAGACTTTCGTTTGTTTTCTGGTAAACGACAACGAGAAAAATAGCTTCACAAGTCGCGGCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTGCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATa                                                                                                                                                        <  4:531485/262‑1 (MQ=255)
                                                                                                        gTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCAc                                                                                                                                                                                  >  2:29126/1‑249 (MQ=255)
                                                                                                         ttGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCtt                                                                                                                                                                              <  4:33394/252‑1 (MQ=255)
                                                                                                          tGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATAc                                                                                                                                                       <  4:290414/274‑1 (MQ=255)
                                                                                                                 aCTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCACAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCag                                                                                                                                                                                                                                                                              >  8:1086461/1‑148 (MQ=255)
                                                                                                                  cTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCaga                                                                                                                                                                                                                                                                             >  5:633652/1‑148 (MQ=255)
                                                                                                                   tgtcGTGTCAAAAGGCCAACACTACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCCCAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGAc                                                                                                                                             <  4:333080/272‑1 (MQ=255)
                                                                                                                   tGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCTGATAGTTGCTCATCATACGGGTTTCGa                                                                                                                                              >  2:419456/1‑274 (MQ=255)
                                                                                                                   tGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGAc                                                                                                                                             >  2:482848/1‑275 (MQ=255)
                                                                                                                   tGGCGTGTCAAAGGGCCAACCATACTTGCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAATGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGAc                                                                                                                                             <  1:142482/275‑1 (MQ=255)
                                                                                                                     ggttgtCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTGTCAGCAAATGCGGTATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGt                                                                                                                                           <  1:116934/272‑1 (MQ=255)
                                                                                                                      cGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTcc                                                                                                                                                                                                                                                                         >  7:217459/1‑148 (MQ=255)
                                                                                                                           cAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCg                                                                                                                                                                                                                                                                    <  7:819954/148‑1 (MQ=255)
                                                                                                                            aaaGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGa                                                                                                                                                                                                                                                                   <  8:217459/148‑1 (MQ=255)
                                                                                                                                  ccAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCaa                                                                                                                               >  4:236924/1‑274 (MQ=255)
                                                                                                                                       cATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtg                                                                                                                                                                                                                                                        >  6:614020/1‑148 (MQ=255)
                                                                                                                                          aCTTTCGCTTTTTTTCTGATACACGCAACCGGGCAAAATAACTTCACAAGGCGCGTCTTTCCGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACa                                                                                                                       <  2:484345/274‑1 (MQ=255)
                                                                                                                                              tCGCTTTTTTTCTGATACACGCCACCGAGCGAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCTTCACTATCAAGCGAAACAGCaaa                                                                                                                  >  2:513819/1‑275 (MQ=255)
                                                                                                                                               cgttttttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCGTTAAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCaaa                                                                                                                  <  2:182162/271‑1 (MQ=255)
                                                                                                                                               cGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAg                                                                                                                                                                           >  2:176205/1‑217 (MQ=255)
                                                                                                                                                gCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACCCGGCAAACGCGGCGAGTTCCCGCTTCAGCGTTTAGTTGCTCATGATACGGGTTTCGACGTCATCCCTATCAAGCCAAACAGAAAAGc                                                                                                                >  4:37424/1‑275 (MQ=255)
                                                                                                                                                  tttttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAtt                                                                                                                                                                                                                                             <  5:614020/148‑1 (MQ=255)
                                                                                                                                                   ttttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAttt                                                                                                                                                                                                                                            <  7:1086461/148‑1 (MQ=255)
                                                                                                                                                    tttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGATGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAATATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACCTCATACCTATCAAGCGAACCAACCAAGcccag                                                                                                            >  4:484344/1‑272 (MQ=255)
                                                                                                                                                            ttACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCTTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCAcggc                                                                                                    <  2:480802/274‑1 (MQ=255)
                                                                                                                                                            gatacagctcccGAGAAAAATAACTTCACGTGGGCCGTCTTTCAGCAAAGTCGGCTGGTAGCTGCTGACCGTGAGGATAGTGACATTTTTCTTCTGCTCGCCCCGTGTCCGGGCGATTAGGGCATTGCTGCTGGCATTTTCGATAAAGATGGTTTCGCGCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCAcggc                                                                                                    <  4:5489/265‑1 (MQ=255)
                                                                                                                                                            aacacacgacccGAGCAGAATACGTTGACAAGGGGCGTCTTTGAGAAAGTGGGCGATGTAGCTGCTGACCGGTATGATAGTGAGATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGGAGAGCATTGCTGCTGCCATTGTCGATAAAGCTGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTGCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCAcggc                                                                                                    <  4:2198/265‑1 (MQ=255)
                                                                                                                                                              acacGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCg                                                                                                                                                                                                                                                                        >  7:857173/1‑109 (MQ=255)
                                                                                                                                                              acacGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCg                                                                                                                                                                                                                                                                        <  8:857173/109‑1 (MQ=255)
                                                                                                                                                                 cGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCGGACTCTGGGCGAg                                                                                                                                                                                                                                                                  <  6:163238/112‑1 (MQ=255)
                                                                                                                                                                 cGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCGGACTCTGGGCGAg                                                                                                                                                                                                                                                                  >  5:163238/1‑112 (MQ=255)
                                                                                                                                                                 cGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGt                                                                                                                                                                                                                              >  5:314551/1‑148 (MQ=255)
                                                                                                                                                                      ccGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACt                                                                                          >  2:186503/1‑275 (MQ=255)
                                                                                                                                                                         aGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCg                                                                                                                                                                                                                      <  7:993891/148‑1 (MQ=255)
                                                                                                                                                                         aGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAg                                                                                       <  1:84494/275‑1 (MQ=255)
                                                                                                                                                                             aaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGa                                                                                                                                                                                                                                  >  6:336929/1‑132 (MQ=255)
                                                                                                                                                                                   cTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATGGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATa                                                                                                                                                                                                                                      >  6:717266/1‑122 (MQ=255)
                                                                                                                                                                                     tCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCgg                                                                                                                                                                                                                                                                       <  5:210054/87‑1 (MQ=255)
                                                                                                                                                                                     tCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCgg                                                                                                                                                                                                                                                                       >  6:210054/1‑87 (MQ=255)
                                                                                                                                                                                           ggCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTTCTCATCATACGGTTTTCGACGTCATCACTTTCAAGCGAAACAGCAAAGCCCTGTTCACGGTGGAGGTCACTCACTTTTTCGGatgtgtctag                                                                     >  4:437199/1‑271 (MQ=255)
                                                                                                                                                                                                gTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCaa                                                                                                                   >  2:200946/1‑224 (MQ=255)
                                                                                                                                                                                                    ttCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGcc                                                                                                               <  4:200945/224‑1 (MQ=255)
                                                                                                                                                                                                      cAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTACTCGCCCGGACTCTGGGCGAg                                                                                                                                                                                                                                                                  <  7:618120/75‑1 (MQ=37)
                                                                                                                                                                                                      cAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTACTCGCCCGGACTCTGGGCGAg                                                                                                                                                                                                                                                                  >  8:618120/1‑75 (MQ=37)
                                                                                                                                                                                                       tgCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGa                                                                                                                                                                                        >  5:535576/2‑148 (MQ=255)
                                                                                                                                                                                                            aTGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGa                                                                                                                                                                                                                                                                   >  8:452979/1‑68 (MQ=255)
                                                                                                                                                                                                            aTGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGa                                                                                                                                                                                                                                                                   <  7:452979/68‑1 (MQ=255)
                                                                                                                                                                                                                 cgATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTAACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCCCTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACCCAGTTTTTCAAGTGGGTGTAACTCCTGGCGAAGGGTAACTTc                                                >  4:322960/1‑274 (MQ=255)
                                                                                                                                                                                                                   aTGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGAt                                                                                                                                                                                                                                       >  5:640175/1‑89 (MQ=255)
                                                                                                                                                                                                                   aTGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGAt                                                                                                                                                                                                                                       >  5:245210/1‑89 (MQ=255)
                                                                                                                                                                                                                   aTGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGAt                                                                                                                                                                                                                                       <  6:245210/89‑1 (MQ=255)
                                                                                                                                                                                                                   aTGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGAt                                                                                                                                                                                                                                       <  6:640175/89‑1 (MQ=255)
                                                                                                                                                                                                                   aTGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCa                                                                                                                                                                            >  6:226375/1‑148 (MQ=255)
                                                                                                                                                                                                                     gTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGa                                                           <  2:322961/259‑1 (MQ=255)
                                                                                                                                                                                                                         ctgctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTg                                                                                                                                                                                                                  >  5:614783/1‑104 (MQ=255)
                                                                                                                                                                                                                         ctgctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTg                                                                                                                                                                                                                  <  6:614783/104‑1 (MQ=255)
                                                                                                                                                                                                                         ctgctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGtat                                                                                                                                                                      <  6:535576/148‑1 (MQ=255)
                                                                                                                                                                                                                         ctgctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATTTGCACGGCGGGGGTCACTCAGTTTTTCGCGTTTGTTGAGTTCCTGGCCAACGTTACCTTCAGAAACTcc                                       >  4:534147/1‑275 (MQ=255)
                                                                                                                                                                                                                                 cGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgctg                                                                                                                                                                                                                                                  <  8:981040/64‑1 (MQ=255)
                                                                                                                                                                                                                                 cGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgctg                                                                                                                                                                                                                                                  >  7:981040/1‑64 (MQ=255)
                                                                                                                                                                                                                                 cGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCtc                                                                                                                                                              <  5:226375/148‑1 (MQ=255)
                                                                                                                                                                                                                                  gTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCa                                                                                                        <  4:108229/201‑1 (MQ=255)
                                                                                                                                                                                                                                   tgatgaTAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGAt                                                                                                                                                                                                                                 <  8:395083/79‑1 (MQ=255)
                                                                                                                                                                                                                                   tgatgaTAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGAt                                                                                                                                                                                                                                 >  7:395083/1‑79 (MQ=255)
                                                                                                                                                                                                                                    gatgatAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCgg                                                                                                                                                                                           >  6:443594/1‑116 (MQ=255)
                                                                                                                                                                                                                                    gatgatAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCgg                                                                                                                                                                                           <  5:443594/116‑1 (MQ=255)
                                                                                                                                                                                                                                           gTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCAc                                                                                                                                                                                  >  1:10429/1‑118 (MQ=255)
                                                                                                                                                                                                                                             gACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGcc                                                                                                                                                                                                                        <  6:1055583/78‑1 (MQ=255)
                                                                                                                                                                                                                                             gACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGcc                                                                                                                                                                                                                        >  7:183780/1‑78 (MQ=255)
                                                                                                                                                                                                                                             gACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGcc                                                                                                                                                                                                                        <  8:183780/78‑1 (MQ=255)
                                                                                                                                                                                                                                             gACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCgg                                                                                                                                                                                           <  7:990500/107‑1 (MQ=255)
                                                                                                                                                                                                                                             gACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTg                                                                                                                                                                 <  7:1075000/133‑1 (MQ=255)
                                                                                                                                                                                                                                             gACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTg                                                                                                                                                                 >  8:1075000/1‑133 (MQ=255)
                                                                                                                                                                                                                                             gACATTTTTCTGCTGCTCGCCCCGATTCCGTGCGAGTAGGGCATTTCTGTTGACATTTTCGCTTAAGATTGTTTCGCCCTGTTGAACCTGCGACGCGGCAAACTCGGCGAGTTCACGCTTCATCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCATAAACTCCAGTGGCTTTTGCCAGATCgg                   <  2:281133/275‑1 (MQ=255)
                                                                                                                                                                                                                                              aCATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCNTACGGGTTTCGACGTCATCACTATCAAGCGAAACACCAAAGCCATGTGCACGGAGGAGGTAAATCAGTTTTTCGNGTGTGTTGAGATCCTTGCGAATGGGAAATTNNAGAACTACAGTGGCcttttgaagatc                     >  3:73894/1‑264 (MQ=255)
                                                                                                                                                                                                                                               cATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCtt                                                                                                                                                                              <  3:10429/118‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGt                                                                                                                                                                                                                              >  7:1118656/1‑67 (MQ=255)
                                                                                                                                                                                                                                                  ttttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGt                                                                                                                                                                                                                              <  8:1118656/67‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCa                                                                                                                                                                                   <  6:807773/110‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCa                                                                                                                                                                                   >  5:807773/1‑110 (MQ=255)
                                                                                                                                                                                                                                                  ttttCTTCTGCTCGCCCAGAGGCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAg                        <  1:73894/265‑1 (MQ=255)
                                                                                                                                                                                                                                                        tctGCTCGCCCAAGGTCCGGGGGAGTAGGGAATTGCTGTTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAAcc        <  4:25698/275‑1 (MQ=255)
                                                                                                                                                                                                                                                        gcggccctcccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAAcc        <  4:33034/267‑1 (MQ=255)
                                                                                                                                                                                                                                                         ctacTCGCCCGGACTCTGGGCGAGCAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGa                                                                                                                                                                                                                 <  8:418751/70‑1 (MQ=25)
                                                                                                                                                                                                                                                              cGCCCGGACTCTGGGCGAGCAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACg                                                                                                                                                                                 <  6:63851/100‑1 (MQ=255)
                                                                                                                                                                                                                                                              cGCCCCGAGTCCGGGCGAGTAAGGCATTGCTGGTGCCATGTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGAccc  <  2:186380/275‑1 (MQ=255)
                                                                                                                                                                                                                                                              cGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATc                                                                                                                                 >  7:228747/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                cccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTg                                                                                                                                                                 >  7:563284/1‑114 (MQ=255)
                                                                                                                                                                                                                                                                 ccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTc                                                                                                                                                                                    >  7:649183/1‑94 (MQ=255)
                                                                                                                                                                                                                                                                 ccAGAGTCCGGACGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTc                                                                                                                                                                                    <  8:649183/94‑1 (MQ=255)
                                                                                                                                                                                                                                                                   agagTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCg                                                                                                                                               <  8:865979/129‑1 (MQ=255)
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                             
CAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGACCC  >  NC_000913/1343774‑1344302

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