Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,287,009 A→T 7.1% intergenic (+483/+57) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,0090AT 7.1% 176.4 / 3.2 66intergenic (+483/+57)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (19/40);  new base T (2/3);  total (22/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

AAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGC‑TGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGAT  >  NC_000913/1286722‑1287290
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                         
aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGC‑TGGCGCTCAAGTCGACCATTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCa                                                                                                                                                                                                                                                                                           <  2:583613/288‑1 (MQ=255)
    cgtACGCAGCGCCGCGCGCTGCTCCTTGTGGGTCGTTGCTGTGGGCAACGCTCGATGGT‑GGGCGCTCAAGTCGAACCTTGGTGGAAGCTTCTCAGCGTTACCCGTTTGGGCAGAGTGTGTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATggg                                                                                                                                                                                                                                                                                       <  1:458022/286‑1 (MQ=14)
    atctctgaGCGCTGGGCGTTTGTGCCTGCGGGGCGGTGCGTGCGGCAACGCTCTCTCGC‑TGGCGCTCCAGCCGAACCGTGGTCGAAGGTTCTCAGCCGTCCCCGCTTGCGCAGAATATTTGATTGCGGAGTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATggg                                                                                                                                                                                                                                                                                       <  1:115411/281‑1 (MQ=14)
    aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTGCTTGGCTGCGGCAACGCTCTCTCGC‑TGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATggg                                                                                                                                                                                                                                                                                       <  2:227861/288‑1 (MQ=255)
    aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGTAACGCTCTCTCGC‑TGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATggg                                                                                                                                                                                                                                                                                       <  2:925418/288‑1 (MQ=255)
     ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGC‑TGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCa                                                                                                                                                                                                                                                                                     <  2:417165/289‑1 (MQ=255)
     ttACTCAGCGCGGCGCGCTTCGCCCTTCGGGGCGTTGCCTGCGGAAACGCTCTCTCGC‑TGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCa                                                                                                                                                                                                                                                                                     <  2:784104/289‑1 (MQ=255)
                  cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGC‑TGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGAtt                                                                                                                                                                                                                                                                         >  2:228645/1‑288 (MQ=255)
                       ttggacctgtggtgcGTTGCCTCGGGGAAGGCTCTC‑CGCGGGGTGCGCAAGTTGAAAGTTGAGCGAAGTGTATCACCCGTCCCCGTTTGGGCATAATATTTGATTGCGTAGTCGT‑TTGAGAATTCGGTGGTTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAt                                                                                                                                                                                                                                                                   <  1:741126/275‑1 (MQ=14)
                               tCGTGTCGTTGCCTGCGGCAACGCTCTCTCGC‑TGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTC‑ATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCgttggg                                                                                                                                                                                                                                                           <  2:839444/289‑6 (MQ=11)
                               tCGGGTCGTGGGGTTGGGCAACGCTACCGCGC‑CGGGCGTGAAGTAGAACCTTGGTCGAGGCGTCTACTCCTTCGCGTATAGGGGAGAAGATTTGATTGCGGAGTCGT‑TTCAGATTTCGGGGGCTTTTGAAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                                                                                                                                                                                                           <  1:683961/289‑1 (MQ=14)
                                gggttCGTTGCCTGCGGAAACCCCCTCTCGC‑TTGCGCTCAAGAGGAACCGTGGTCGATGGTTCGCATCCGTCACCGCTTGGGCAGAATATTAGATTGGGGAGTCGT‑TTCAGGATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                                                                                                                                                                                                           <  1:564202/284‑1 (MQ=14)
                                        tGCCTGCGGCAACGCTCTCTCGC‑TGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGTATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATTGGCAGAATATTTTATTCCGGATTCGCTTGAGAGTTCAgg                                                                                                                                                                                                                                                  >  2:740607/1‑289 (MQ=255)
                                            tGCGGCAACGCTCTCTCGC‑TGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGCCt                                                                                                                                                                                                                                              >  2:775626/1‑289 (MQ=255)
                                                       tctctcGC‑TGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                                                                                                                                                                                                                   <  2:249475/289‑1 (MQ=255)
                                                        cttcgctttggcgctgcgacctacgcgtGGGCGAAGCTTCTCAGCCTTACCCGCATTGGCATCATATTATATTGGCGCTTCGT‑TTGAGGGTTCAGGTACTTTTGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                                                                                                                                                                                                                                   <  1:482265/262‑1 (MQ=14)
                                                           tcGC‑TGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggt                                                                                                                                                                                                                               <  2:50955/289‑1 (MQ=255)
                                                                tgggtcTCAAGGCGAACCGTGGTGGAAGTTTCTGATCCTTACCCGCTGGGGCAGAGTACATGATTGGGGAGTCGT‑TTGAGGTTTCCGGTGTTTTTGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                                                                                                                                                                                                           <  1:748050/284‑1 (MQ=14)
                                                                tGGCGCTCCAGTCGAACCTTGGTCGGAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                                                                                                                                                                                            <  2:386162/288‑1 (MQ=255)
                                                                    gcTCAAGTCGAACCTTGGTCGAAGCTTCTCATGCTTCCCCGCTTGGGCAGAATAGTTGATTGCGGATTCGT‑TTGAGAATTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                                                                                                       <  2:77500/289‑1 (MQ=255)
                                                                    gcTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                                                                                                       <  1:228645/289‑1 (MQ=255)
                                                                     ctcaggcctaaactgggTGGAAGTTTTGCATCCTTCCACGCTTGGGCATAAGATTTGATTGCGGATTCGT‑TTGAGAATTCCGTGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTATC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                                                                                                       <  1:706519/273‑1 (MQ=14)
                                                                     ctcaggacgaacctcgtttgggtgttCTCGTCCTTCGTGGATAGGGGAGAATGTTAGATTAGGGATTCGT‑ATGAGGATTCCGGGGGTTTTGGAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                                                                                                       <  1:455599/264‑1 (MQ=14)
                                                                     cTCAAGTCGAACCTTTGTTGGAGCTTCTCATCCTTCCCCGCTTTGGCAGGATAATTGGTTGCGGATTCGT‑TTGAGAATTCCGGGGCGTTTGAAAGTGATGGTTGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                                                                                                       <  2:255076/288‑1 (MQ=255)
                                                                        aaGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCATGCCTCTTCTATACATCTTACGGTTCGTTCCATGTGCAGACTCTCTCTCCCGGCGGTCCGAGCCGAACTTTGTTCAAAGC‑TCTCAATCCTTCCCCGCCTGGTCAGATTTTTGAATCACGGTTTCGCTGGTGATTTCCGGGCCTTTTGAAGGTGATGGTGTTGGGGGAAgc                                                                                                                                                                                                                   >  1:255076/1‑288 (MQ=2)
                                                                         aGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAAAAAATTGGATTGCGGATTCGT‑TTAGGAATTCCGGGGCTTTTGAGAGAGATGGTGGTGGGAGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                                                                                                                                                  <  2:233731/289‑1 (MQ=255)
                                                                                            cTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCt                                                                                                                                                                                               >  2:864783/1‑289 (MQ=255)
                                                                                                     cgtaatcgtttgggcagcggagtggATTGCGGGGTAGT‑GTGAGAATTTCGGTGGTTTTGTAAGAGATGGTGGTGGGGGAAGGATTAGTCTGGGCTGCGTGCTTCGCCCTGCGGGTCGTTGCCTGCGGCAACGCTCGCGCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                                                                                                                                                                      <  1:915816/265‑1 (MQ=14)
                                                                                                       tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGT‑TTCAGACTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGCAGGATTACTCACCGCGCCGCGCCTCGCGCTTAGGGTCGTGTCCCCCGTCTACGCTCGCCCGGCGGCGCCCCACCCAAACCCCGCTCGACCCTCCTC‑ATCTTTGCCCGCTTCGGCAGAGTCGTGGGTTCCGGTCTCGCTTCACAGTTCAGGGACTGTCGCATGTGAAGGTTGTTGCGGTAGGATTACttctccctcagcgtctcgccgctcga                                                                                                                                                                                    >  1:816004/1‑264 (MQ=2)
                                                                                                        ccccGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTTggg                                                                                                                                                                                   >  2:788664/1‑289 (MQ=255)
                                                                                                                  gCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCt                                                                                                                                                                          >  2:144434/1‑288 (MQ=255)
                                                                                                                        tatTGGATTCCGGATGCTTGTTGA‑AATTCTGGGTGTTTTGAAAGTGTTGGTGGTGGGGGAAGGATTCCTCATCGTTGCGCTCTGCGCCCTTCGGGTCGTGGCCTACGCCAACGCTGTATCGCTGACGGTGGAGCCGGACCTTAAGCGAAACTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                                                                                                                                                   <  1:788664/289‑1 (MQ=14)
                                                                                                                           ttGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACg                                                                                                                                                                 >  2:593051/1‑288 (MQ=255)
                                                                                                                                tGCGGATTCGT‑TTTAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCCCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTCCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGTTCGTTGCTCGCGGCAACGctctc                                                                                                                                                            >  2:281558/1‑288 (MQ=255)
                                                                                                                                  cGGATTCGT‑TTGAGAACTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCcctctc                                                                                                                                                          >  2:673368/1‑288 (MQ=255)
                                                                                                                                           tttggaaatgtgtggtcTTTTGAAAGTGTTGGTGGTGGGGGAAGTATTACTCGGCGGTGCTGGCTTCGTCCTTCGGGTCGTTGCCTGCGGCAACGCTCTATCGCTGACGCTCGAGACGATCCTTAGTGGATGCGTGTC‑ATCGTTACCCGTGTGGGCAGAATATTTGATTGGGGATTCGTTTGAGATTTCAGGTACTTTTGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                                                 <  1:90371/273‑1 (MQ=14)
                                                                                                                                           attcagtgtttcggggcggctggatgggATGGGTGTGGGGGAAGGATTAGCCTCGGCGGCGTTCTTCGCCCTTCGGGTCGTGGCCTGCGGCGTCGCTCGCGCGCCGCCGCTCGGGCCGAACCTTAGTCGAAGCGTATC‑ATCCTTACCCGTATGGGCAGAAGATTTGATTGGGGAGTCGTTTGAGATTTCGGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                                                 <  1:653183/262‑1 (MQ=11)
                                                                                                                                            ttgagtgtttggtgggtgtatgaggTGATAGTGGTGGGCGAAGGATTTCTCAGCGCCGTGCGCTTCGCCGGTCGGGTCGTGGCCTGCGGTAACGCTGTCGCGCTGCCGATCCGGCCGAACCTTAGTCGAAGCGTCTG‑ATCCTTCCCCGCATGGGCAGAATGTTTGATTGCGGATTCGTTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                                                 <  1:609273/264‑1 (MQ=11)
                                                                                                                                                               tGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTTCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGaa                                                                                                                             <  2:549568/289‑1 (MQ=255)
                                                                                                                                                                              ggggggAAGGTTTACTCGCGGCGGTACGCTTCGCCCTTTGGGTCGTTGCCAGCGGCAACGCTCTCGCGATGTCGCCCGTGCATAACGTTGGTCGCGTCTGATG‑TGTCTAAGCGGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTc                                                                                                               <  1:593051/287‑1 (MQ=11)
                                                                                                                                                                                        aTTAGTCATCGGTTCGCGCTTCGCTGGTCGGGTCGTTGTCTTCGGAAGCTCTCGCTCGCTGACGCTCGAGCCGAATCTTAGTCGACGGTTGTG‑TTCCTTCCCTGTATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                     <  1:673368/288‑1 (MQ=14)
                                                                                                                                                                                         ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                                    <  2:738709/288‑1 (MQ=255)
                                                                                                                                                                                          tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTTCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAATGTTTACTCTGCGCTGCTTGCTTCGCCCTTCGGATCGTTGCTTGCGGCACCCCTCCCTCGCTTGCACTCGACTCGAACCTTTCTCGTAGATTCGCATCCTTCCCTGCTTgcgtag                                                                                                  >  1:737452/1‑284 (MQ=14)
                                                                                                                                                                                                 cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                                                                                                                                                                                                           <  1:767259/129‑1 (MQ=255)
                                                                                                                                                                                                 cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTCAGCGCTGCTCATTTCTCTCTTAGGGTCGCTGCCTGCGGCAACGCTTTCTCGCTGGCGCTCGACTCGAACCTTGGGCGAAGCTTCTCATCTCTACCCGCCTGGGACGCATAttt                                                                                           >  1:408084/1‑289 (MQ=14)
                                                                                                                                                                                                      cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTCAtt                                                                                       >  2:468548/1‑288 (MQ=255)
                                                                                                                                                                                                           ttggctcttgggggtgtgtccTCCGGCAACTTTCACTCGCTGCCTCTCGTGCCGAACGTTAGTTGAGGCTTGTA‑ATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAt                                                                                 <  1:144434/270‑1 (MQ=14)
                                                                                                                                                                                                                ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                                        <  2:755951/261‑1 (MQ=255)
                                                                                                                                                                                                                  ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGCATATTTTATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAACGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGCGTCGAACATTGGTCGAAGCTTCTCATCCTTCCCCTCATGTGCAGCATATTTAATTGCGTAGTCGTCgg                                                                          >  2:838209/1‑289 (MQ=14)
                                                                                                                                                                                                                          gTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGAATCCGCTTGAGCGTTCAGGAACTTTTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCACAGCTCTCTCGCTGGCCCTCGAGTCGAACCTTGGCTGAAGCTTCCACTCCTTCACCGACTGGCCAGAataactgtctctgatacactactgactgcg                                                                   >  1:755951/1‑261 (MQ=14)
                                                                                                                                                                                                                               cTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGAc                                                             <  2:710946/289‑1 (MQ=255)
                                                                                                                                                                                                                                  cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCCGTGCTGCGCGCTTCGCCCTTCGGTTCGTTGCCTGCGGCAACGCCCTCTCGCTGGCGCCCCAGTCGAACCTTGGTCGACGCTTCTCATCCTTCCCCGCATGGGTAACAGATTTACATTCGTCTCCGTTGGGCAGTTCCggcatgtt                                                          >  1:297593/1‑283 (MQ=14)
                                                                                                                                                                                                                                             tctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtg                                               <  2:737452/289‑1 (MQ=255)
                                                                                                                                                                                                                                                        cTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCtt                                                                                                                         <  2:848056/204‑1 (MQ=255)
                                                                                                                                                                                                                                                         tCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGa                                    <  2:408084/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                gAACCTTAGTCTAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATAGTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCg                            <  2:835386/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                  acctgggTCGAAGTTTCTG‑ATCCTTCCCCGTATGGGCAGAATATTTGATTTCGGCTTCGCTGGAGAGTTCAGGTACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa                           <  2:828375/282‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aCCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCATGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaa                                                      <  2:482433/261‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aCCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTc                                                                                                               >  1:848056/1‑204 (MQ=255)
                                                                                                                                                                                                                                                                           gAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTTGGCAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGGTTCGAACCTTCGAAg                 >  2:597841/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                              gCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCTCATTGGCAGAATATTTAATTGCGGCTTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGCTGTTGGTGGGGCAAGGATTCGACCCTTCGAAGTCg              >  2:864938/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                tgctc‑ATCCTTCCCCGCATGGGCAGAATATTGGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGAt            <  2:316192/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        ccTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCa     <  2:255237/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        ccTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCa     <  2:295744/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        ccTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCGTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGTGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCa     <  2:797279/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGa   >  2:854989/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                           tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCACGCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGGCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATTGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTTAAAGTGGTGGGGGTGGGGGAAGGATTCGacactccaaatacgcagaaggcaggt  >  1:255237/1‑264 (MQ=11)
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                         
AAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGC‑TGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGT‑TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTC‑ATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGAT  >  NC_000913/1286722‑1287290

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