Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 850,993 T→C 5.6% intergenic (+28/+21) 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,9930TC 5.6% 293.3 / 6.1 94intergenic (+28/+21)ompX/opgEouter membrane protein X/OPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base T (41/46);  new base C (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 = 4.96e-02

GGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGAAATGACCATAAAAGGCACCTGGA  >  NC_000913/850720‑851258
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                         
ggTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCctc                                                                                                                                                                                                                                                                          >  2:371917/1‑275 (MQ=255)
                         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)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                         
GGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCTCCTACGGTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGAAATGACCATAAAAGGCACCTGGA  >  NC_000913/850720‑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: