Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 747,900 G→C 5.6% intergenic (‑131/+21) abrB ← / ← ybgO regulator of aidB expression; inner membrane protein/putative fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913747,9000GC 5.6% 432.5 / 12.0 124intergenic (‑131/+21)abrB/ybgOregulator of aidB expression; inner membrane protein/putative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (52/65);  new base C (3/4);  total (55/69)
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.89e-01

Only 100 of 125 total aligned reads displayed.

AGTAATTCCGCGCATACTAAAGATGATCCCGGCGATCATCGGCCCAAGTAATAGCGCCGCCGGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGTAGCTTCAGACCATTACCAATCAACAATGCTTTATCTTCTTCAACCAACGTTTCTTCGGTATAAAAGGAG  >  NC_000913/747632‑748173
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                 
          tgcATACTAAAGATGATCCCGGCGATCATCGGCACAAGTAATAGCGCCGCCGGGAAATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCATCATTATTTTCTCTTTGGTGTATCGATCCAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAAGTTTACGAAAGAGATGCCCTTTACGGCGGCAAAGGGCAACATgg                                                                                                                                                                                                                                                                   >  4:407324/2‑275 (MQ=255)
               aCTAAAGATGATCCCGGCGATCATAGGCCCAAGGAAAAGCGCCGCCGGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAAt                                                                                                                                                                                                                                                              <  4:40842/275‑1 (MQ=255)
                 tcaaGATGATCCCGGCGATCATCGGCCCAATTAATAGCGCCGCCGGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTAATAGCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAAt                                                                                                                                                                                                                                                              <  2:434920/271‑1 (MQ=255)
                                     accGGCCCAAGTAAAAGGGCCGCCGGGGGATGTACCGGGAGGAAACCAAAAGAAAGGAGAAGGGATAGCACACATAACATTCCCCACTGCAAAACTGGCATCCGTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCTAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaa                                                                                                                                                                                                                                          <  4:515945/271‑1 (MQ=255)
                                      taggCCCAAGAAATAGGGCAGCCGGGAGATGTACCGGGAGGAAACCAATAGAAAGGAGACTTGATAGCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCCCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAg                                                                                                                                                                                                                                       <  4:109547/273‑1 (MQ=255)
                                         gaccAAATAATAGCGCCCGCGGGAGAATGACCGCGAGGAAACCAATAGAAAGCGGAAGTGATAGCAAACCTAACATTNCCCACTGCAAAACTGGCATCCTTCGCTCCCACATTATTTTCTCTTTGGTGTATCGATAAAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCACTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                                                                                                                                                                    <  3:132459/273‑1 (MQ=255)
                                               gTAATAGCGCCGCCGGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCCATCAAGGAAAAGAGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATCGTTAATCTTAAAAACCACCACTg                                                                                                                                                                                                                               >  4:140667/1‑274 (MQ=255)
                                                        ccgccgGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                                                                                                                                                                    >  4:419008/1‑260 (MQ=255)
                                                            cgGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACt                                                                                                                                                                                                                                <  2:419009/260‑1 (MQ=255)
                                                                         gcgAGGAAACCAATTGAAAGGAGAAGTGATAGTACACATAACATTCCCCACTGGAAAGCTGGCATCCTTCGCTCCCTCATTATGTTCTCGTTGGTGGATCGATACAGGATAAAGAAAGATGGGCCAGGAAAAAAGAACAGCAAGCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGGGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                                                                                                                                                                     <  4:216391/274‑1 (MQ=255)
                                                                                       aGAAAGGAGAAGTGATAGCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGAc                                                                                                                                                                                      <  2:187409/275‑1 (MQ=255)
                                                                                             gagaAGTGATAGCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTTTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATTCCCTTTCCGGCGGCAAAGGGCAACATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGAGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGt                                                                                                                                                                                >  2:234123/1‑275 (MQ=255)
                                                                                                 aGTGATAGCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAg                                                                                                                                                                            <  2:480373/275‑1 (MQ=255)
                                                                                                       aGCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAgg                                                                                                                                                                                                                                                                            >  4:148389/1‑173 (MQ=255)
                                                                                                        gCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtat                                                                                                                                                                      >  2:141827/1‑274 (MQ=255)
                                                                                                           cacaTAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcat                                                                                                                                                                                                                                                                        <  2:148390/173‑1 (MQ=255)
                                                                                                            acaTAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGc                                                                                                                                                                                                                   >  4:368499/1‑225 (MQ=255)
                                                                                                             caTAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAgg                                                                                                                                                                                                                                                                            >  4:401616/1‑167 (MQ=255)
                                                                                                                aaCATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGc                                                                                                                                                                                                               <  2:368500/225‑1 (MQ=255)
                                                                                                                 aCATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcat                                                                                                                                                                                                                                                                        <  2:401617/167‑1 (MQ=255)
                                                                                                                       cccACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGAc                                                                                                                                                      >  4:77599/1‑275 (MQ=255)
                                                                                                                       cccACTGCAAAAATGGCAACCTTCGCTCCCTCATTATTTTCGCTTTGGTGTATAGATACAGCATAAAGAAAGCTTGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCGCTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGAc                                                                                                                                                      <  4:141826/275‑1 (MQ=255)
                                                                                                                            tGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                                                                                                                                                        >  4:265893/1‑236 (MQ=255)
                                                                                                                            tGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGATAATCTTACAAAGCACCACGGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCAGGAGAGGTTTTTTCGTAGACCaaaa                                                                                                                                                 >  4:508764/1‑273 (MQ=255)
                                                                                                                                aaaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAgg                                                                                                                                                                                                                                                                            >  8:799064/1‑148 (MQ=255)
                                                                                                                                aaaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGa                                                                                                                                                                                    <  2:265894/236‑1 (MQ=255)
                                                                                                                                aaaCTGGCATCCTTCGCTCCCTCATTATTTTCGCTTTGGTGGATCGATACAGCATAAAGAAAACTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAAGTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGAGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGAc                                                                                                                                                      <  2:508765/266‑1 (MQ=255)
                                                                                                                                   cTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAggg                                                                                                                                                                                                                     >  3:138925/1‑200 (MQ=255)
                                                                                                                                       cATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCa                                                                                                                                                                                                                 <  1:138925/200‑1 (MQ=255)
                                                                                                                                        aTCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATCGGTTTTTTCGTAACCCAACATGTCCATGGTGa                                                                                                                                     >  4:389838/1‑275 (MQ=255)
                                                                                                                                              cGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATATAGAAAGCTGGGCCAGTAAAACATAACAGCAAACAAGGAAAAGGGGCAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCCAAGGGGCCGGCCTCGATGGGGTTGCCTGGATTGGACGAAAGTTCAGCCGTATAGGTTTTTTGGTAGACCAACATGTCCCGGGTTAACTCga                                                                                                                               >  4:51190/1‑274 (MQ=255)
                                                                                                                                                 tCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcacca                                                                                                                                                                                                                                  <  4:466122/173‑1 (MQ=255)
                                                                                                                                                      cATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTAt                                                                                                                                                                                                                            >  4:317482/1‑174 (MQ=255)
                                                                                                                                                        ttattTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGc                                                                                                                     <  2:51191/275‑1 (MQ=255)
                                                                                                                                                          attTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                                                                                                                                                        <  2:317483/174‑1 (MQ=255)
                                                                                                                                                                  tttGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAgg                                                                                                                                                                                                                                                                            >  4:77675/1‑114 (MQ=255)
                                                                                                                                                                  tttGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcatc                                                                                                                                                                                                                                                                       >  4:569774/1‑119 (MQ=255)
                                                                                                                                                                  ttgtgtgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTg                                                                                                           <  4:370454/271‑1 (MQ=255)
                                                                                                                                                                   ttGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGc                                                                                                                                                                                                                                      >  2:547839/1‑151 (MQ=255)
                                                                                                                                                                      gtgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcat                                                                                                                                                                                                                                                                        <  2:77676/114‑1 (MQ=255)
                                                                                                                                                                      gtgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATgg                                                                                                                                                                                                                                                                   <  2:569775/119‑1 (MQ=255)
                                                                                                                                                                       tgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcacca                                                                                                                                                                                                                                  <  4:547838/151‑1 (MQ=255)
                                                                                                                                                                         tATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcatca                                                                                                                                                                                                                                                                      >  2:168153/1‑113 (MQ=255)
                                                                                                                                                                              aTACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAg                                                                                                                                                                                                                                                                             >  5:1066682/1‑101 (MQ=255)
                                                                                                                                                                                 aaaaaTAAAAAAAGCTGGACCAAGAAAAAAAAACAGCAATCAAGGCAAAGGGGAAAATCAGCAATTTTCGTAAAGAGNTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAaggcgag                                                                                            <  3:117321/271‑1 (MQ=255)
                                                                                                                                                                                   gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGa                                                                                                                                                                                                                         <  7:799064/148‑1 (MQ=255)
                                                                                                                                                                                    cATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCaa                                                                                                                                                                                                                                                                               <  7:728137/93‑1 (MQ=255)
                                                                                                                                                                                    cATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCAtattta                                                                                                              >  4:504998/1‑254 (MQ=255)
                                                                                                                                                                                        aagaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGa                                                                                                          <  2:504999/254‑1 (MQ=255)
                                                                                                                                                                                           aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATg                                                                                                                                                                                                                                                                    >  5:509413/1‑97 (MQ=255)
                                                                                                                                                                                           aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                                                                                                                                                                         <  5:186716/124‑1 (MQ=255)
                                                                                                                                                                                           aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                                                                                                                                                                         >  6:186716/1‑124 (MQ=255)
                                                                                                                                                                                             agCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTc                                                                                                                                             <  1:62746/214‑1 (MQ=255)
                                                                                                                                                                                                tGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGt                                                                             <  2:77600/275‑1 (MQ=255)
                                                                                                                                                                                                  ggCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCg                                                                                                                                                                         >  3:121903/1‑181 (MQ=255)
                                                                                                                                                                                                      aGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtata                                                                                                                                                                     <  1:121903/181‑1 (MQ=255)
                                                                                                                                                                                                      aGGAAAAAATAACAGCAAGCCATGAAAAGGGGTAAAGCAGCAATTGTCGGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGt                                                                       <  2:389839/275‑1 (MQ=255)
                                                                                                                                                                                                         aaaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGc                                                                                                                                                                                                                                                                          <  8:679965/77‑1 (MQ=37)
                                                                                                                                                                                                         aaaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAggg                                                                                                                                                                                                                                                                           >  5:735154/1‑76 (MQ=255)
                                                                                                                                                                                                         aaaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAggg                                                                                                                                                                                                                                                                           <  6:735154/76‑1 (MQ=255)
                                                                                                                                                                                                           aaaaaaACAGCAAGCAAGGGAAGGGGGAAAATCAGCAATTTTCAGAAAGGGATGCCCTTTCCGGCGGCAAAGGGCAGCATGGTAAATCAGTTATAGTTAAACTTAAAGAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCAGCCAGTAAACGTAGCtt                                                                  <  1:116822/275‑1 (MQ=255)
                                                                                                                                                                                                           aaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATAGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGTAGCtt                                                                  <  2:428611/275‑1 (MQ=255)
                                                                                                                                                                                                           aaaaTAACAGCAATCAAGGAAAAGGGGAAAAACAGCAATTGTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGAATCATGGTAAACCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGGGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGTAGCtt                                                                  <  2:215276/275‑1 (MQ=255)
                                                                                                                                                                                                                 aCAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgg                                                                                                                                                                                                       >  8:810613/1‑136 (MQ=255)
                                                                                                                                                                                                                 aCAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgg                                                                                                                                                                                                       <  7:810613/136‑1 (MQ=255)
                                                                                                                                                                                                                    gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATc                                                                                                                                                                                                                                                             <  6:991130/79‑1 (MQ=255)
                                                                                                                                                                                                                    gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATc                                                                                                                                                                                                                                                             >  5:991130/1‑79 (MQ=255)
                                                                                                                                                                                                                    gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                                                                                                                                                                          >  5:306382/1‑82 (MQ=255)
                                                                                                                                                                                                                    gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                                                                                                                                                                          <  6:306382/82‑1 (MQ=255)
                                                                                                                                                                                                                    gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                                                                                                                                                                      <  5:932131/134‑1 (MQ=255)
                                                                                                                                                                                                                    gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                                                                                                                                                                      >  6:932131/1‑134 (MQ=255)
                                                                                                                                                                                                                      aaTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcaccac                                                                                                                                                                                                                                 <  6:461465/105‑1 (MQ=255)
                                                                                                                                                                                                                        tCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                                                                                                                                                        >  8:109085/1‑112 (MQ=255)
                                                                                                                                                                                                                        tCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                                                                                                                                                        <  7:109085/112‑1 (MQ=255)
                                                                                                                                                                                                                         cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATgg                                                                                                                                                                                                                                                                   <  6:862312/68‑1 (MQ=255)
                                                                                                                                                                                                                           aGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATg                                                                                                                                                                                                        >  7:436919/1‑125 (MQ=255)
                                                                                                                                                                                                                           aGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATg                                                                                                                                                                                                        <  8:436919/125‑1 (MQ=255)
                                                                                                                                                                                                                             gAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGa                                                                                                                                                                                                          >  5:912091/1‑121 (MQ=255)
                                                                                                                                                                                                                               aaaGGGGAAAATCAGCAATTTTCTGAAAGTGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGTAGCTTCAGACCATTACCAATcaaca                                              <  2:407325/275‑1 (MQ=255)
                                                                                                                                                                                                                                 aGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCt                                                                                                                                                                                               >  5:538301/1‑128 (MQ=255)
                                                                                                                                                                                                                                 aGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCt                                                                                                                                                                                               <  6:538301/128‑1 (MQ=255)
                                                                                                                                                                                                                                      aaaaTCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAg                                                                                                                                                                                                                       <  8:569128/99‑1 (MQ=255)
                                                                                                                                                                                                                                      aaaaTCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAg                                                                                                                                                                                                                       >  7:569128/1‑99 (MQ=255)
                                                                                                                                                                                                                                             gCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAg                                                                                                                                                                                                                       <  8:790658/92‑1 (MQ=255)
                                                                                                                                                                                                                                             gCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAg                                                                                                                                                                                                                       >  7:790658/1‑92 (MQ=255)
                                                                                                                                                                                                                                                            agaTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAg                                                                                                                                                                                                                                       <  8:149913/61‑1 (MQ=255)
                                                                                                                                                                                                                                                              aTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCGTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGTAGCTTCAGACCATTACCAATCAACAATGCTTTATCTTCTTCAACCAACGTttcttc               <  2:40250/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                gCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGa                                                                                                                                                                                                          <  8:836753/86‑1 (MQ=38)
                                                                                                                                                                                                                                                                 cccTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                                                                                                                                                                      <  8:212358/89‑1 (MQ=39)
                                                                                                                                                                                                                                                                 cccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAACGCACCACTGTATCGAAGGGGCCAGCCTCGATggggg                                                                                                                                                                                                    >  5:984104/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                 cccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggggg                                                                                                                                                                                                    <  6:984104/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                                                                                                                                                        >  5:640374/1‑103 (MQ=255)
                                                                                                                                                                                                                                                                 cccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                                                                                                                                                        <  6:640374/103‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtat                                                                                                                                                                      >  7:342309/1‑120 (MQ=255)
                                                                                                                                                                                                                                                                  ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                                                                                                                                                                      >  7:63591/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                  ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                                                                                                                                                                      <  8:63591/88‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaa                                                                                                                                                                                   <  5:288351/107‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTAGCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaa                                                                                                                                                                                   >  6:288351/1‑107 (MQ=255)
                                                                                                                                                                                                                                                                   cTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACAt                                                                                                                                                <  7:943323/141‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cggcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATATCACTGGTGAAACCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGTAGCTTCAGACCATTACCAAACAACAATGCTTTATCTTCTTTAACCAACGTTTTTGTGGTATAAAagggg  >  4:21591/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                 
AGTAATTCCGCGCATACTAAAGATGATCCCGGCGATCATCGGCCCAAGTAATAGCGCCGCCGGGAGATGTACCGCGAGGAAACCAATAGAAAGGAGAAGTGATAGCACACATAACATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGTATTCAGCATATTTATTGAAGGTATAAGGCGAGGCGTTTTCATCCAGTAAACGTAGCTTCAGACCATTACCAATCAACAATGCTTTATCTTCTTCAACCAACGTTTCTTCGGTATAAAAGGAG  >  NC_000913/747632‑748173

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