Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 850,995 T→G 5.9% intergenic (+30/+19) ompX → / ← opgE outer membrane protein X/OPG biosynthetic transmembrane phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913850,9950TG 5.9% 282.0 / 10.0 94intergenic (+30/+19)ompX/opgEouter membrane protein X/OPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base T (40/46);  new base G (3/3);  total (45/49)
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.05e-01

GGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGAAATGACCATAAAAGGCACCTGGA  >  NC_000913/850745‑851258
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                       
gggCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCCCCTGGATTGCCGGTGTTGGTTACCGCTTCTAATACCTTTGGTGATATAAACAATCCCCCTCGCGGGGGCGATTTTTGTTTTTAAgt                                                                                                                                                                                                                                                 >  4:372243/1‑274 (MQ=255)
                tgtAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGGCGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAata                                                                                                                                                                                                                                 >  2:397921/1‑275 (MQ=255)
                 gtAGTGGGTGTGGGTTATGGTACCTTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTTTTACTATCAGATGCGTATTCGTTGCCTCGACCTTGGCCCCCGGATTTCCCGTGTTCGTTCCCGCTTCTAATCACTTTGTTTTACTCAACATTCCGCCTCACGGGGCGGACTTTTGTTTTTTAGGTTTGGTTTCGaaactta                                                                                                                                                                                                                                >  4:323660/1‑271 (MQ=255)
                  tAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAatat                                                                                                                                                                                                                                <  2:372248/274‑1 (MQ=255)
                   agtgtggtagtgtTATGGTTAATTCCANNCCACTGAGTACCGGACATACAAACGCGACANCATCGACTACGTTTTATCCTGCGGTGCGGGTCTTTAGGTCAACCCGATGGAACACGTTGCTCTTGACTTCTCTTNCGAGCAGAGACGTATTCGTAGCGTGGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGGGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTACGGTTTCGGGTCGAAAatat                                                                                                                                                                                                                                <  3:60454/262‑1 (MQ=255)
                        gtgtgggtTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTcc                                                                                                                                                                                                                         >  2:74971/1‑274 (MQ=255)
                          ttgggtTACGAAAAATTCCAACCCNNTGAATCACCGACCTACAAACAAGAAACCAGNGACTACGGTTTCGCCTGCGGGGCGGGTCGGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGANCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGt                                                                                                                                                                                                                            <  3:59481/269‑1 (MQ=255)
                          gtgggtTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAACACGTTGCTCTGGACTTCTCTTANGAGCATAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCANTTTGGTGATATAAAAAATCCGCCCCTCGGGGCGGATTTTTTTTTTNTNGTGTTCGGGTCGAATATATCGttttcta                                                                                                                                                                                                                       >  3:42797/1‑271 (MQ=255)
                             ggtTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCgaga                                                                                                                                                                                                                     <  2:53358/274‑1 (MQ=255)
                               ttATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCCGGATTTCCGGAGTTGGTTCCCTCTTCTAATCACTTTCGTGCGTTAAAAATTCCGCCGCCCGGGCGCGATTCTTTTTTTGAAGTGTTCGGCTCGAAAATTACGGTCCCggcattg                                                                                                                                                                                                                  >  4:148187/1‑269 (MQ=255)
                                       acttCCAGCCCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAAt                                                                                                                                                                                                           <  1:104301/272‑1 (MQ=255)
                                       aaTTCCAGACCACTGAATACCTGACCTACAAACACGACACCAGCGAATACCGTCTCTCCTACGGTGCGTGTCTGCAGTTCAACCCGTTGGAAAACGTTGCTCTGGACTTCTCTTACGATCAGAGCCGTATTCGTCGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTCGGTGATATACAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAAt                                                                                                                                                                                                           <  2:422569/274‑1 (MQ=255)
                                       aaTTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAg                                                                                                                                                                                                                                                  >  2:17671/1‑235 (MQ=255)
                                           ccAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGttt                                                                                                                                                                                                                                              <  4:17671/235‑1 (MQ=255)
                                              tcccACTGAATACCCCACCTACAAACCGGACTCCTGCGCTATCGGTTTTTCCTCCGGTGCGGGTCGGCAGTTCAACCCGATGGAAACGTTGTTCTGGACTTCTCTTACGGGCAGAGCCGTATTCGTAGCGTGGACGTAGGCACCTGGATTGCCGGGGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCt                                                                                                                                                                                                  <  4:74970/273‑1 (MQ=255)
                                               aCCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTCCCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTATCGGGGCTGATTTTTCTTTTTCAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCTACCTTCt                                                                                                                                                                                                  >  1:87307/1‑275 (MQ=255)
                                               aCCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCTCGTTGGTGATATAAAAAATCCGCCTCTCGGGGTGGATTTTTGTTTTTAAGGTTTTGGGTTGAAATANNNTGTTCCGGGGTGGTTGTTATTCAACttat                                                                                                                                                                                                  >  3:141361/1‑273 (MQ=255)
                                                 cACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGt                                                                                                                                                                                                <  1:42797/275‑1 (MQ=255)
                                                   cTGAATACCCGCCCTACAAACACGACACCAGCGACTCCGGTTTCTCCTACGTGTCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAGTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAAGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGt                                                                                                                                                                                                <  2:369495/273‑1 (MQ=255)
                                                   cTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTaa                                                                                                                                                                                              >  2:41621/1‑275 (MQ=255)
                                                   cTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTaa                                                                                                                                                                                              >  2:451558/1‑275 (MQ=255)
                                                           ccGACCTACAAACACGACTCCCGGTACTACGTATTCTCCTACGGTTCGTGTCTGCACTTCAACCCGATGGAAAACGTTGCTCTGTGCTTGTCTTACGAGCATGGCCGTATTCGTAGCGTTGACGTAGTCACCTGGATTGCCGGTGTTGGTTACGGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGt                                                                                                                                                                                      <  4:156665/275‑1 (MQ=255)
                                                           ccGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGTTCGAAAATATCGGTTCCGAGTTGGTTGTAATCCACCTTCTGTAACTGGaa                                                                                                                                                                                        >  2:428362/1‑273 (MQ=255)
                                                           ccGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCATCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGaaa                                                                                                                                                                                                                                    >  4:447542/1‑229 (MQ=255)
                                                             gaccgacgaaaccGACACCAGCGACTCCGGTTTCTCCGACGGGGGGGGTCTGGAGTTCAACCCGATTGAAAACGTTGCTCTGGACTTCTCTTACGGGCAACGCCGTATTCGTAGCGTTGACGTAGTCACCTGTATTGCCGGTGTTGGTTACCGCTTATATTCACTTTGGTGAGATAAAAAATCCGCCTCTCGGGGCGGATTTTTGGTTTTAAGGTTTCGGGTCGAAAAAATCGGTTCCGAAATGGTTGTAATCCACCTTCTGTAACTGGAAGTTg                                                                                                                                                                                    <  4:428357/263‑1 (MQ=255)
                                                              aCCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGTTTTTTGTTTTTTAGGTTTCGGGCCGAAAATATCGGTTCCGGGATGNNNGAAATCCACGTTCTGTAAGTGGAAGTTgt                                                                                                                                                                                   >  3:152934/1‑274 (MQ=255)
                                                               ccTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCATCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAatat                                                                                                                                                                                                                                <  2:447547/229‑1 (MQ=255)
                                                                                  cgccTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCt                                                                                                                                                               <  1:12543/272‑1 (MQ=255)
                                                                                    aCTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCttt                                                                                                                                                             >  1:109032/1‑275 (MQ=255)
                                                                                          gtttttTCCTACGGTGCGGGTATGCAGTTCAACCCGATGGAAAACGTTGCTCTGGAATTCTCTTACGTGCAGAGCCTNATTCGTAGCGTTGACGTAGGCACCTTGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTtcttct                                                                                                                                                        <  3:109032/273‑1 (MQ=255)
                                                                                          ggTTTCCCTTACGGTGCGGGTCTGCAGTGCAACCCGATTGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTGGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTGCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTtcttct                                                                                                                                                        <  4:41621/275‑1 (MQ=255)
                                                                                           gTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTtcttct                                                                                                                                                        <  1:152934/274‑1 (MQ=255)
                                                                                           gTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTtcttct                                                                                                                                                        <  2:323664/274‑1 (MQ=255)
                                                                                                         gCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCg                                                                                                                                                                                                                                                                       <  6:250125/149‑1 (MQ=255)
                                                                                                               cTGCAGTGCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCAACTGGATGGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGa                                                                                                                                   <  1:154438/275‑1 (MQ=255)
                                                                                                                   aGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttt                                                                                                                                                                                                                                                             <  6:140022/149‑1 (MQ=255)
                                                                                                                       cAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCggg                                                                                                                                                                                                                                                                     >  8:257258/1‑137 (MQ=255)
                                                                                                                       cAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCggg                                                                                                                                                                                                                                                                     <  7:257258/137‑1 (MQ=255)
                                                                                                                       cAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGt                                                                                                                                                                                                                                                         >  4:528456/1‑149 (MQ=255)
                                                                                                                         aCCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGttt                                                                                                                                                                                                                                                       <  5:454028/149‑1 (MQ=255)
                                                                                                                           ccGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGttttt                                                                                                                                                                                                                                                     <  2:528461/149‑1 (MQ=255)
                                                                                                                             gATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTaa                                                                                                                                                                                                                                                   <  6:346154/149‑1 (MQ=255)
                                                                                                                               tGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTCTTTGATGTTAATTTCCTTCg                                                                                                                    >  2:412538/1‑274 (MQ=255)
                                                                                                                                    aaCGTTGCTCTGGACTTCCCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTaa                                                                                                              <  2:218466/275‑1 (MQ=255)
                                                                                                                                       gTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGt                                                                                                                                                                                                                                         >  7:450891/1‑149 (MQ=255)
                                                                                                                                        ttGCTCTGGACTTCTTTTACGAGCAGAGCCGGATTCGTAGCGTTGGCGTAGGCACCTTGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATccc                                                                                                          <  1:141361/275‑1 (MQ=255)
                                                                                                                                              tGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAat                                                                                                                                                                                                                                  >  8:420939/1‑149 (MQ=255)
                                                                                                                                                gACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGACGTTGGTGATGTATATCGGCCCGGCTTTCTTCCCAGatt                                                                                                                                                  >  4:457955/1‑226 (MQ=255)
                                                                                                                                                   ttCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCgg                                                                                                                                                                                                                             >  6:460049/1‑149 (MQ=255)
                                                                                                                                                     ctctTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGgag                                                                                             >  2:217427/1‑275 (MQ=255)
                                                                                                                                                         tACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCga                                                                                                                                                                                                                       <  5:194395/149‑1 (MQ=255)
                                                                                                                                                            gAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGAt                                                                                                                                                                                                                    <  5:279816/149‑1 (MQ=255)
                                                                                                                                                            gAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttttta                                                                                                                                                                                                                                                         <  7:240376/112‑3 (MQ=38)
                                                                                                                                                            gAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttttta                                                                                                                                                                                                                                                         >  8:240376/1‑110 (MQ=38)
                                                                                                                                                                      gTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGa                                                                                                                                   <  1:47334/220‑1 (MQ=255)
                                                                                                                                                                            gTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGc                                                                      >  2:421172/1‑275 (MQ=255)
                                                                                                                                                                              aGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCg                                                                                                                                                                                                                        >  2:267217/1‑127 (MQ=255)
                                                                                                                                                                                  ttGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGAt                                                                                                                                                                                                                    <  4:267213/127‑1 (MQ=255)
                                                                                                                                                                                             aCCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGa                                                                                                                                                                                         <  6:364084/143‑1 (MQ=255)
                                                                                                                                                                                             aCCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGa                                                                                                                                                                                         >  5:364084/1‑143 (MQ=255)
                                                                                                                                                                                                       ccGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGAt                                                                                                                                                                                                                    >  2:532017/1‑106 (MQ=255)
                                                                                                                                                                                                        cGGTGTTTGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATc                                                                                                                                                                        >  6:260572/1‑149 (MQ=255)
                                                                                                                                                                                                           tgtTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGtt                                                                                                                                                                                                                <  4:532012/106‑1 (MQ=255)
                                                                                                                                                                                                                     cGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttct                                                                                                                                                           >  6:230995/1‑149 (MQ=255)
                                                                                                                                                                                                                         tCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTtctc                                                                                                                                                       <  5:230995/149‑1 (MQ=255)
                                                                                                                                                                                                                            aaTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGAt                                                                                                                                                                               >  8:374238/1‑122 (MQ=255)
                                                                                                                                                                                                                            aaTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGAt                                                                                                                                                                               <  7:374238/122‑1 (MQ=255)
                                                                                                                                                                                                                            aaTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCa                                                                                                                                                      <  8:9797/147‑1 (MQ=255)
                                                                                                                                                                                                                            aaTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCa                                                                                                                                                      >  7:9797/1‑147 (MQ=255)
                                                                                                                                                                                                                             aTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGt                                                                                                                                                                             <  6:140858/123‑1 (MQ=255)
                                                                                                                                                                                                                             aTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGt                                                                                                                                                                             >  5:140858/1‑123 (MQ=255)
                                                                                                                                                                                                                             aTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGat                                                                                                                                                   <  5:460049/149‑1 (MQ=255)
                                                                                                                                                                                                                                aCTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGAtata                                                                                                                                                >  7:472548/1‑149 (MQ=255)
                                                                                                                                                                                                                                  tttGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATaaa                                                                                                                                              <  7:420939/149‑1 (MQ=255)
                                                                                                                                                                                                                                      gTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACggg                                                                                                                                          <  5:260572/149‑1 (MQ=255)
                                                                                                                                                                                                                                       tGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTcc                                                                                                                                                                                                                         >  3:128216/1‑69 (MQ=255)
                                                                                                                                                                                                                                         atatAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACgg                                                                                                                                           >  6:197943/1‑145 (MQ=255)
                                                                                                                                                                                                                                         atatAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACgg                                                                                                                                           <  5:197943/145‑1 (MQ=255)
                                                                                                                                                                                                                                          tataAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTAtt                                                                                                                                      <  8:450891/149‑1 (MQ=255)
                                                                                                                                                                                                                                           ataAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCgaga                                                                                                                                                                                                                     <  1:128216/69‑1 (MQ=255)
                                                                                                                                                                                                                                            taAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCttt                                                                                                                                                             >  5:188833/1‑124 (MQ=255)
                                                                                                                                                                                                                                            taAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCttt                                                                                                                                                             <  6:188833/124‑1 (MQ=255)
                                                                                                                                                                                                                                            taAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttc                                                                                                                                                            >  2:545555/1‑125 (MQ=255)
                                                                                                                                                                                                                                            taAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTg                                                                                                                                    >  5:32489/1‑149 (MQ=255)
                                                                                                                                                                                                                                            acaaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAAc                                                                                                                                             <  6:186230/138‑1 (MQ=255)
                                                                                                                                                                                                                                            acaaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAAc                                                                                                                                             >  5:186230/3‑140 (MQ=255)
                                                                                                                                                                                                                                              aaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGaa                                                                                                                                                                                        <  5:56807/95‑1 (MQ=255)
                                                                                                                                                                                                                                              aaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGaa                                                                                                                                                                                        >  6:56807/1‑95 (MQ=255)
                                                                                                                                                                                                                                                caaTCCGCCTCTCGGTGCGGATTTTTGTTTTTAAGGTTTCGGGTCTAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTtcttct                                                                                                                                                        <  4:545550/124‑1 (MQ=255)
                                                                                                                                                                                                                                                aaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGACAAAGCCTAAGAAGTCAATGGCAGAGCGGCGGGCTTTAATCAAACGCTTCCCTTTATGTTCCCTGGAAATGACAAAAAAAGGCCCCGGGa  >  4:407488/1‑275 (MQ=255)
                                                                                                                                                                                                                                                 aaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGtt                                                                                                                               >  7:429797/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      gCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTa                                                                                                                                        <  8:131951/135‑1 (MQ=255)
                                                                                                                                                                                                                                                      gCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTa                                                                                                                                        >  7:131951/1‑135 (MQ=255)
                                                                                                                                                                                                                                                       ccTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTgg                                                                                                                                                                                   >  2:335730/1‑91 (MQ=255)
                                                                                                                                                                                                                                                         tctcGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCt                                                                                                                       >  8:225652/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                       
GGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGAAATGACCATAAAAGGCACCTGGA  >  NC_000913/850745‑851258

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