Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,631,245 C→T S171S (AGC→AGT)  slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,631,2450CT97.5% 230.2 / ‑2.3 83S171S (AGC→AGT) sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (41/38);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

CGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGC  >  NC_000913/4631000‑4631490
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
cgcgCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTg                                                                                                                                                                                                                                                    <  1:29728/251‑1 (MQ=255)
  cgcTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTggg                                                                                                                                                                                                                                                  >  1:29197/1‑250 (MQ=255)
               cTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTa                                                                                                                                                                                                                                      >  1:428555/1‑250 (MQ=255)
                           cGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGGGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCt                                                                                                                                                                                                                         >  1:129965/1‑251 (MQ=255)
                             cACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGCACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGCACACCGGGCCGAGTTAAGAAGCCTGGCAAGGGGCGAAAGAACTGTGGCTAACCGGC‑CAGAGCCAGCCTAACGCCTGTGACCAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTgg                                                                                                                                                                                                                       >  1:415633/1‑251 (MQ=255)
                                  tGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGGCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAAATACTACTATGCAAAATGGAACACCGGGCAAAGAGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGAGGCTAACCGGCCAAG‑GCCAGCCTAACGCCTTGGACAAGTTTTTTAGTGTCTGGCGTGGGGCCGGTAAAAAAAATACGCTGGCGTAt                                                                                                                                                                                                                  >  2:184589/1‑251 (MQ=255)
                                    aaagccagttgcGTCAAGGAACCGGCGCGGCGTGAAGACTGGCTGGGCTGGGTCGCCTTTAGCCCGTAAAAGCCCGGAACTACCGAAGCGCAATTGTA‑TTCCTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑ACGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGAGTGCGTCAGGTAAACAAGATCCGCttgcttagtt                                                                                                                                                                                                                <  2:223674/240‑10 (MQ=255)
                                      aTCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACCAGATCCGCTGGCGTAttt                                                                                                                                                                                                                >  1:224920/1‑249 (MQ=255)
                                           gTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGAAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGGGTGAAAAAGGCTGGCAAGGGGAGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGGGACAAGTTATTTAGTGTCTGGAGTGCGTCAGGTCAACAAGATCCGCTGGCGTGTTTAAAGCGt                                                                                                                                                                                                         >  1:403164/1‑251 (MQ=255)
                                                 tCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGATAGCCTTTTGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCCAAATGGCACACCGGGCAGAGTGAAGAAGCCAGGCAAGGGGCGAAAGAGATGTGGCTAACCGGC‑AAGAGCCCGCCTAACGCCGGTGACCAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGt                                                                                                                                                                                                   >  1:414283/1‑251 (MQ=255)
                                                         cTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGa                                                                                                                                                                                            >  1:484003/1‑250 (MQ=255)
                                                         cTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAACCCGGGACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGG‑CGAAAGAGCTGTGGCTAACCGGC‑AAGAGGCCGCCTAACGCCTGTGACCAGTTATTTTGTGTCTGGCGTGCGTCAGGGAAAAAAGATACGCTTGCGTATTTAGAGCGTGTCCGTCTTGCGATg                                                                                                                                                                                          >  1:397503/1‑251 (MQ=255)
                                                          gcgcgcTGCGTGAAGACTGGCGTGGCTGGTTGGCCGTTAGCCCGGAAAAGCCCGGAAGTACCGAAGCGCAAT‑GGAATGTCTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGATCCAGCCTAGCGCCTGTGACAAGTTATTTAGTGTCTGGCGGGCGGCAGGTAAACAAGATCCGCTGTAGTATTTAGAGCGTATCCGTCTGGTGATg                                                                                                                                                                                          <  2:318062/249‑1 (MQ=255)
                                                            gcgcGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGaa                                                                                                                                                                                        <  1:184589/251‑1 (MQ=255)
                                                                ggcgTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCg                                                                                                                                                                                    <  1:499833/251‑1 (MQ=255)
                                                                     gAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GGAATTACTACTATGCAAAATGGAACACCGGGCAGAGAGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCCACCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGGAAACGGGTa                                                                                                                                                                                >  1:74645/1‑250 (MQ=255)
                                                                       aGACTGGCGTGGCTTTTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGt                                                                                                                                                                                                                                       <  2:279982/193‑1 (MQ=255)
                                                                       aGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGt                                                                                                                                                                                                                                       >  1:279982/1‑193 (MQ=255)
                                                                        gACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAaca‑c                                                                                                                                                                           <  1:569744/251‑1 (MQ=255)
                                                                              ggttgCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCt                                                                                                                                                                     <  2:224920/247‑1 (MQ=255)
                                                                                  gCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCCAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTgga                                                                                                                                                                  >  2:528498/1‑249 (MQ=255)
                                                                                     tGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAAc                                                                                                                                                               <  1:533931/250‑1 (MQ=255)
                                                                                          gCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTATGGCGTGCGTCAGGTAAACAAGATCCGATGGCGTATTTAGAACGAATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCATGGTAACAGGGCt                                                                                                                                                         >  2:316262/1‑251 (MQ=255)
                                                                                              ttAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTgg                                                                                                                                                       <  1:627777/249‑1 (MQ=255)
                                                                                                  cccGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTagag                                                                                                                                                                                                            >  1:100392/1‑193 (MQ=255)
                                                                                                  cccGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAAAACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCGGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTagag                                                                                                                                                                                                            <  2:100392/193‑1 (MQ=255)
                                                                                                   ccGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATGACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTggcagggca                                                                                                                                                <  1:430180/251‑1 (MQ=255)
                                                                                                         acgCCCGGAACTACCGAAGCGCAAT‑GTATTTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGACGCCTGGCAAGGGGCGAAAGAGCTGTAGCTAACCGGC‑AAGAGCACGCCTAACGCCTGTGCCAAGTTATTTAGTGTCTGGCGTCCGTCAGGTAAACAAGATCCGCTGGCGTATTTTGAGCGTATCCGTCTGGCGATGAAAGCGGGAAACA‑CTGCCCTGGATACAGTGCTggcagggcagtcgcc                                                                                                                                          >  2:267762/3‑246 (MQ=255)
                                                                                                         aaGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGAtgcc                                                                                                                                          >  1:270930/1‑251 (MQ=255)
                                                                                                          aacccGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTCGTGTCAGGCATGCGTCAGGTAAACCAGATCCGCTGGCGTATTTCGAGCGTATCCGGCTGGCGATGAAAGCGGGTA‑CAGCTGGCCGTGTAACGTTGCTGGCAGGGCCGAtgcct                                                                                                                                         >  2:243136/3‑251 (MQ=255)
                                                                                                            cccGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGAtgcctgc                                                                                                                                       <  1:59492/251‑1 (MQ=255)
                                                                                                                  ccTACCGAACCGTAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCATAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCGAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTa                                                                                                                                 <  2:559521/250‑1 (MQ=255)
                                                                                                                      ccGAAGCGCAAT‑GTAATTACTACTAGGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAg                                                                                                                             <  1:651609/251‑1 (MQ=255)
                                                                                                                          aGCGCAAT‑GTAATTACTACTCTGCAAAATGGAACACCGGGCAGAGTGAAGACGCCTAGCAAGAGGCGAAAGAGCTGAGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGCATGGCGTGCGTCAGGAAAACAAGATACGCTGGCGTATTTAGAGCGTAACCGTCTGGCGATAAAATCGGGTAACA‑CTGGCCTGGCAACAGTGCTGGCAGGGGAAGTGACTGCCGATTACCacacg                                                                                                                          >  2:206521/1‑246 (MQ=255)
                                                                                                                          aGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACt                                                                                                                          >  1:512929/1‑250 (MQ=255)
                                                                                                                          aGCGCAAT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTATGGCTAACCGGC‑AAGAGCACGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGAGGAGGCCTGCCGATTACCAGACtt                                                                                                                         >  2:86515/1‑250 (MQ=255)
                                                                                                                                         actacgttgcaaaagggagtacCGGGCATAGTGAATAAGCCTGCCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AATATCCAGCCTAACTCCTGTGACAAGTTAGTTAGTTTCTGGCGTGCGTCAGGTAAACAAGTTCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGATAGCGGGTAACA‑CTGGTCTGGTAACAGTTCTGGCAGGGCAGATTCCTGCCGAGTTCCAGACTATCGCC‑TCGGCAAt                                                                                                           <  2:279981/230‑1 (MQ=255)
                                                                                                                                            actaTGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCa                                                                                                         <  1:287449/249‑1 (MQ=255)
                                                                                                                                                 gcaaaaTGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACt                                                                                                  <  1:366454/248‑1 (MQ=255)
                                                                                                                                                  caaaaTGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTg                                                                                                 >  1:537822/3‑251 (MQ=255)
                                                                                                                                                  caaaaTGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTg                                                                                                 <  1:447467/249‑1 (MQ=255)
                                                                                                                                                  caaaaTGCAACACCGGGCAGAGTGAAGCAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGACAGCCTAACGCCTGAGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCAGCTGGCGTATTTAGAGCGTATCCGATTGGCGGTGAAAGCGGGTAACA‑CTGGCCTGGTAACCGTGCAGGCAGGGCAGATGCCTGCAGCATACCAGAATATCGCC‑GCGTCAATCAATTCCCTg                                                                                                 >  2:294919/3‑251 (MQ=255)
                                                                                                                                                       tGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGTAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAACACGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGCCTATCTCC‑TCGGCAATCATTTTACTGGCgg                                                                                             >  2:97026/1‑249 (MQ=255)
                                                                                                                                                       tGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGa                                                                                             >  1:464228/1‑250 (MQ=255)
                                                                                                                                                            cgcctggCAGAGTTAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGTC‑ATGTGCCAGCCTATCTCCTGTGACAAGTTATTTATTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCGGGCTATGAAAGCGGGTACCT‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACAGGAGAACAAcc                                                                                       <  2:358856/246‑1 (MQ=255)
                                                                                                                                                              ccGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCt                                                                                     <  1:71199/251‑1 (MQ=255)
                                                                                                                                                               cGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCAAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGCGCGTCAGGTAAACACGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCGGCCGATTACCAGGATATTGCC‑TCGGCAATCagtgtaatggcgaatacacct                                                                                     >  2:223829/1‑230 (MQ=255)
                                                                                                                                                                gggCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACCACCCTa                                                                                    >  1:410396/1‑250 (MQ=255)
                                                                                                                                                                   cAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATAc                                                                                <  1:119223/251‑1 (MQ=255)
                                                                                                                                                                    agagTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGaacaa                                                                                         <  1:461247/241‑1 (MQ=255)
                                                                                                                                                                    agagTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACg                                                                               >  1:521152/1‑251 (MQ=255)
                                                                                                                                                                    agagTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGAGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGCTGCCTGCCGATTACCAGACTATCGCC‑CCGGCAATCATTTCCCTGGAGaaccaaactcgga                                                                                 >  2:461247/1‑239 (MQ=255)
                                                                                                                                                                     gagTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACgg                                                                              >  1:8359/1‑251 (MQ=255)
                                                                                                                                                                     gagTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTA‑CGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCGGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGt                                                                             <  1:407059/251‑1 (MQ=255)
                                                                                                                                                                      agTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGATGGCAGGGCAGATGCCGTCCGGTTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCCAACAACCCTAATACggg                                                                             >  2:477564/1‑250 (MQ=255)
                                                                                                                                                                         gaagaaGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACt                                                                          >  1:565571/1‑251 (MQ=255)
                                                                                                                                                                           agacgCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑CAAAGCCAGCCTAACGCCTGAGACTAGTTATTTAGTGTCTGGCGTGCGTCAGGAAAACACGATCAGCTGGCGTATTTCGAGCGTATCCGTCTGGCGCTGCCAGCGGGTAACA‑CTGGCCTGGTAACCGTGCTGGAAGGGCAGATTCCTGCCGACTACCAGCCCAGCGCC‑TCGGCACTCCTTTCACAGGCGAAAAACCCTaaaaaggtcc                                                                           >  2:672321/5‑240 (MQ=255)
                                                                                                                                                                            gaagACTGGCAAGGGGCGAAAGAGCGGTGGCTAACCGGC‑AAGCGCCAGCCTACCGCCTGTGACAAGTTAATTAGGGTCTGGCGTGCGTCAGGTAAACAAGATCCGCGGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CGGGCCTGGTAACGGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGAc                                                                       <  2:165437/251‑1 (MQ=255)
                                                                                                                                                                                   ggCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AGGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCTTGCGTCAGGTAAACAAGATCCGCTGGCGTATTGAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTcg                                                                  <  2:566074/249‑1 (MQ=255)
                                                                                                                                                                                            cGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACt                                                       <  1:335787/251‑1 (MQ=255)
                                                                                                                                                                                            cGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACt                                                       <  1:528498/251‑1 (MQ=255)
                                                                                                                                                                                                     gggtgCTAACCGGC‑AATAGCCAGCCTAACGACTGCGGCAAGTTATTTAGTGACTGGCGTGCGTAAGGGAATCAAGATCCGTGGGCGTATGTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACG‑CTGGCCTGGTAATAGTGCGGTCAGGGCAGATGCCTTCCGATTACCAGACTTTCGCC‑TCGGCACTCATTTCACTGGCGAACAACCCGAATACGGTACTTACCTTCGCGCGTACAACTGGCGTGAca                                              <  2:51851/247‑2 (MQ=255)
                                                                                                                                                                                                       tgGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACCAGTTATATAGCGTCTGGCGTGCGTCCGGTAAACAAGATCCGCAGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGAAACA‑CAGGCCTGGTAACAGTGCTGGCAGGGAAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTAACTGGCGAACAAACCTAAACCGGACCTGACCTTCGCGCGTACAAcagtcacgaccgc                                            >  2:257397/1‑239 (MQ=255)
                                                                                                                                                                                                           tAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGAttt                                         <  1:600263/250‑1 (MQ=255)
                                                                                                                                                                                                            aaCCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTGTTGTCTGGCTTGCGTCAGGTAAACAAGAGCCGCTGGCGTATTTAGAGCGTATCCTACTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTTCCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTTACCTTCGCGCGTACAACTGGCGCGACCGATTTTa                                       <  2:412412/251‑1 (MQ=255)
                                                                                                                                                                                                                  c‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTc                                 >  1:236470/1‑251 (MQ=255)
                                                                                                                                                                                                                    aaGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGCTTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGACCAACCCTAATACGGTACCTACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTc                                 >  2:509511/1‑250 (MQ=255)
                                                                                                                                                                                                                    aaGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGAGTCTGACGTGCGTCAGGTAAACCAGATCCGATGGCGTATTTAGCGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCGGGCAGGGCAGATGCCTGCCGATTAACAGACAATCGCA‑TCGGCAATCATTTCACTGGCGAACAACCCTCATACGGTACTGACCTTCGCGCGGACACATGGCGACACCGGTTTTAcccatc                                 >  2:441524/1‑247 (MQ=255)
                                                                                                                                                                                                                    aaGAGCCAGCCTAACGCCTGAGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAAAGCGTAACCGACTGACGATGAAAGCGGCAAACA‑ATGGCCTGGGACCAGACCTGGCAAGGCATAAGGCTGGCGATTACCAGACGATCGCCACCGGCAA‑CATTTCACAGGCGGACAAACCTAAAACGTCGCTGACCTTGGCGCGTATGACTGCAGCGACCGAatttcacagac                                 >  2:194549/1‑243 (MQ=255)
                                                                                                                                                                                                                       ggccagccTAACGCCTGTGACAATTTATTTAGTTTCTGGCTTGCGTAAGGTAAACAAGCTCCGCTGGCGTATTTAGATCGTATCCGTCTTGCGATGAACGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAAt                             <  2:618920/250‑1 (MQ=255)
                                                                                                                                                                                                                       agccatccTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAAt                             <  1:97887/251‑1 (MQ=255)
                                                                                                                                                                                                                       agccagccTACCGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGAATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAAt                             >  1:680231/1‑251 (MQ=255)
                                                                                                                                                                                                                       agccagccTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCaact                             >  1:183180/1‑249 (MQ=255)
                                                                                                                                                                                                                        gccagccTACCGCCTGTGACAAGTTATCTAGTGTCGGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGTCGCGACCGATTTTACCCGTCAAATg                            <  2:464228/251‑1 (MQ=255)
                                                                                                                                                                                                                         ccagctTATCGCCTGTGTCAAGTTATTTATTGTCTGGTGGGCGTCAGGTAAACAAGTTCCGCTAGCGTAGTTAGAGCGTTTCCGTCAGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGGCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACTGTACTGACCTTCGCGCGTTCAACTGGCGCGACCGATTTTAccctgtaaatgg                           <  2:488034/251‑10 (MQ=255)
                                                                                                                                                                                                                           agccTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGTCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAg                                                                                                                             <  1:584519/152‑1 (MQ=255)
                                                                                                                                                                                                                           agccCAACGCCTGGGACTAGTTTTTTAGTGGCTGTCTGGCGTCAGGTAAACAAGAACCGCAGGCGTATTTAGAGCGTATCCGTCTGGCGAGGAAAGCGGGTAACA‑CTGGCCTGGTAACATTGCTGGCAGGGCAGATGCCTTCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATggcg                         <  2:102375/251‑1 (MQ=255)
                                                                                                                                                                                                                            tccTAACGCCTGAGACAAGTTATTTGGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAACCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCCGCCGAGTACCAGACTAACGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACTGTACTAAACTTCGCGCGTCCCACTAGCGCGCCCTATTTAACCCGTCACATggcgg                        >  2:551289/2‑251 (MQ=255)
                                                                                                                                                                                                                              cTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGTGTTTTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CAGGACTGGTAACAGTGCGGGCAGGGCAGGTGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATAAGGGACTGACCTTCCCGCGAACCACTGGCGCGACCGAttttcaccgtccaaagtcggcg                      >  2:523354/1‑233 (MQ=255)
                                                                                                                                                                                                                                    ccTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCg                >  1:178693/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gtgACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGCTTACCCGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGAACTGACCTTCGAGCGTACAACTGGCGCGACCGATTTTACCCGTCCAATGGCGGGGGTgggggtt             >  2:296564/1‑246 (MQ=255)
                                                                                                                                                                                                                                              gTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAgtg      <  1:114705/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   ttAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGgc  <  1:477564/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
CGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAAT‑GTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCC‑TCGGCAATCATTTCACTGGCGAACAACCCTAATACGGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGC  >  NC_000913/4631000‑4631490

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑