Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,160,733 C→T G20S (GGC→AGC)  sthA ← pyridine nucleotide transhydrogenase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,160,7330CT99.4% 550.9 / NA 156G20S (GGC→AGC) sthApyridine nucleotide transhydrogenase, soluble
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (68/87);  total (69/87)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 156 total aligned reads displayed.

GAGCGTTTCCCTGCAATATTTCACAGTGATTACGTTCGTAAAATCCCTGACGCATGCGCGTTTGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTG  >  NC_000913/4160439‑4161030
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                         
          cTGCAATATTTCACAGTGATTACGTTCGTAAAATCCCTGACGCATGCGCGTTTGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCAAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCCCCTTGCTTAACCAAGCCCATTGCAGCGCTTTCGCCGCCGGGGCCg                                                                                                                                                                                                                                                                                           >  2:1669/1‑301 (MQ=255)
               atatTTCACAGTGATTACGTTCGTAAAATCCCTGACGCATGCGCGTTTGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATACCGCTCGATAAATCCGACGCGCGCACCTTGATTAACCAGGCCCATTGCAGCGCTTTCGCCTCCGGGGCCgtacc                                                                                                                                                                                                                                                                                       >  2:23152/1‑296 (MQ=255)
               atatTTCACAGTGATTACGTTCGTAAAATCCCTGACGCATGCGCGTTTGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAAc                                                                                                                                                                                                                                                                                       >  2:779300/1‑300 (MQ=255)
                 atTTCACAGTGATTACGTTCGTAAAATCCCTGACGCATGCGCGTTTGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTa                                                                                                                                                                                                                                                                                    >  2:721185/1‑301 (MQ=255)
                     cacaGTGATTACTTTCGTAAAATCCCTGACGCATGCGCGTTTGTTGAGTAATCACGTTACCGGCATGGTTAAGGATATCGGCCAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAAGTCTATAAGTCGGCTGCCGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGAGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTAc                                                                                                                                                                                                                                                                                <  2:1283523/301‑1 (MQ=255)
                        aaaggtgtataagagacagtCCCTGACGCATGCGCGTTTGTTTATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTa                                                                                                                                                                                                                                                                              <  2:1235442/281‑1 (MQ=255)
                          tGATTACGTTCGTAAAATCCCTGACGCATGCGCGTTTGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTa                                                                                                                                                                                                                                                                              <  2:436151/298‑1 (MQ=255)
                          tGATTACGTTCGTAAAATCCCTGACGCATGCGCGTTTGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTa                                                                                                                                                                                                                                                                              <  2:511378/298‑1 (MQ=255)
                          tGATTACGTTCGTAAAATCCCTGACGCATGCGCGTTTGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTa                                                                                                                                                                                                                                                                              <  2:1196825/298‑1 (MQ=255)
                             ttACGTTCGTAAAATCCCTGACGCATGCGCGTTTGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCat                                                                                                                                                                                                                                                                        <  2:639350/301‑1 (MQ=255)
                                         aaTCCCTGACGCAGGCGCGTTTGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAAc                                                                                                                                                                                                                                                                                       <  2:874450/274‑1 (MQ=255)
                                           ccccTTACGCATGCGCGTTGGTTGATTAATCACTTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTTACGGCGTGACGGGGAGCTTTCGACGGTATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGaa                                                                                                                                                                                                                                                          <  2:1094821/300‑1 (MQ=255)
                                                 aCGCGTGCGCGTTTGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGc                                                                                                                                                                                                                                                    <  2:471617/301‑1 (MQ=255)
                                                          cgTTTGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAggg                                                                                                                                                                                                                                           >  2:1323905/1‑301 (MQ=255)
                                                            tttGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGc                                                                                                                                                                                                                                          <  2:603985/300‑1 (MQ=255)
                                                            tttGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGc                                                                                                                                                                                                                                          <  2:518663/300‑1 (MQ=255)
                                                            tttGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGc                                                                                                                                                                                                                                          <  2:1048849/300‑1 (MQ=255)
                                                            tttGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGc                                                                                                                                                                                                                                          <  2:1434372/300‑1 (MQ=255)
                                                            tttGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGCAATCGTAGGAATGTGGCATGGTAggg                                                                                                                                                                                                                                           >  2:67978/1‑299 (MQ=255)
                                                                ttgATTAATAACGTTATCGGCATTGTTAAGGATATCGGCAAAATAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGAGTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTAcc                                                                                                                                                                                                                                     <  2:440592/301‑1 (MQ=255)
                                                                 tgATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCt                                                                                                                                                                                                                                    <  2:1024216/301‑1 (MQ=255)
                                                                   aTTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTg                                                                                                                                                                                                                                   <  2:1261070/300‑1 (MQ=255)
                                                                     aaagcACGTGATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACGGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTc                                                                                                                                                                                                                                <  2:48698/297‑1 (MQ=255)
                                                                           cGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCTGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATAc                                                                                                                                                                                                                          <  2:572308/301‑1 (MQ=255)
                                                                            gataggggaaggggtAAGGGTATCGGCAAGAGGAGGGGGGAGCAGGCGGGAATGGTCGCTGTAAAGTGGGTGTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGGAGGGATGGTGCACCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGAGCGCACCTTGCGTAACCAGGCGCAGTGAAGCGCTTTCGCCGCCGGGGCCGGAACCAATTACTATGGCATCGTAATAGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATAc                                                                                                                                                                                                                          <  2:62236/286‑1 (MQ=255)
                                                                                      tGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGaa                                                                                                                                                                                                                                                          >  2:508501/1‑258 (MQ=255)
                                                                                      tGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTAACCATGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGTCTTACCTGTTCTTATACATAAAAGCAACa                                                                                                                                                                                                              >  2:428030/1‑301 (MQ=255)
                                                                                       ggTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTcgccgc                                                                                                                                                                                                                                                                                                   >  1:966812/1‑216 (MQ=255)
                                                                                       ggTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGGTCTTATACATAAAAGCAACa                                                                                                                                                                                                              >  2:542383/1‑301 (MQ=255)
                                                                                       ggTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGCCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGAGCCGGAACCTTTTACTATCGCATCGTAATTGTAGGCATGTGGCATGGTAGGGCTTACCTTTTCTTATACATACAAGCAACa                                                                                                                                                                                                              >  2:1434099/1‑301 (MQ=255)
                                                                                        gTTAGGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTc                                                                                                                                                                                                                                >  2:316115/1‑282 (MQ=255)
                                                                                        gTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGGGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACCTTTTGATAACGCTCGCTAACTGCGCCGCGCGCGCCTTCCTTAACCAGTCTCCTTTCAGCTATGTCGccgccgn                                                                                                                                                                                                                                                                                                >  1:1325154/1‑217 (MQ=255)
                                                                                        gTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCgg                                                                                                                                                                                                                                                                                                >  1:482227/1‑218 (MQ=255)
                                                                                        gTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCAt                                                                                                                                                                                                                                                  <  2:482236/264‑1 (MQ=255)
                                                                                        gTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTCTACATAAAAGCaaaaa                                                                                                                                                                                                             >  2:26294/1‑298 (MQ=255)
                                                                                             ggATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATgg                                                                                                                                                                                                        >  2:64791/1‑301 (MQ=255)
                                                                                                 agcGGCAAAAGAAGAGCGTAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATGCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACAGCGACGCGCGCACCTTGCTTAACCAGGACCATTGCAGCGCTTTCGCCGCCGGAGCCGGaa                                                                                                                                                                                                                                                                                        <  1:462669/215‑1 (MQ=255)
                                                                                                        aaaGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTAc                                                                                                                                                                                                                                                                                <  1:779286/218‑1 (MQ=255)
                                                                                                        aaaGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTAc                                                                                                                                                                                                                                                                                <  1:721172/218‑1 (MQ=255)
                                                                                                        aaaGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTAc                                                                                                                                                                                                                                                                                <  1:424183/218‑1 (MQ=255)
                                                                                                        aaaGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTAc                                                                                                                                                                                                                                                                                <  1:1669/218‑1 (MQ=255)
                                                                                                         aagaagAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATc                                                                                                                                                                                            <  2:1319999/301‑1 (MQ=255)
                                                                                                                   gAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGTCGGGGCCGGAACCTATTACTATGGCatcgn                                                                                                                                                                                                                                                                     >  1:1246320/1‑217 (MQ=255)
                                                                                                                        gTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAAtt                                                                                                                                                                             >  2:51606/1‑301 (MQ=255)
                                                                                                                          cGGGAATGGTCTCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGa                                                                                                                                                                           <  2:1081457/301‑1 (MQ=255)
                                                                                                                          cGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACTGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTg                                                                                                                                                                            <  2:1246371/300‑1 (MQ=255)
                                                                                                                          cGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTg                                                                                                                                                                            <  2:997259/300‑1 (MQ=255)
                                                                                                                            ggAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGAGACGGAGAGCTTTCGACGGGATGGTGCACCAGTGGGTGCAACAGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGGACATTGCTTAACCAGGCACATTGCAGCGCTTTCGCCGACGGGGCCGGAACCTATTACTATGGCATCGACATCGTAGTAATGGGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTCACATTTTATCGCGGGTAAGTCAATTGAta                                                                                                                                                                         >  2:1350742/1‑300 (MQ=255)
                                                                                                                              aaTGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGTAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGaa                                                                                                                                                                                                                                                          <  1:307847/218‑1 (MQ=255)
                                                                                                                                      cTGTAAAGTGGGTTTTTATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTGATCGTAGGAATGTGGCATGGTAGGGGTTACCTTGTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCAtt                                                                                                                                                                <  2:1424346/300‑1 (MQ=255)
                                                                                                                                      cTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTAACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATCACTATGGCATCGTAATCGTAGGAATGTGGCAt                                                                                                                                                                                                                                                  >  1:182098/1‑218 (MQ=255)
                                                                                                                                        gTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTa                                                                                                                                                                                                                                                                              >  2:1001152/1‑188 (MQ=255)
                                                                                                                                        gTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATgg                                                                                                                                                                                                                                                >  1:102110/1‑218 (MQ=255)
                                                                                                                                        gTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTAcc                                                                                                                                                                                                                                     >  2:1178252/1‑229 (MQ=255)
                                                                                                                                         tAAAGTGGGTTTTGATTGAATTCTCTAATGCGGCTGACGGCGGGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATc                                                                                                                                                            <  2:1187194/301‑1 (MQ=255)
                                                                                                                                                   nttgattgaATTCTATAATGCGGCGGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTAcc                                                                                                                                                                                                                                     <  1:1178210/217‑1 (MQ=255)
                                                                                                                                                    ttgattgaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATcc                                                                                                                                                 <  2:217362/301‑1 (MQ=255)
                                                                                                                                                       attgaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTg                                                                                                                                                                                                                                   <  1:1095014/216‑1 (MQ=255)
                                                                                                                                                          gaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCCACCGCCGCCAACATTTTGCTAACTCTCGATAACTGCGACGCGCGCACCTTGCTTAACCATTCCCATTTCAGCGCTTTCGCCGCCTGGGCCTGACCCTATTACTATTGCATCGTAATCTTCGGAATGTGGCATGGTAGGGCTTACCTGCTCtt                                                                                                                                                                                                                              >  1:649723/1‑218 (MQ=255)
                                                                                                                                                          gaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCAt                                                                                                                                                                                                                                                  >  2:458822/1‑198 (MQ=255)
                                                                                                                                                          gaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCAt                                                                                                                                                                                                                                                  <  1:458813/198‑1 (MQ=255)
                                                                                                                                                            aTTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGaa                                                                                                                                                                                                                                                          >  1:281040/1‑188 (MQ=255)
                                                                                                                                                            aTTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGaa                                                                                                                                                                                                                                                          <  2:281043/188‑1 (MQ=255)
                                                                                                                                                             ttCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGGAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATGTATCTGGCTATCTCAGGAGCGGc                                                                                                                                         >  2:282110/1‑300 (MQ=255)
                                                                                                                                                               cTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTTCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCACCGTACTCGTAGGAATGTGGCATGGTAGGGCTCACCTGTTCTTATACa                                                                                                                                                                                                                         >  1:889605/1‑218 (MQ=255)
                                                                                                                                                               cTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctt                                                                                                                                      >  2:855840/1‑301 (MQ=255)
                                                                                                                                                               cTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAACCGTAGGAATGTGGCATGGTAGGGCTTACCTGGTCTtata                                                                                                                                                                                                                           >  1:1069985/1‑216 (MQ=255)
                                                                                                                                                                                   tGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGTACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTAc                                                                                                                   <  2:1025257/300‑1 (MQ=255)
                                                                                                                                                                                     aCGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATcc                                                                                                                                                 <  2:74887/268‑1 (MQ=255)
                                                                                                                                                                                               ttCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACAt                                                                                                                                                                                                                        >  2:351236/1‑187 (MQ=255)
                                                                                                                                                                                               ttCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCAAACATGTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACg                                                                                                      <  2:85530/301‑1 (MQ=255)
                                                                                                                                                                                               ttCGACGGGATGGTGCCCCAGTGGGGGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGGACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACAt                                                                                                                                                                                                                        <  1:351230/187‑1 (MQ=255)
                                                                                                                                                                                                 cGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTTCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGACTGTGGCATGGTACGGCGTACCTGTTc                                                                                                                                                                                                                                >  1:1235179/1‑177 (MQ=255)
                                                                                                                                                                                                 cGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGTGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGCCGtt                                                                                                    >  2:382833/1‑301 (MQ=255)
                                                                                                                                                                                                 cGACGGGATGGTGCCCCAGTGGGTGCAACCGCAGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTc                                                                                                                                                                                                                                <  2:1235227/177‑1 (MQ=255)
                                                                                                                                                                                                      ggATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTAt                                                                                                <  2:1070015/300‑1 (MQ=255)
                                                                                                                                                                                                       gATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGANTGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTg                                                                                              >  2:1403752/1‑301 (MQ=255)
                                                                                                                                                                                                        aTGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGaa                                                                                                                                                                                                                                                          >  2:1328809/1‑144 (MQ=255)
                                                                                                                                                                                                        aTGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGaa                                                                                                                                                                                                                                                          <  1:1328756/144‑1 (MQ=255)
                                                                                                                                                                                                               cccAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATc                                                                                                                                                                         <  1:245758/218‑1 (MQ=255)
                                                                                                                                                                                                                       gTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCAATTCCAGCGCCTTCGCCGCCGGGGCCGGCACCTATTACTATAGCCTCGTCATCGTAGGAATGTGGTCGGGTAGGGCTTACCTGTTCTTATACATACAAGCCACACACTGGTAACATTTTATCGCGGTTAATTCAATTGATCCCCGTCa                                                                                                                                                                  >  1:841432/1‑217 (MQ=255)
                                                                                                                                                                                                                          cAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATAc                                                                            <  2:889624/300‑1 (MQ=255)
                                                                                                                                                                                                                             ccgccgCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACtctc                                                                        <  2:1148576/301‑1 (MQ=255)
                                                                                                                                                                                                                                ccgccAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATAGACGATACTCTCCGt                                                                     >  2:407373/1‑301 (MQ=255)
                                                                                                                                                                                                                                   ccAACAGTTTGATAACGCTCAATAACTGCGACGCTCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGTCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCta                                                                                                                                                      <  1:610043/217‑1 (MQ=255)
                                                                                                                                                                                                                                       naTTTTGGTAACGCGCGATAACGGCGACGCGCGCACCTTGCTTAACCAGGCGCATTGGAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTGTGGCATCGTAATCGTAGGAATGGGGCATTGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTATCATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATcc                                                                                                                                                 <  1:1263662/217‑1 (MQ=255)
                                                                                                                                                                                                                                                    cTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCg                                                   >  2:617733/1‑299 (MQ=255)
                                                                                                                                                                                                                                                    cTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGc                                                  <  2:348165/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cgcACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTAc                                                                                                                   <  1:51605/217‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cgcACCTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCAGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTTTTATACATAAAAGCAACAGAATGGTAACATTTTAGCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCGTTGTCTGTGTCTTAc                                                                                                                   <  1:449145/217‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cgaaCCTTGCATAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTAc                                                                                                                   <  1:23152/214‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ccTTGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTACAAGAACCGTCCaaa                           >  2:1096837/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                             tGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTTcc                                                                                                            >  1:482201/1‑217 (MQ=255)
                                                                                                                                                                                                                                                                             tGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGc                                                                                                             >  1:1029940/1‑216 (MQ=255)
                                                                                                                                                                                                                                                                             tGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTTTCCGTTTAAGCGGAATGTTTCCGCTGTACGTAAAAGAAACGGCCAAAGaa                        >  2:788875/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                             tGCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGaa                        >  2:192711/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                              gCTTAACCAGGCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCATAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGcc                                                                                                            >  1:73246/1‑216 (MQ=255)
                                                                                                                                                                                                                                                                                        gCCCATTGCAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGCCGTTTTan                                                                                                >  1:1314888/1‑217 (MQ=255)
                                                                                                                                                                                                                                                                                           cATTGCAGCGCTTTCGCCGCCGGGGCAGGAACCTATCACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATg          <  2:182098/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           aaTTGAAGCGCTTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATg          <  2:257640/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                         
GAGCGTTTCCCTGCAATATTTCACAGTGATTACGTTCGTAAAATCCCTGACGCATGCGCGTTTGTTGATTAATCACGTTATCGGCATGGTTAAGGATATCGGCAAAAGAAGAGCGGAGCAGTCGGGAATGGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTG  >  NC_000913/4160439‑4161030

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