Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 867,717 G→C 5.0% L66F (TTG→TTC)  gsiA → glutathione ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913867,7170GC 5.0% 584.5 / 16.9 159L66F (TTG→TTC) gsiAglutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (75/76);  new base C (4/4);  total (79/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-02

Only 100 of 160 total aligned reads displayed.

ATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTC  >  NC_000913/867450‑867985
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                             
    aTCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcggg                                                                                                                                                                                                                                                                    <  2:406600/275‑1 (MQ=255)
    aTCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcggg                                                                                                                                                                                                                                                                    <  4:320162/275‑1 (MQ=255)
    aTCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGGATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGGAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACAAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcggg                                                                                                                                                                                                                                                                    <  2:301036/275‑1 (MQ=255)
        cgcCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGCTGCACTTGATGCCGGTAATGTGCTGGCGGTTGAAAACCTGAATATTGCCTTTATGCAGGAACAGCAGCAAATAGCTGCGGTCCGCAATCCCTCTTTTAGTCCGCAACGCGGTGAGACGCTGGCAACTGTTGGCCAATCCGGCTCCGGTAAGTCAGGTACTGCCGTTGCATTGATGCCCCTGTTGAAAcccgtcgcgta                                                                                                                                                                                                                                                                 >  4:153505/1‑264 (MQ=255)
             ggggCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGt                                                                                                                                                                                                                                                          >  2:484694/1‑275 (MQ=255)
                                           gTATCTATCGTAAAAAAGGGGACACCGTTGCCAAACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAgcggcg                                                                                                                                                                                                                            <  1:112684/275‑1 (MQ=255)
                                           gTATCTACCGTGAAAAAGGGGAAACCGTTNNNAAACAGTGATGAACTTGAGGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAgcggcg                                                                                                                                                                                                                            <  3:160066/275‑1 (MQ=255)
                                                    gTGAAAAAGGGGAGACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGa                                                                                                                                                                                                                   <  4:327526/275‑1 (MQ=255)
                                                     tGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGaa                                                                                                                                                                                                                  >  2:106787/1‑275 (MQ=255)
                                                     tGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTCGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCCCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTTTTGGCGAATCCGGCTCCGTTAAGTCAGTGACTCTGTTGGCATGGATGCGCCTGTTGGGACAGGCGGGG‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGTCGCCGTggcgg                                                                                                                                                                                                                   >  4:51131/1‑273 (MQ=255)
                                                                 cacCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAAc                                                                                                                                                                                                        >  2:255495/1‑273 (MQ=255)
                                                                                    tgaAGTTGTGGCCGGGAAGGGGCGGGCGGTGGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGca                                                                                                                                                                                    <  4:463076/274‑1 (MQ=255)
                                                                                    tgaACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGtt                                                                                                                                                                                                                                     >  4:171512/1‑225 (MQ=255)
                                                                                      aaCTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAgg                                                                                                                                                                                                                                                                        >  4:199930/1‑189 (MQ=255)
                                                                                        cTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAg                                                                                                                                                                                                                                 <  2:171513/225‑1 (MQ=255)
                                                                                         ttGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATg                                                                                                                                                                              >  4:183599/1‑275 (MQ=255)
                                                                                                  gTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGtt                                                                                                                                                                     >  2:174930/1‑275 (MQ=255)
                                                                                                         gCTGTCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGTTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTg                                                                                                                                                              <  2:487168/275‑1 (MQ=255)
                                                                                                                             atatTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGa                                                                                                                                                                                                                                                                             <  5:97633/145‑1 (MQ=255)
                                                                                                                                gtGCCTGTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTTTCTTTTAGTCTGCAACGCGGGGAGACGCTGGGAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGGTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGC‑GTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAg                                                                                                                                       <  2:76971/273‑1 (MQ=255)
                                                                                                                                   ccTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcgggc                                                                                                                                                                                                                                                                   >  2:269449/1‑149 (MQ=255)
                                                                                                                                   ccTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTTTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCACCATGTTCGCGTGGCGGAGTTGGCGAGGTTATTTTAGGAg                                                                                                                                    >  2:433806/1‑275 (MQ=255)
                                                                                                                                     tttATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcgggc‑g                                                                                                                                                                                                                                                                 >  8:330080/1‑148 (MQ=255)
                                                                                                                                       tATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGtt                                                                                                                                                                                                                                                              <  4:269448/149‑1 (MQ=255)
                                                                                                                                        aTGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAgc                                                                                                                                                                                                                                >  4:286910/1‑178 (MQ=255)
                                                                                                                                         tGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGttt                                                                                                                                                                                                                                                             <  7:330080/148‑1 (MQ=255)
                                                                                                                                            aGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAgcggcg                                                                                                                                                                                                                            <  2:286911/178‑1 (MQ=255)
                                                                                                                                             ggACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTCCAACAgg                                                                                                                                                                                                                                                                        <  5:390603/134‑1 (MQ=255)
                                                                                                                                                  agcagAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGtt                                                                                                                                                                                                                                                              <  6:428031/138‑1 (MQ=255)
                                                                                                                                                   gcagAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCCTGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGAGGAACTCGc                                                                                                                       >  1:105821/1‑272 (MQ=255)
                                                                                                                                                    cagAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAg                                                                                                                                                                                                                                                      >  5:49443/1‑144 (MQ=255)
                                                                                                                                                    cagAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAg                                                                                                                                                                                                                                                      <  6:49443/144‑1 (MQ=255)
                                                                                                                                                        aaaggccTGCGGTCCGCAATCCTGTTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCAGTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAAccc                                                                                                                <  2:263721/268‑1 (MQ=255)
                                                                                                                                                                 cGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCCGAAAGCTGCACAAATGCGCCATGTTCGCGGTGCCGATATGGCGGTGATATTTCCGGAGCCGATGGCATCCCCTAACCCGGTTTTTACt                                                                                                      >  4:260632/1‑275 (MQ=255)
                                                                                                                                                                   gTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACtg                                                                                                     <  2:308881/274‑1 (MQ=255)
                                                                                                                                                                      cGCAATCTCTCTTTCAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGGTTTAGTACAGTGCGATAAAATGATGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGACCCGATGACATCGCTGAACCCGGTATTTACTGTgg                                                                                                  >  2:190672/1‑275 (MQ=255)
                                                                                                                                                                       gCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAAt                                                                                                                                                                                                                                           >  5:1091887/1‑136 (MQ=255)
                                                                                                                                                                       gCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAAt                                                                                                                                                                                                                                           <  6:1091887/136‑1 (MQ=255)
                                                                                                                                                                              ctctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAg                                                                                                                                                                                                                                                      >  8:1019990/1‑118 (MQ=255)
                                                                                                                                                                              ctctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAg                                                                                                                                                                                                                                                      <  7:1019990/118‑1 (MQ=255)
                                                                                                                                                                                ctTTTAGTCTGCAAGGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTcg                                                                                                                                                                                                                      >  6:774553/1‑148 (MQ=255)
                                                                                                                                                                                ctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGc                                                                                                                                                                                                                                         >  7:214826/1‑129 (MQ=255)
                                                                                                                                                                                ctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGc                                                                                                                                                                                                                                         <  8:214826/129‑1 (MQ=255)
                                                                                                                                                                                ctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTcg                                                                                                                                                                                                                      >  8:734684/1‑148 (MQ=255)
                                                                                                                                                                                  tttAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCg                                                                                                                                                                                                                                                  <  5:553466/118‑1 (MQ=255)
                                                                                                                                                                                       tCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGt                                                                                                                                                                                                                                                               >  6:195160/1‑100 (MQ=255)
                                                                                                                                                                                          gCAACGCGGTGAGGCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGCCTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCaa                                                                             <  2:323960/275‑1 (MQ=255)
                                                                                                                                                                                          gCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGAtt                                                                                                                                                                                                            >  7:786644/1‑148 (MQ=255)
                                                                                                                                                                                            aaCGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                                                                                                                                                                          <  7:734684/148‑1 (MQ=255)
                                                                                                                                                                                             aCGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAg                                                                                                                                                                                                                                                           >  5:1086587/1‑98 (MQ=255)
                                                                                                                                                                                                  gTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAg                                                                                                                                                                                                    <  8:786644/148‑1 (MQ=255)
                                                                                                                                                                                                    gagaCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGttt                                                                                                                                                                                                                                                             >  8:36289/1‑89 (MQ=255)
                                                                                                                                                                                                    gagaCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGttt                                                                                                                                                                                                                                                             <  7:36289/89‑1 (MQ=255)
                                                                                                                                                                                                     agaCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGa                                                                                                                                                                                                                   <  6:729217/130‑1 (MQ=255)
                                                                                                                                                                                                     agaCGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGa                                                                                                                                                                                                                   >  5:729217/1‑130 (MQ=255)
                                                                                                                                                                                                         gCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAg                                                                                                                                                                                                                                                      <  8:152687/91‑1 (MQ=255)
                                                                                                                                                                                                         gCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAg                                                                                                                                                                                                                                                      >  7:152687/1‑91 (MQ=255)
                                                                                                                                                                                                              cAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATaa                                                                                                                                                                                                                                              >  8:829938/1‑94 (MQ=255)
                                                                                                                                                                                                              cAATCGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATaa                                                                                                                                                                                                                                              <  7:829938/94‑1 (MQ=255)
                                                                                                                                                                                                                aTTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCCTTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGa                                                                                                                                                                                                                   <  8:743923/119‑1 (MQ=255)
                                                                                                                                                                                                                aTTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGa                                                                                                                                                                                                                   >  7:743923/1‑119 (MQ=255)
                                                                                                                                                                                                                     tgGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGAt                                                                                                                                                                                                             <  6:100348/120‑1 (MQ=255)
                                                                                                                                                                                                                     tgGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGAt                                                                                                                                                                                                             >  5:100348/1‑120 (MQ=255)
                                                                                                                                                                                                                        cGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGc                                                                                                                                                                                                                                         <  5:457798/89‑1 (MQ=255)
                                                                                                                                                                                                                           aTCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCTCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGa                                                                                                                             <  4:29162/194‑1 (MQ=255)
                                                                                                                                                                                                                               ggCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGc                                                                                                                                                                                                   <  5:231296/120‑1 (MQ=255)
                                                                                                                                                                                                                                 cTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                                                                                                                                                                          >  8:573519/1‑111 (MQ=255)
                                                                                                                                                                                                                                 cTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                                                                                                                                                                          <  7:573519/111‑1 (MQ=255)
                                                                                                                                                                                                                                  tCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAg                                                                                                                                                                                                                                                      <  6:267940/66‑1 (MQ=255)
                                                                                                                                                                                                                                   ccGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGttt                                                                                                                                                                                                                                                             >  5:352762/1‑58 (MQ=255)
                                                                                                                                                                                                                                   ccGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGttt                                                                                                                                                                                                                                                             <  6:352762/58‑1 (MQ=255)
                                                                                                                                                                                                                                      gTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGCAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCAGGTTCGCGGTGCGGAAATGGCGATGATATTTCAGGAGCCGATGACATAGCTGAACCCTGTATTTACTGAGGGTGACACGATTGACGCATGACTTCGTCTGCATCAGAAGTCCAGTCGTGAAGAAGCGATGGta                                    >  1:155320/1‑271 (MQ=255)
                                                                                                                                                                                                                                      gTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTc                                    >  4:184256/1‑272 (MQ=255)
                                                                                                                                                                                                                                       tAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGc                                                                                                                                                                                                                                         >  7:1128411/1‑74 (MQ=255)
                                                                                                                                                                                                                                       tAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGc                                                                                                                                                                                                                                         <  8:1128411/74‑1 (MQ=255)
                                                                                                                                                                                                                                            cAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAg                                                                                                                                                                                                                                                      <  5:743372/56‑1 (MQ=255)
                                                                                                                                                                                                                                              gTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCg                         >  2:346050/1‑275 (MQ=255)
                                                                                                                                                                                                                                                 aCTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGAt                      <  2:154468/275‑1 (MQ=255)
                                                                                                                                                                                                                                                  cTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCg                                                                                                                                                                                                  <  8:165721/102‑1 (MQ=255)
                                                                                                                                                                                                                                                  cTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCg                                                                                                                                                                                                  >  7:165721/1‑102 (MQ=255)
                                                                                                                                                                                                                                                     cGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGc                    <  2:233243/273‑1 (MQ=255)
                                                                                                                                                                                                                                                     cGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTg                  >  2:421972/1‑275 (MQ=255)
                                                                                                                                                                                                                                                       ttGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGGTAAGCGGAATCTGGa                >  4:233242/1‑275 (MQ=255)
                                                                                                                                                                                                                                                       ttGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTgg                 >  2:19853/1‑274 (MQ=255)
                                                                                                                                                                                                                                                         ggCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGAAATCGCTGa                                                                                                                    >  3:145296/1‑173 (MQ=255)
                                                                                                                                                                                                                                                             ttGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                                                                                                                                                                          <  8:907265/83‑1 (MQ=255)
                                                                                                                                                                                                                                                             ttGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                                                                                                                                                                          >  7:907265/1‑83 (MQ=255)
                                                                                                                                                                                                                                                               gATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCg                         >  4:289632/1‑258 (MQ=255)
                                                                                                                                                                                                                                                               gATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAAATTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCATTCGTGAAAAAGCGAGGGTCGAGGCGAAGCGGGAGCTGGCGCAGGTa         >  4:472006/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                aTGCGCCTGTTGGAACAGGCGGGC‑GGTGTAGTACAGTGCGATAAAAGGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACg       <  2:543360/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cgcCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                                                                                                                                                                          >  8:1004380/1‑77 (MQ=255)
                                                                                                                                                                                                                                                                   cgcCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGa                                                                                                                                                                                                          <  7:1004380/77‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cgcCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATg                     <  2:289633/258‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cgcCTGTTCCAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGa                                                                                                                                                                                            >  8:127750/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                   cgcCTGTTCCAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGa                                                                                                                                                                                            >  6:50853/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                   cgcCTGTTCCAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGa                                                                                                                                                                                            <  7:127750/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ccTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGCAGGTCGAGGCGAAGCGGATGGCGGATCAGGTACGCATTc  >  4:362609/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                     ccTGTTCCAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGtt                                                                                                                                                                                                                                     >  6:689243/1‑48 (MQ=38)
                                                                                                                                                                                                                                                                       tGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTa                                                                                                                                                                                                     <  6:1085580/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGca                                                                                                                                                                                    <  5:1033705/92‑1 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                             
ATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGC‑GGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTC  >  NC_000913/867450‑867985

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑