Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,262,255 G→C 10.1% S196T (AGC→ACC)  dusA → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,262,2550GC10.1% 550.4 / 43.6 170S196T (AGC→ACC) dusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (79/71);  new base C (8/9);  total (89/81)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-01

Only 100 of 173 total aligned reads displayed.

ATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCG  >  NC_000913/4261990‑4262520
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                          
accAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCGGGTGACGGGGAAAGCGCGTATTGGCATCGCGGACCAGGCCAGCTATTAATTTCTCTGCGATTTCATCATTACCGTTTCCGGCAAAGGCGAGGGTGAGATGTTCATCATCCCAGCCCGTAAACCCTGGCTTTCGGGGTAAGGCCCaaaaa                                                                                                                                                                                                                                                                    >  4:255871/3‑273 (MQ=255)
aTCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGTGTCTTGGCATCCACGGCCAGGACAGCTATGAATTGCTATGCGATTTCATCAATACCGTTTCCGGCCAAGGCGAGTGTGAGATGTTCATCATCCACCCACGTAAAGCCTGGCTTTCGGTGTTAATCCCgaacg                                                                                                                                                                                                                                                                    >  4:320213/1‑272 (MQ=255)
aTCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaag                                                                                                                                                                                                                                                                    >  2:387851/1‑274 (MQ=255)
   aaCCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaagaaa                                                                                                                                                                                                                                                                 >  2:458279/1‑274 (MQ=255)
    aCCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGATTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGTTTAAGCCCGaaaaaaaa                                                                                                                                                                                                                                                                >  4:441424/1‑274 (MQ=255)
                    cccGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCcgccgc                                                                                                                                                                                                                                               <  2:206391/275‑1 (MQ=255)
                                            cATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGGTAAGCCCGAACGAACACCGGGAAATCCCGCCGCTCGATTATCCGCGTGTGTAGCaat                                                                                                                                                                                                                      >  2:461155/1‑274 (MQ=255)
                                                                                        aCTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGCTTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGTAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggctg                                                                                                                                                                           >  2:505916/1‑271 (MQ=255)
                                                                                           gCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTa                                                                                                                                                                         >  2:344613/1‑273 (MQ=255)
                                                                                              tGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCaa                                                                                                                                                                     <  4:344612/274‑1 (MQ=255)
                                                                                                       tGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGaa                                                                                                                                                           <  2:371350/275‑1 (MQ=255)
                                                                                                                        aTTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACACTGTCGATTAACGGTGGTATCAAGTCCGCGGAAGAGGCAAAAGCACAccc                                                                                                                                          >  4:525475/1‑274 (MQ=255)
                                                                                                                           tcggtgagggttaAAAAGCGGTTTTGCGTCGACGGAGAGGACAGGTATGAATTTTTCGGGGAGTTCTGTAATACGGTTACGGGCAAGGGGGGGGGTGAGATGTTCATCATCGACGGTCGTAAAGTCTGGCTTTCGGGGTTAAGCCCGGAAGAAAACCGTGAAATCCCGCAGCTCGATTATCCCCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                                                                                                        <  4:79306/273‑1 (MQ=255)
                                                                                                                           ccGGTGACGGGGAAAACGGGTATTGGAATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATAAATCCCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGAACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                                                                                                        <  4:282568/274‑1 (MQ=255)
                                                                                                                           ccGGTGACGGGGAAAACGCGTATTGGCATCGACGACCAGTACAGCTAGGAATGGCTCTTCGATTTCAGCAATGCCGGTTCGTGCAAATGCGAGTGTGAGATGTTCAGCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                                                                                                        <  4:120581/274‑1 (MQ=255)
                                                                                                                            cggtgacggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaa                                                                                                                                                                                                                                                                      <  5:653262/148‑1 (MQ=255)
                                                                                                                             gggtgcggtgaAAGCGCGTTATGGAATCGACAACCAGAACGGCTATGAATTTCTCTGCGGTTTCAACAATGCCGTTACAGGCAAAGGCGGGTGGGAGATGTTCATCAGCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCACCGCTGGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCaa                                                                                                                                      <  4:448198/269‑1 (MQ=255)
                                                                                                                                 acggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaagaaa                                                                                                                                                                                                                                                                 <  7:1045122/148‑1 (MQ=255)
                                                                                                                                      gaAAACGCTGATTGGCATCGTCGACCAGAAAAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGACAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatggat                                                                                                                             <  2:133648/274‑1 (MQ=255)
                                                                                                                                      gaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaa                                                                                                                                                                                                                                                                      >  7:721477/1‑138 (MQ=255)
                                                                                                                                       aaaaCGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCg                                                                                                                                                                                                                                                             <  4:399684/146‑1 (MQ=255)
                                                                                                                                            gcgTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                                                                                                                                                                                                      <  6:242414/148‑1 (MQ=255)
                                                                                                                                             cgTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGcc                                                                                                                                                     >  2:202968/1‑243 (MQ=255)
                                                                                                                                                 ttGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAg                                                                                                                                                 <  4:202967/243‑1 (MQ=255)
                                                                                                                                                      aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaa                                                                                                                                                                                                                                                                      >  5:355616/1‑122 (MQ=255)
                                                                                                                                                      aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaa                                                                                                                                                                                                                                                                      <  6:355616/122‑1 (MQ=255)
                                                                                                                                                      aTCGACGACAAGGCCAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGgcg                                                                                                            <  2:445838/275‑1 (MQ=255)
                                                                                                                                                        cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGcgcg                                                                                                          <  2:255872/275‑1 (MQ=255)
                                                                                                                                                                 ggACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaaga                                                                                                                                                                                                                                                                   >  7:206456/1‑114 (MQ=255)
                                                                                                                                                                 ggACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaaga                                                                                                                                                                                                                                                                   <  8:206456/114‑1 (MQ=255)
                                                                                                                                                                  gACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAg                                                                                                                                                                    >  1:71196/1‑207 (MQ=255)
                                                                                                                                                                     aGCTATGAATTTTTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGaa                                                                                             <  2:497998/275‑1 (MQ=255)
                                                                                                                                                                        tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                                                                                                                                                                                                                            <  7:285672/130‑1 (MQ=255)
                                                                                                                                                                          tGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgcg                                                                                                                                                                                                                                  <  6:76151/138‑1 (MQ=255)
                                                                                                                                                                          tGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgcg                                                                                                                                                                                                                                  >  5:76151/1‑138 (MQ=255)
                                                                                                                                                                           gAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCga                                                                                                                                                                                                                                                                       <  7:1003564/100‑1 (MQ=255)
                                                                                                                                                                           gAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCga                                                                                                                                                                                                                                                                       >  8:1003564/1‑100 (MQ=255)
                                                                                                                                                                                 ctctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTa                                                                                                                                                                         >  2:555152/1‑187 (MQ=255)
                                                                                                                                                                                    tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCAAACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTGTCAGAATCCGGGTATTCTggc                                                                              >  1:95952/1‑275 (MQ=255)
                                                                                                                                                                                     gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGa                                                                                                                                                                                                            >  5:334677/1‑148 (MQ=255)
                                                                                                                                                                                        aTTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgccg                                                                                                                                                                                                                                                <  8:864928/110‑1 (MQ=255)
                                                                                                                                                                                          ttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaa                                                                                                                                                                                                                                                                      <  6:118560/86‑1 (MQ=255)
                                                                                                                                                                                             atcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaa                                                                                                                                                                                                                                                                      <  5:409640/83‑1 (MQ=255)
                                                                                                                                                                                             atcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaa                                                                                                                                                                                                                                                                      >  6:409640/1‑83 (MQ=255)
                                                                                                                                                                                              tcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgtg                                                                                                                                                                                                                             >  6:170926/1‑122 (MQ=255)
                                                                                                                                                                                               caAGACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCAt                                                                                                                                                                                                  >  8:553486/1‑148 (MQ=255)
                                                                                                                                                                                                   aactacaccGGTAAAGGTGGGTGGAACACGTTCGTCAGCGACGTATGTAAAGCCTGGCTTTGGGGGTTAATCCCGAAAGTAAACCGTGAAATCCCCGATCTCGAATTTCGGTGGGGGTAGCAAATGGAGCGGTACTTTCCGCATCTTACAAGGTCCATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCAAACCTGCGACATAGGGATGGCGTGATGGTCGGGCGGGAGGCGTATCAGAATCCGGGGATCCTTGCGGCGGTAGACCGGga                                                               <  4:534349/268‑1 (MQ=255)
                                                                                                                                                                                                       tttCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCTTGACTTTCCGCATCTTACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCCTGATGTCCGGGCGCGAGGCGTCTCAGAATCCGGGTATTCTGGCGGCGGTAGTCCgggtgc                                                             >  1:109223/1‑270 (MQ=255)
                                                                                                                                                                                                         tCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTAACAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTGGGGCGAGAGGCGTATACGAATCCGGGTATGAGGGCGGCGGTGGACCg                                                                  >  1:8798/1‑266 (MQ=255)
                                                                                                                                                                                                            ggCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa                                                                                                                                                                                     <  6:334677/148‑1 (MQ=255)
                                                                                                                                                                                                            ggCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTgg                                                      >  1:90667/1‑275 (MQ=255)
                                                                                                                                                                                                             gCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGc                                                                                                                                                                                                                <  5:1045826/120‑1 (MQ=255)
                                                                                                                                                                                                             gCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGc                                                                                                                                                                                                                >  6:1045826/1‑120 (MQ=255)
                                                                                                                                                                                                                    cGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggt                                                                                                                                                                             <  7:553486/148‑1 (MQ=255)
                                                                                                                                                                                                                      aGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTACCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACg                                                                                                                                                                               >  8:359264/1‑144 (MQ=255)
                                                                                                                                                                                                                          tgAGCTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATACCGCCGCTCGATCACCCGCGAGTGTAACAACTGAAGCGTGACTTTCCTCACCTGACAATGTCGCTTAACGTCGGTGTCAAGTCGCTGGAAGAGGCCACAGCACACCGGACACATATTGTCGGGGTGAGGGTCGGGCGCGAGCCGTTACAGACCGCGGGTAATCTGGAGGAGGTAGACCGGGGGATCGTGGGTTCATCGGacagcg                                        >  4:204118/1‑270 (MQ=255)
                                                                                                                                                                                                                             gaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                                                                                                                                                                                               <  8:1049943/89‑1 (MQ=255)
                                                                                                                                                                                                                                   catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggt                                                                                                                                                                             >  6:712638/1‑133 (MQ=255)
                                                                                                                                                                                                                                   catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggt                                                                                                                                                                             <  5:712638/133‑1 (MQ=255)
                                                                                                                                                                                                                                   catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGATCGAATATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTCTTCTGGCGGGTGTAGACCGGGAGTACTTTGGTTTCTCGGTTAcggtagacggta                               >  4:432250/1‑264 (MQ=255)
                                                                                                                                                                                                                                    atcatcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGt                                                                                                                                                                                        >  6:571212/1‑121 (MQ=255)
                                                                                                                                                                                                                                    atcatcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGt                                                                                                                                                                                        <  5:571212/121‑1 (MQ=255)
                                                                                                                                                                                                                                    atcatcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTgg                                                                                                                                                             >  7:951389/1‑148 (MQ=255)
                                                                                                                                                                                                                                        tcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAaga                                                                                                                                                         >  7:1011578/1‑148 (MQ=255)
                                                                                                                                                                                                                                        tcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAaga                                                                                                                                                         >  6:828711/1‑148 (MQ=255)
                                                                                                                                                                                                                                            cgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTg                                                                                                                                                                                                             <  7:934007/92‑1 (MQ=255)
                                                                                                                                                                                                                                            cgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCa                                                                                                                                                                      >  7:693061/1‑131 (MQ=255)
                                                                                                                                                                                                                                            cgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCa                                                                                                                                                                      <  8:693061/131‑1 (MQ=255)
                                                                                                                                                                                                                                            cgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCt                                                                                                                                                               <  5:190923/138‑1 (MQ=255)
                                                                                                                                                                                                                                               acgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCaaa                                                                                                                                                  <  8:951389/148‑1 (MQ=255)
                                                                                                                                                                                                                                                 gTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCAGCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGATCCGGTGGCGGTAGTg                 <  2:525476/275‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTg                                                                                                                                                                                                             <  7:1068668/86‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt                                                                                                                                                                                                         <  5:429132/90‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACtt                                                                                                                                                                                                         >  6:429132/1‑90 (MQ=255)
                                                                                                                                                                                                                                                  tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTc                                                                                                                                                                                                       >  6:119746/1‑92 (MQ=255)
                                                                                                                                                                                                                                                  tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTc                                                                                                                                                                                                       <  5:119746/92‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAaga                                                                                                                                                         >  6:971889/1‑138 (MQ=255)
                                                                                                                                                                                                                                                  tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAaga                                                                                                                                                         <  5:971889/138‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaaGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCaa                                                                                                                                                   >  5:1088648/1‑143 (MQ=255)
                                                                                                                                                                                                                                                   aaaGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTAGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGATCCGGTGGCGGTAGTgcg               <  2:42536/275‑1 (MQ=255)
                                                                                                                                                                                                                                                         tGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                                                                                                        >  7:996521/1‑148 (MQ=255)
                                                                                                                                                                                                                                                         tGGCTTTCGGGGAAAAGCCAGAAAGAAAACAGTGAAAACCCGCCGCTAGAGTGTCCGCGTAGTGATCAACGGAAGCGTGACTTTCCGCATCTGAGAATGTCGAATAACGGTGGTATCAAGTCGCGGGAAGAGGCCAAAGCACACCTGCAACATAGGGATGGCGGGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGGATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGATCCGGTGGCGGTAGTGCGCGCCa          <  4:384392/274‑1 (MQ=255)
                                                                                                                                                                                                                                                          ttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCa                                                                                                                                       >  8:70824/3‑148 (MQ=255)
                                                                                                                                                                                                                                                          gtcTTGCGGGGTTAAGCCACAAAAAAACCCGGGAAACCCGGCCGCTAGATTTGCCGCGTTTGTAGCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATAAACGGTGGGATCAAGTCGCTGGAGGAGGCCAAAGCACACCGGCAGCATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGATCCGGTGGCGGTAGTGCGCGCCATg        <  4:458278/273‑1 (MQ=255)
                                                                                                                                                                                                                                                          ggCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCa                                                                                                                                       >  5:313683/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           gCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCacac                                                                                                                                            >  1:111704/1‑142 (MQ=255)
                                                                                                                                                                                                                                                              ttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAAcgt                                                                                                                                   <  8:1011578/148‑3 (MQ=255)
                                                                                                                                                                                                                                                              ttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACat                                                                                                                                   >  5:836124/1‑148 (MQ=255)
                                                                                                                                                                                                                                                              ttCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                                                                                                                                                                                                                          >  7:1072604/1‑46 (MQ=37)
                                                                                                                                                                                                                                                              ttCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                                                                                                                                                                                                                          <  8:1072604/46‑1 (MQ=37)
                                                                                                                                                                                                                                                              ttCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGa                                                                                                                                                            >  7:4638/1‑123 (MQ=255)
                                                                                                                                                                                                                                                               tCGNGGTTAAGCCCGAAAGAAAACCGTGACATCCCGCCGCTAGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTGCCGCATCTGACAATGTCGATTAAAGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                                                                                                        <  3:111704/142‑1 (MQ=255)
                                                                                                                                                                                                                                                                cGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACatat                                                                                                                                 <  5:828711/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tggtttAAGCCCGGAATAAAACCGTGAACTCCCGCCGCTCGATTATCCGCGTGGGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCg  <  2:87565/270‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ggggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCAGGAGATCTTTGGTTCCTCGGATACCGATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCg  <  2:290795/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGca                                                                                                                                               <  8:268244/128‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATc                                                                                                                                                                                                 <  6:513627/76‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATc                                                                                                                                                                                                 >  5:513627/1‑76 (MQ=255)
                                                                                                                                                                                                                                                                         gCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggtg                                                                                                                                                                            >  6:472338/1‑96 (MQ=255)
                                                                                                                                                                                                                                                                         gCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggtg                                                                                                                                                                            <  5:472338/96‑1 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                          
ATCAACCTGAATGTCGGCTGCCCGTCTGACCGGGTGCAGAACGGCATGTTTGGTGCGTGTCTGATGGGTAATGCGCAGCTGGTTGCCGACTGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGATCCGGTGGCGGTAGTGCGCGCCATGTATCCG  >  NC_000913/4261990‑4262520

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