Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,440,768 G→A L150L (CTG→TTG)  paeA ← putative cadaverine/putrescine exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,440,7680GA97.7% 259.8 / ‑2.9 87L150L (CTG→TTG) paeAputative cadaverine/putrescine exporter
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (38/45);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

CATCTTCGTGGAGATCAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCGCC  >  NC_000913/4440525‑4441012
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
cATCTTCGTGGAGATCAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTg                                                                                                                                                                                                                                                   >  1:465277/1‑251 (MQ=255)
     acgtgggaaccaacccAAATCACGTTTTCAGGCGGTGTCTTTGAAGACGGACCGGTACGTGATTCCAGCTCAAAGCCGTTTTCAATCAGCCCGTGTTCCTGTTTACGTA‑ACACCCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑c                                                                                                                                                                                                                                                 <  2:399961/233‑1 (MQ=255)
      cGTGGAGATCAAACCAAATCACGTTTTCCCGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCCCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTA‑ACCCGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAAc                                                                                                                                                                                                                                             <  1:569471/250‑1 (MQ=255)
                   ccAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGTAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAgc                                                                                                                                                                                                                               >  1:508985/1‑251 (MQ=255)
                                   cgcgGTGTAATTGAAGACGGAAAGGAACGGGATTCCAGCTAAACGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACAGGCTTCAACTACCGCGTAGCAGTCATCAGAAGTGACGTCATCTTT‑ACGTACCATTGGCAGTTTGATCATCCGGAAAATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGAG‑CACACGAACAGGAAGaggaccctcgggttgacgac                                                                                                                                                                                                               >  2:572311/1‑232 (MQ=255)
                                    gcgGTGACATTGAAGACGGAACGCTCCGGGATTCCAGCCCAAAGACGTTTTCAATCAGCGCGTTTGCCTGTTTACGGA‑ACACGACCGCCAGCGCACCGGCTTCCACTACCGCGGAGATGTCATCAGAAGTGATGGCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATa                                                                                                                                                                                                              <  1:266922/251‑1 (MQ=255)
                                      ggTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACg                                                                                                                                                                                                            >  1:28974/1‑251 (MQ=255)
                                                     ggAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTc                                                                                                                                                                                             <  1:218228/251‑1 (MQ=255)
                                                              taggctgccAGCACAAAGACGTATTAAATCAGCTCGTGTTCCTGTTTACGTA‑ACAAACCCGCAAGAGAACCGGCTTCCACTACCTCGTAGCTGTCATCAGGAGGGATGTCATCTTT‑ACGTACCATTGGCAGTTTGACAATACGGGAGTTTATGTTTGCCAGGCCGCTGGATAACCACCCCAAAGGAGTG‑CCAACGCACAGGCAGCAGCACATCGGGGTGATTATACTCCACGGCACAGCTTCGGGCGCact                                                                                                                                                                                    >  2:591106/8‑249 (MQ=255)
                                                              cGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTA‑CCACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGACAGGCCGTTGAAGAACCACACCAAAGGGGAG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGAAACGCAAAGCCACGGCTTTTGGAGCaac                                                                                                                                                                                    >  2:402906/1‑250 (MQ=255)
                                                                ggATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATc                                                                                                                                                                                   >  1:574498/1‑250 (MQ=255)
                                                                  aTTCCAGCTCAAAGCCGTTTCCAATCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATa                                                                                                                                                                                <  1:517949/251‑1 (MQ=255)
                                                                     ccAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCg                                                                                                                                                                             <  1:62817/251‑1 (MQ=255)
                                                                     ccAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTGACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATAcc                                                                                                                                                                              <  1:365932/250‑1 (MQ=255)
                                                                       aGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAACCATACCgct                                                                                                                                                                           >  2:268285/1‑249 (MQ=255)
                                                                        cccaaaaGGCGTTTTCAACCAGCCCGTTTCCCTGTTTCCGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATg                                                                                                                                                                          <  2:71022/247‑1 (MQ=255)
                                                                          acAAAGACGTGTTCAATCACCACGGTTTCTTGTTAACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGGGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCg                                                                                                                                                                        <  2:607746/250‑1 (MQ=255)
                                                                            aaaGACGATTTCAATCAACTCGTGTTCCTGTTTACGTA‑AAACGCACGCCAGCGCACCGGCTTCCACTACCGCGTTGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGATGAAGAACCACACCACCGGGGTG‑CAAACGTACAGGCAGAAGCGCATTGGGTTGGTGATACGCAAAGCCACAGCTTCTGGCGCATTCATACCGATGGGtt                                                                                                                                                                      >  2:81826/1‑251 (MQ=255)
                                                                                  gTTTTCAACCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCgggg                                                                                                                                                                <  1:505245/251‑1 (MQ=255)
                                                                                   ttttCAATCTGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCATCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCAGACCGATGCGTTTAgggg                                                                                                                                                                <  1:674394/250‑1 (MQ=255)
                                                                                       cAATCAGCTGCTTTTCCTGTTTACGTAGAC‑CGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGGTTGAAAATACGTAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTaaa                                                                                                                                                           <  2:465277/251‑1 (MQ=255)
                                                                                         aTCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCCCACCAACGGGGTG‑CAAACGTACAGGCCGAAGCGCCTCGGGTTGATGCTACGCAAAGCCACAGCTTCTGGCGCCATCATACCGATGCGTTTc                                                                                                                                                                    >  2:176978/1‑240 (MQ=255)
                                                                                         aTCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTc                                                                                                                                                                    <  1:176978/240‑1 (MQ=255)
                                                                                           catcacGTGTGCCTTTTTACGTA‑ACACGACCGCCAGCGCACCGGCTTCCACTACCTCGTAGATGTCAGCAGAAGTTATGTCATCGTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTT‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGGTTCGGGGTTAAATcc                                                                                                                                                        <  1:434169/245‑1 (MQ=255)
                                                                                           cAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCg                                                                                                                                                       <  1:277071/251‑1 (MQ=255)
                                                                                           agatgcgtgttcattgttACGTA‑ACCCCCACGCCAGCTAATCGGCTCCCCCTACCGAGTAGCTGTCATCTGCATTGATGTCATCTTT‑ATGTCCCATTGTCAGTCTGAAAATACGGAAGATTATGTTCGTCAAGCCGTTGAAGATCCACCTCATCGGGGTG‑CATACGTACAGGCAGAAGCGCATCGGGTTGATGAGACGCAAATCCACAGCTTCTGGTGCCATCATCCCGATGCGTTTCGGTGTTAAATcc                                                                                                                                                        <  2:371937/234‑1 (MQ=255)
                                                                                             gCTCGTGGTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑Ca                                                                                                                                                    <  1:423085/251‑1 (MQ=255)
                                                                                                    ttacTGTTTACGTC‑ACAAGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCAACAGAAGTGAGGTCATCTTT‑ACGTACCATTGACAGTTTGAAAATACGGAAGATTATGTTCGCCCGGACGGTGAAGAACCACACCAACGGGGAG‑CAAACGTACAGGCAGAAGCGCATCGGGATGATGAAACGAACAGCCACAGCATCTGGCGCAATCATAGCGAGGCGTATCAAGGTTAAATCAG‑CAAACAg                                                                                                                                               >  2:150090/4‑249 (MQ=255)
                                                                                                    ttCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAgg                                                                                                                                              >  1:638907/1‑250 (MQ=255)
                                                                                                     tCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGAt                                                                                                                                            >  1:522546/1‑251 (MQ=255)
                                                                                                       cTGTTTACGTA‑ACACGCACGCCAGCGCACCGGCTTCCACTACCGCGTAGATTTCAACAGAAGTGAAGTCATCTTT‑ACGTACCATTGGCAGTTTGCAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGAG‑CAAACGAACAGGCAGAAGCTCATCGGGTTGATGATAAGCAAAGACACAGCTTCTGGCGAAATCAAACCGATGCGTTTCGGGGTTAAAGCCGACAAA‑AGGATa                                                                                                                                           >  2:448671/1‑250 (MQ=255)
                                                                                                        tGTTTACGTA‑CCACGACCGCCAGAGCACAGGATTCCAATACCTCGTAGATGTCATCAGAAGTGATGTCCTCTTT‑ACGTACCCTTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGACGTTGGAGAACCACACCAACGGGGTGACAAA‑GTACAGGCAGAAGCGCATGGGGTTGCTGATAAGCAAGGCCAACGCTTCAGGAGAAATCATACCGAGACGGTTCGGGGTTAAATACG‑AAAACAGGATaa                                                                                                                                          >  2:679513/1‑250 (MQ=255)
                                                                                                        tGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGAACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAACCAGGATaac                                                                                                                                         >  1:495841/1‑250 (MQ=255)
                                                                                                         tttaccGTA‑ACCCCCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCAGTGGCAGTTTGAAAATACGGAAGATTATGTGCGCCAGGCCTTTGAAGAACCACACCAACGGTGTG‑CAAAc                                                                                                                                                                                                                                             <  2:663606/146‑1 (MQ=255)
                                                                                                         gTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAAc                                                                                                                                                                                                                                             >  1:663606/1‑151 (MQ=255)
                                                                                                                    cccgcccgccAGCGCACCGGCTTCCACTTCCGCGGAGATGGCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTAGGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAACCGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCa                                                                                                                              <  2:522546/249‑1 (MQ=255)
                                                                                                                    cacGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCa                                                                                                                              <  1:262383/251‑1 (MQ=255)
                                                                                                                       ccccgccATCGCACCGGCTCCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGGACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCTCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCAGCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAGCATGCCAGCCACTa                                                                                                                           <  2:110571/250‑1 (MQ=255)
                                                                                                                          cgccAGCGCACCGGCTTCCACCACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACg                                                                                                                        <  1:679513/251‑1 (MQ=255)
                                                                                                                           gccAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACg                                                                                                                        >  1:661569/1‑250 (MQ=255)
                                                                                                                            ccAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACgag                                                                                                                      <  1:2194/251‑1 (MQ=255)
                                                                                                                              agcgacccGGTTTCCTCTACCGAGGAGATGGCAACAGAAAGGATGTCATCTTTTAGGT‑CCATTGGCAGTTTGAAAATACGGAATATTATTTTCGCCAGGCCGGTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑aa                                                                                                                   <  2:495841/245‑1 (MQ=255)
                                                                                                                              aGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑aa                                                                                                                   <  1:85118/251‑1 (MQ=255)
                                                                                                                                  cACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAACGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGGTAAACATGCCAGTCACTatagggaaaggg                                                                                                                >  2:156362/1‑240 (MQ=255)
                                                                                                                                      ggCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGaa                                                                                                           <  2:638435/251‑1 (MQ=255)
                                                                                                                                           ccACTACCGCTTAGATGTCATCAGAAGTGATGTCATTTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGAGTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATaaa                                                                                                       <  2:213749/250‑1 (MQ=255)
                                                                                                                                                              atcaGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGct                                                                                   <  1:545540/251‑1 (MQ=255)
                                                                                                                                                                      ggATGTCATCTGT‑ACGTACCAGTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCCGGCCGTTGAAGATCCACACCAACGGGGTG‑CAAACGGACAGGCCGAAGCGCATCGGGTGGATGATACGCAAAGCCACAGCTTCTGGCGCAAACATACCGATGCGTTTCGGGGTCAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGATCAGTTGCTCAGAGAGTTCTGCCGAc                                                                           <  2:98700/250‑1 (MQ=255)
                                                                                                                                                                        atttCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACat                                                                         <  2:635359/248‑1 (MQ=255)
                                                                                                                                                                        aTGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCCGGCCGTTGAAGAACCACACCAACCGGGCG‑CAAACGTACAGGCAGAAGCGCATCGGGATGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGATTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTGACGAGCAAAG‑GAGAATAAAGCTCAGTTGCTCAGAGAGCGCTGCCGACat                                                                         >  2:153154/1‑251 (MQ=255)
                                                                                                                                                                            cATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGATCCCCACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGATAGCTCTGCCGAAAtata                                                                      <  1:150090/250‑1 (MQ=255)
                                                                                                                                                                                  t‑ACGTACCATTGGCAGTTTGAAAATCCGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGga                                                                <  2:574498/250‑1 (MQ=255)
                                                                                                                                                                                    aCGTACCAATGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGgag                                                               <  1:3224/250‑1 (MQ=255)
                                                                                                                                                                                     cgtcccATTGGCAGTTTGAAAATACGGAAGCTTATGTTCGCCAGGCCGTTGAAGAAACACACCAACGGGGTG‑CAACCGTACAGGCAGAAGCACATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCACTCATACCGATGCGTTTCGGGGTTAAATCCG‑AAAACAGGATCAACATGCCAGACACTAACGAG‑CAAGAGAGAATAAAGCTCAGTTGCtagg                                                                                         >  1:661713/5‑220 (MQ=255)
                                                                                                                                                                                        aCCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACa                                                           <  1:238918/250‑1 (MQ=255)
                                                                                                                                                                                              ggCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCtgtg                                                    <  1:344982/251‑1 (MQ=255)
                                                                                                                                                                                                  gTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGaaa                                                 >  1:236074/1‑250 (MQ=255)
                                                                                                                                                                                                  gTTTGAAAAGACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGAGTG‑CAAACGTACAGGCAGAAGCGCATCGGGATGATGATCCGCAAAGCCCCAGCTTCTGGCGCAGTCATACCGATGAGTTTCGAGGTTAAACCCG‑CAAACAGGAGAACCATGCCAGCCACTAACGAA‑AAAGAGAGAATAAAGCACAGTTGCTCAGCGAGCTCTGCAGCCATATACAGGGAACACAGGCGGTGaaa                                                 >  2:502794/1‑250 (MQ=255)
                                                                                                                                                                                                      taaaaTACGGACGCTTATGTTCGCCAGGCCGTTGCAGCACCCCACCAACGGGGGC‑CAAACGTACAGGCAGAAGCGCATAGGGTTGATGATACACAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTACATCCG‑CAACCAGGATAAACATGCACGTCACTACCGAG‑GAAGAGAGAATAAAGCTCATTTGCTTCGAGCGCTCTGCCGACATATAGCGGGAGAACATGCGGTGTAAAgatt                                            >  2:164726/2‑248 (MQ=255)
                                                                                                                                                                                                      gAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAAt                                                                                                                                                          <  2:567838/143‑1 (MQ=255)
                                                                                                                                                                                                      gAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAAt                                                                                                                                                          >  1:567838/1‑143 (MQ=255)
                                                                                                                                                                                                       aaaaTACGGAAGATTATGTTCGCCACGCCGTTGAAGAACCACACCAACGCGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGCTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTAGCTCAGCGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTgg                                          >  2:317503/1‑251 (MQ=255)
                                                                                                                                                                                                          aTACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTgg                                          >  1:655732/1‑248 (MQ=255)
                                                                                                                                                                                                            aCGGAAGATTATTGTCGCCATGTTGTTGAAGATCCACACCAACGGGGTT‑CAAACGTTCAGGCATAAGCGCATCGGTTTGATGATACGCAAAGCCACAGCGTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CATACAGGATAAACATGCCAGTCACTAACTAG‑ATAGAGAGAATAAAGCTGAGTTGCTCCGAGTGCTCTGCCGACACATAGCGGGAGAACAGGCTTGGAAAAT‑CTGGAGaaa                                     <  2:533954/251‑1 (MQ=255)
                                                                                                                                                                                                            aCGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGaa                                      >  1:506353/1‑250 (MQ=255)
                                                                                                                                                                                                            aCGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTAACTAACGAG‑AAAGAGAGAATAAAGCTCAg                                                                                                 >  1:584235/1‑192 (MQ=255)
                                                                                                                                                                                                              ggcaGATTAAGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CGAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGATCAGCGAGCTCTGCCGACATATAGCGGGAGAACAGGCGGTGAAAAG‑CTGG‑Gaaaaat                                  >  2:295819/4‑249 (MQ=255)
                                                                                                                                                                                                               gAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCAGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATg                                  <  1:268285/251‑1 (MQ=255)
                                                                                                                                                                                                               gAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATg                                  >  1:491444/1‑251 (MQ=255)
                                                                                                                                                                                                                aaGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGc                                 >  1:623118/1‑251 (MQ=255)
                                                                                                                                                                                                                  gaatATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAATCAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCg                               <  2:669678/248‑1 (MQ=255)
                                                                                                                                                                                                                   aTTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGc                              <  1:580580/251‑1 (MQ=255)
                                                                                                                                                                                                                   aTTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGAAACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGCG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACCTATAGCGGGAGAACAGGCTGTGAAAGG‑CTGGAGAAAAAGCCGc                              >  2:7644/1‑251 (MQ=255)
                                                                                                                                                                                                                   aTTATGTGCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGAGCATCGGGTGGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGGTGCTCAGAGAGCGCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGc                              <  2:661569/251‑1 (MQ=255)
                                                                                                                                                                                                                    ttATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGc                              >  1:467258/1‑250 (MQ=255)
                                                                                                                                                                                                                     tATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGAGACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGCa                             <  2:3237/250‑1 (MQ=255)
                                                                                                                                                                                                                      aTGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGCAt                            <  1:317503/250‑1 (MQ=255)
                                                                                                                                                                                                                           cGCCAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGCATCACcg                       >  1:557110/1‑250 (MQ=255)
                                                                                                                                                                                                                            gcccggCCGTTGAAGAACCACACCAACGGGGTG‑CAAAAGTTCAGGAAGAAGCGCATCGGGTTGATGATACGCAAAGCCTCAGCTTCTGGCGCAATCATACCGATGCGCTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCGCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGCATCTCcgac                     <  2:377633/247‑1 (MQ=255)
                                                                                                                                                                                                                              cAGGCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGCATCACcgacg                    >  1:552122/1‑250 (MQ=255)
                                                                                                                                                                                                                                ggCCGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CACACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGCATCACCGACGAt                  >  1:326276/1‑250 (MQ=255)
                                                                                                                                                                                                                                  ccGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGAAGAAAATGCCGCATCACCGACGATAcc               <  1:631393/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ccGTTGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGAAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCGTCTGGGGCAATCATACCGATCGGTTGCGGGGTTAAAGCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTGGCGCAGAGAGCTCTGCCGACATATAGCGGTAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGCATCACCGACGATAcc               <  1:397284/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tttGAAAAACCACACCAACGGGGTG‑CAAACGTCCAGGAAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGATTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑CAAGAGAGCATAAAGCTCAGGTGATAAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGGCGGGA‑AACATACCGCATCACcgaccatcccg              >  2:461519/2‑243 (MQ=255)
                                                                                                                                                                                                                                      tGAAGAACCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGCATCACCGACGATAccgccg           >  1:131792/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gaaCCACACCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTACCGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGCATCACCGACGATACCGCCgaga        >  1:273395/1‑250 (MQ=255)
                                                                                                                                                                                                                                               ccacCAACGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTGCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCGcc  <  1:297198/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
CATCTTCGTGGAGATCAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTA‑ACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTT‑ACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTG‑CAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCG‑CAAACAGGATAAACATGCCAGTCACTAACGAG‑AAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAG‑CTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCGCC  >  NC_000913/4440525‑4441012

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