Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,495,036 C→A 9.7% A343E (GCA→GAA)  nirB → nitrite reductase, large subunit, NAD(P)H‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,495,0360CA 9.7% 275.2 / 9.7 96A343E (GCA→GAA) nirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base C (50/35);  new base A (2/8);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 3.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCA  >  NC_000913/3494749‑3495323
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                               
tAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCa                                                                                                                                                                                                                                                                                                <  2:232273/289‑1 (MQ=255)
              aGTGTGGTCTGGACGTTGCTCCGCGTGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGCAGACGTAGGCGGTATTGGTAATGCGCCTGGTCGCCGCACTGGCGCACGTACTTTCGTTTa                                                                                                                                                                                                                                                                                  >  1:76787/1‑289 (MQ=255)
                     tCTGGACGTTGCTCCGCGTGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTATCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCCCGGTCGCACGCCTGGCGTACGTAGCTACGTTTACCTcgac                                                                                                                                                                                                                                                                           >  2:256035/1‑289 (MQ=255)
                     tCTGGACGTTGCTCCGCGTGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTcga                                                                                                                                                                                                                                                                            <  2:41171/288‑1 (MQ=255)
                      cTGGACGTTGCTCCGCGTGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTcgacg                                                                                                                                                                                                                                                                          >  1:454641/1‑289 (MQ=255)
                        ggACCTTGCTCCGCTGGGGGGTAGTGTCATTAATGAGTCCTGCCAGACTTCCGATACGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATGCTCGGTAGCGAAAACGCCTTGGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGAGGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGAACGTAGCTACGTTTACCTcgacga                                                                                                                                                                                                                                                                         <  1:256035/288‑1 (MQ=255)
                          aCGTTGCTCCGCGTGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAg                                                                                                                                                                                                                                                                      >  2:679333/1‑289 (MQ=255)
                          aCGTTGCTCCGCGTGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAg                                                                                                                                                                                                                                                                      >  2:630373/1‑289 (MQ=255)
                          aCGTTGCTCCGCGTGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAg                                                                                                                                                                                                                                                                      >  1:378465/1‑289 (MQ=255)
                              tGCTCCGCGTGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTaaa                                                                                                                                                                                                                                                                  >  2:230530/1‑289 (MQ=255)
                              tGCTCCGCGTGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTaaa                                                                                                                                                                                                                                                                  >  2:301904/1‑289 (MQ=255)
                                     cgTGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCt                                                                                                                                                                                                                                                           >  2:69651/1‑289 (MQ=255)
                                               ttGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACg                                                                                                                                                                                                                                                    >  2:41857/1‑286 (MQ=255)
                                                tGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTCGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGAACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCt                                                                                                                                                                                                                                                 >  1:346924/1‑288 (MQ=255)
                                                    aTTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCGACGTTTACCTCGACGACAGTAAAGAGATCTACAAACGCCTGATTg                                                                                                                                                                                                                                            >  2:588109/1‑289 (MQ=255)
                                                       aaTGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCa                                                                                                                                                                                                                                         <  2:699803/289‑1 (MQ=255)
                                                            tttcTGCCAGACTTCCGATACGGATATCGACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCACGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGaa                                                                                                                                                                                                                                    <  2:98312/286‑1 (MQ=255)
                                                                              tCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCt                                                                                                                                                                                                                   <  1:630373/288‑1 (MQ=255)
                                                                              tCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCt                                                                                                                                                                                                                   <  1:668957/288‑1 (MQ=255)
                                                                              tCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCt                                                                                                                                                                                                                   <  1:674631/288‑1 (MQ=255)
                                                                              tCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTc                                                                                                                                                                                                                  <  1:69651/289‑1 (MQ=255)
                                                                              tCCGGATATCTACGACATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCt                                                                                                                                                                                                                   <  2:651734/288‑1 (MQ=255)
                                                                              cacgataaccTACGCCATCGTTGAATGCCAAAGCGTGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCGACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCt                                                                                                                                                                                                                   <  1:588109/279‑1 (MQ=255)
                                                                                  gATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGTCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGACGGTGCTGACCTTACCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTTGCACGCGTGTCGCACGTAGCTACGTTTACCTCGACGAACGTAACGAGACCTACAAACGCCTGATTGTCAGCGAAGACAACAAACTCCTCATCGGt                                                                                                                                                                                                               >  1:487602/1‑288 (MQ=255)
                                                                                        tACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGAt                                                                                                                                                                                                                                                             >  2:541979/1‑236 (MQ=255)
                                                                                        tACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGAt                                                                                                                                                                                                                                                             <  1:541979/236‑1 (MQ=255)
                                                                                          cGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCAGCTCGGTGCGGTACtg                                                                                                                                                                                                      >  2:623786/1‑289 (MQ=255)
                                                                                           gCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACtgg                                                                                                                                                                                                     >  2:26271/1‑289 (MQ=255)
                                                                                             cATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTcga                                                                                                                                                                                                                                                                            >  1:514135/1‑216 (MQ=255)
                                                                                             cATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTcga                                                                                                                                                                                                                                                                            <  2:514135/216‑1 (MQ=255)
                                                                                                cGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGc                                                                                                                                                                                                >  2:2850/1‑289 (MQ=255)
                                                                                                       tGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCa                                                                                                                                                                                         <  2:46778/289‑1 (MQ=255)
                                                                                                       tGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCa                                                                                                                                                                                         >  1:46778/1‑289 (MQ=255)
                                                                                                       agggcAATCTGGAACAACCGTGGATTTGGTCTGGTAGGCCCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGAGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCa                                                                                                                                                                                         <  1:539210/286‑1 (MQ=255)
                                                                                                             aGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACt                                                                                                                                                                                   >  1:142199/1‑289 (MQ=255)
                                                                                                                   aacaacCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTa                                                                                                                                                                             >  2:59553/1‑289 (MQ=255)
                                                                                                                        ccggggaTTTGGTCCGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTGGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACctg                                                                                                                                                                        <  1:687708/283‑1 (MQ=255)
                                                                                                                             tATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCa                                                                                                                                                                   >  2:410366/1‑289 (MQ=255)
                                                                                                                                  ggTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCgg                                                                                                                                                                                                                <  1:124146/239‑1 (MQ=255)
                                                                                                                                  ggTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCgg                                                                                                                                                                                                                >  2:124146/1‑239 (MQ=255)
                                                                                                                                    tCTGATAGAACCTGGCTAAAAAATGGCGCAGGGCGCCGGTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACctgct                                                                                                                                                                      <  1:456806/279‑1 (MQ=255)
                                                                                                                                       ggTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTg                                                                                                                                                         >  1:486637/1‑289 (MQ=255)
                                                                                                                                          aGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGaa                                                                                                                                                       <  2:184990/288‑1 (MQ=255)
                                                                                                                                          aGCACCTGACTACAAAATGGCGCAGGTCGCCGTTGACCATAGTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTAATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGaa                                                                                                                                                       <  2:487602/288‑1 (MQ=255)
                                                                                                                                              ccTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTg                                                                                                                                                            >  2:456806/1‑279 (MQ=255)
                                                                                                                                                         aaTGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAGCGCCTTTAAAGATGCTAACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGAGGCATTGGGGATGCGCACGGTGTCACGCATGGCGACCGTAGCTACGTTTACCTCTACGTAAGTCAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACACACAAACCCTGCCATGTGGGGTCGTGGTGGTGGATACGCGCGACTACGGTAACATGCTGAAACTGGTGCTGTAcccactccacatgccg                                                                                                                                       >  1:526932/1‑274 (MQ=255)
                                                                                                                                                              cgcAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGaaaa                                                                                                                                  >  1:170727/1‑289 (MQ=255)
                                                                                                                                                                    tCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCgg                                                                                                                             <  2:378465/288‑1 (MQ=255)
                                                                                                                                                                                  atTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGc                                                                                                              <  2:454641/289‑1 (MQ=255)
                                                                                                                                                                                  atTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGAACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGc                                                                                                              <  2:76787/289‑1 (MQ=255)
                                                                                                                                                                                     ctcGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTATCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAg                                                                                                           <  2:596087/289‑1 (MQ=255)
                                                                                                                                                                                        ggTAGCGAAAACGCCTTTGAGGGTGCTGACCTGATGGCAAAGCTGAAACGGCTGGGGGTGGACGTGGGCGGTATTGGGGATGCGCACGGTCGCACGCCTGGCGAACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCGGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGca                                                                                                         <  1:679333/288‑1 (MQ=255)
                                                                                                                                                                                        ggTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTTCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGcac                                                                                                        >  1:580336/1‑289 (MQ=255)
                                                                                                                                                                                        ggTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGcac                                                                                                        >  2:431/1‑289 (MQ=255)
                                                                                                                                                                                            gCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTc                                                                                                    >  2:563404/1‑289 (MQ=255)
                                                                                                                                                                                             cGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGACAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCg                                                                                                   >  1:95984/1‑289 (MQ=255)
                                                                                                                                                                                             cGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGAt                                                                                                                   >  2:605462/1‑273 (MQ=255)
                                                                                                                                                                                             cGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGAt                                                                                                                   <  1:605462/273‑1 (MQ=255)
                                                                                                                                                                                             cGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCg                                                                                                   >  2:544998/1‑289 (MQ=255)
                                                                                                                                                                                                      cTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGc                                                                                          <  1:2850/289‑1 (MQ=255)
                                                                                                                                                                                                         tGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAg                                                                                       <  2:486637/289‑1 (MQ=255)
                                                                                                                                                                                                          gAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACCGGTTCTGAACGCTTACGAACTGCCGGAAAACCCGGTACCCCTGATCCTGCCACACCACCCGGGTAGCGGCAAg                                                                                       >  1:423471/1‑288 (MQ=255)
                                                                                                                                                                                                                 cTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAAc                                                                                                                                                                 <  1:592566/207‑1 (MQ=255)
                                                                                                                                                                                                                 cTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAAc                                                                                                                                                                 >  2:592566/1‑207 (MQ=255)
                                                                                                                                                                                                                  tgaccttgagcccAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGAACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTAt                                                                               <  1:59553/277‑1 (MQ=255)
                                                                                                                                                                                                                  tGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTAt                                                                               <  2:170727/288‑1 (MQ=255)
                                                                                                                                                                                                                  tGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGAACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTAt                                                                               <  2:205700/288‑1 (MQ=255)
                                                                                                                                                                                                                      cTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGAACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGgtg                                                                          <  2:142199/289‑1 (MQ=255)
                                                                                                                                                                                                                         aGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTTGCGGCAAGCCGTCTATCGGTGTTg                                                                       >  2:266158/1‑289 (MQ=255)
                                                                                                                                                                                                                         aGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACAAGCGACTACGGTCACCTGCTGCAACTGGTGCTGAACGCTATCGAACGGCCGGAAAACCGGGATTCCCTGCTCCTGCCAGCACACCCGGGTAGCGGCAAGCCGTCTACAGGTGTTg                                                                       >  1:616085/1‑289 (MQ=255)
                                                                                                                                                                                                                            gcCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATa                                                                    <  1:26271/289‑1 (MQ=255)
                                                                                                                                                                                                                              caagtcgAAACTGCGGGGCGTAGACGGAGGCGGTATTGGGGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATaaa                                                                  <  1:431/283‑1 (MQ=255)
                                                                                                                                                                                                                                aGCTGCAACTACTGGGCGTAGACGTTGGCGGTATTGGAGATGCGCCCGGTGGCACCCCTGTAGCACGTAGCTAAGTTTACCGCGACGAACGTAAAGAGAACTACAAACGCCTGATTGTCAGAGAAGACAACAAAACCCTGCTCGGGGCGGTACTGGTGGGCGATACTAGCGACCACGGTAACCTGATGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCCGATCCTGCCAGCACACACGGGTAGCGGAAAGCAGGTCATCGGTATTGATaa                                                                   >  1:214937/1‑286 (MQ=255)
                                                                                                                                                                                                                                aGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATaa                                                                   >  2:618125/1‑286 (MQ=255)
                                                                                                                                                                                                                                aGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATaa                                                                   >  2:474541/1‑286 (MQ=255)
                                                                                                                                                                                                                                           ctgGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAgcg                                                     >  1:160552/1‑289 (MQ=255)
                                                                                                                                                                                                                                           ctgGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAgc                                                      >  1:501080/1‑288 (MQ=255)
                                                                                                                                                                                                                                             gggCGGAGAGGTAGGGGGGATTGGGGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAgcgcg                                                   <  1:544998/289‑1 (MQ=255)
                                                                                                                                                                                                                                               gCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAgcgcgc                                                  <  1:563404/288‑1 (MQ=255)
                                                                                                                                                                                                                                                          ggCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCt                                         >  2:127594/1‑286 (MQ=255)
                                                                                                                                                                                                                                                               tATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCACGCAAATCTGCTCCTGCtt                                  >  2:119854/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                aTTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTc                                 >  1:363406/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                  tGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGAc                              <  1:474541/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tGGGGATGCGAACGGTCGAACGCCTGGCGAACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGAc                              <  1:301904/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCgg                                                                                                                             >  1:388460/1‑192 (MQ=255)
                                                                                                                                                                                                                                                                    gTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCgg                                                                                                                             <  2:388460/192‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGAAAGCGCGCACACCTGCCCCTGCTTCGACGTc                           >  1:167850/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                        tGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCAcc                        >  2:321399/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                         gcgcACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGaa                                                                                                                                                       >  1:86287/1‑161 (MQ=255)
                                                                                                                                                                                                                                                                         gcgcACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGaa                                                                                                                                                       <  2:86287/161‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gcgcACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCa                       >  2:53727/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                         gcgcACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCa                       >  1:432096/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                         gcgcACGGTCGCACGCCTGGCGAACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCa                       <  2:346924/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                            cACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAACGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTTCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCCCACCCGGGTAGCGGCAAGCCGTCTCACGGTTTTGATAAACTGCCGGCCAGCGCGAAAACCTGCTCCTGCTCCGACGTCACCAAAg                    >  1:681342/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                 tCGCAGCCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGa                                                          <  2:407586/246‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            ggCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATtgct      >  2:195222/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                              cgcACGTAACTACGCTTACCCCGACGAACGCAACGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAACACTCTCCCCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTCCTGAACGCCACCGAAATGCCGGAAAACCCGGATTCACTGACCCTGCCAGCACACGCGGGTAGCGGCAACCCGTCTACCGGAGTAGACAAACTGCCGGACCGCGACCAAATCTGCCCCTGCTTCCACCTTACCAAAGGTgtttcgcgcgcctc    >  1:475692/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                                               gcACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCa  <  1:321399/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                               
TAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCAGGTCGCCGTTGACCATATTCTCGGTAGCGAAAACGCCTTTGAAGGTGCTGACCTTAGCGCCAAGCTGAAACTGCTGGGCGTAGACGTAGGCGGTATTGGTGATGCGCACGGTCGCACGCCTGGCGCACGTAGCTACGTTTACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCA  >  NC_000913/3494749‑3495323

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