Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,780,060 G→A H889H (CAC→CAT)  ypjA ← adhesin‑like autotransporter YpjA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,780,0600GA96.0% 192.0 / ‑3.3 76H889H (CAC→CAT) ypjAadhesin‑like autotransporter YpjA
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (40/31);  total (45/31)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑GTGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCACACCCTCGACAATTTGTTCACCGCCGGTGAGCGTCGTGCCTGTCGCAGT  >  NC_000913/2779823‑2780302
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
gACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTg                                                                                                                                                                                                                                          >  1:802266/1‑250 (MQ=255)
gACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGACCCACCCTTTCCCACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCCTTTTCCAGCACCCTATTGGTCGCTAAATTTCCGGAAATGGCAAAAGAACCTTACTG‑ATGAGTACCGCtcg                                                                                                                                                                                                                                         >  2:48195/1‑249 (MQ=255)
          tcttctGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCCTTGCCCCCCCCTTTCCAACCGTGGCGTCGCC‑AGCTTCGGTACCTGCCAATACCAATAAATTACCGCCATTTTCCAGCCACATCTTGGTCGCTAAATTGCCGGAAATGGACAAAGTACCCTCCTGCGT‑AGAACCGCTGATTTCAATAcc                                                                                                                                                                                                                               >  2:67361/1‑251 (MQ=255)
           cttctGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGAG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCGGATTTCAATACCg                                                                                                                                                                                                                              >  1:84072/1‑251 (MQ=255)
                                             tGACCGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGt                                                                                                                                                                                           >  1:212838/1‑251 (MQ=255)
                                              gACCGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTG‑GCGGCGCCCTGACCCAGGATTTGCAATGCCCCACCCTTGCCAACCGGGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCCTTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGAACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTCGCCGTGCTCGTCGGGGGAGCG‑GCCCCGCCGtt                                                                                                                                                                                          >  2:733944/1‑251 (MQ=255)
                                                     cAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACCTATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACg                                                                                                                                                                                   >  2:316490/1‑251 (MQ=255)
                                                     cAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACg                                                                                                                                                                                   >  1:261320/1‑251 (MQ=255)
                                                      aaGGGTATATTGCCCGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑CGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACg                                                                                                                                                                                   >  1:528029/1‑250 (MQ=255)
                                                      aaGGGAATATTGCCCGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCACG‑AGCTTCGGTACCTGCTAAAACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTCAATTGCCGGAAATGGAAAAAGTACCATACTG‑ATGAGTACCGCTGATTTCAATACCGTTAACCGTGGTCGTCTTGAGTGAG‑GCACCGCGGTTCTGGAc                                                                                                                                                                                    >  2:247941/1‑249 (MQ=255)
                                                        gggTATATTGCCCGCCAGAGTGACCCTTTGTG‑GCGGAGTCCTGACCCGGGTTTTGCATTGCCCCACCCTTGCCAGCCGTGGAGTCGTG‑AGCTGCGGTCCCTGCTAATACAAAAAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGata                                                                                                                                                                                <  1:650525/251‑1 (MQ=255)
                                                            atatTGCCCGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGGGGGGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACAGATTGGTAGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑CTGAGTACCGCTGATTTCAATACCGTTAGCCGTACTCGTCTGGAGAGCG‑GGCCCCTTGTTCTGGACGTTATcat                                                                                                                                                                            >  2:229824/1‑249 (MQ=255)
                                                               ttGCCCGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCCGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGGGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTTCCGAAACATATTGGTCGCTAAATTGCCGGAAATGGTAAAAGTACCGTACTG‑GTGAGTACCGCTGATTTCAAt                                                                                                                                                                                                                                  >  2:1216/1‑195 (MQ=255)
                                                               ttGCCCGCCAATGTTAACCTTTGTG‑GCGGAGTCCTGCCCCAGGTTTTGCATTCCCCCCCCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTCCCTGCTAATCCCAATAAAGTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTc                                                                                                                                                                          <  1:824271/250‑1 (MQ=255)
                                                                   ccGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGCCCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGGGGAGTCGCG‑AGTTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCttt                                                                                                                                                                     <  1:612624/251‑1 (MQ=255)
                                                                   ccGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTAACGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTCATG‑ATTGGGACCACTGATTTCACTACCCTTAGCCGTGCTCCTCTGGAGA‑CGCGCACCGCTGTTCTGGACGATATCTGTCGCttt                                                                                                                                                                     >  2:28485/1‑251 (MQ=255)
                                                                     gCCAGAGTTAACTTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACCAATTGGTCGCTAAATTGCCGGAAATGGAAAAAGGACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCTCTTTGc                                                                                                                                                                   >  1:655321/1‑251 (MQ=255)
                                                                     gCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑CGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGc                                                                                                                                                                   >  1:332695/1‑251 (MQ=255)
                                                                      ccAGAGTTACCCTTT‑TGTGCGGAGGCCTGCCCTGGGTTTTGCCTTCCCCCCCCCGTGCCAACCGGGGTGGCGCG‑AGCTTAGGTTCCTATAAAAAACAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCTCTTTGcc                                                                                                                                                                  <  2:858581/251‑1 (MQ=255)
                                                                       cAGAGTTAAACTTTGTG‑GCGGAGTCCTGACCCAGGGTTTGCATTGCCCCACACTAGACAAACGTGGAGTCGCG‑AGCTTCGGAACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATAACGTTAGCCGAGCTCGTCTGTAGCGAG‑GCACCGCTGGTCTGGACGATATCTTACGCTTTGcc                                                                                                                                                                  >  1:775312/1‑250 (MQ=255)
                                                                         gagTTACCCTTTGGG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCCCCCTTGCCAACCGTGGAGCCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCAtt                                                                                                                                                               <  1:273206/251‑1 (MQ=255)
                                                                         gagTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCCACCGTGGAGTCGCG‑AGCTTCGGGACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGCCGCTAAATTGCCGGTAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGCTTTCAATACAGTTAGCCGTGCTCGTCTGGAGGGCG‑GCACCGATGTTCTGGAAGATATCTGTCGCTTTGCCAtt                                                                                                                                                               >  2:782771/1‑251 (MQ=255)
                                                                             tAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCCACCGTGGCGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCCCCGCTGATCTGGGCGATATCTGTCGATTTGCCATTCTCg                                                                                                                                                           >  2:284227/1‑251 (MQ=255)
                                                                             tAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTAACGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCg                                                                                                                                                           >  1:116573/1‑251 (MQ=255)
                                                                               cccTTTTTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCCCCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTCCCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGtt                                                                                                                                                         <  1:768156/250‑1 (MQ=255)
                                                                                    tgtg‑GCGGAGGCCTGCTCCAGGTTTTGCCTTGCCCACCCCTTGCCACCCGTGGAGTCGCG‑AGCCTCGGTACCTGCTAATCCCAATAAATTGCCGCCATTTTCCAGCAACCTATGGGCCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTTATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGATGATATCTGTCGCTTTGCCATTATCGTTGACt                                                                                                                                                     <  1:161781/250‑1 (MQ=255)
                                                                                       g‑gCGGAGTCCGGCCCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTCCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACAGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATATGTCGCTTTGCCATTATCGTTGACTGTCa                                                                                                                                                 <  1:458556/251‑1 (MQ=255)
                                                                                         gCGGATTCCTGCCCCAGGTTTTGCATTGCCCCCCCCTTGCCAACCGTGGAGTCGCG‑AGTTTCGGTACCTGCTAATCCCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAg                                                                                                                                                <  1:408955/251‑1 (MQ=255)
                                                                                                cctaccccAGGTTTTGCTTTGCCCCCCCCTGCCCAACGGGGGTG‑CGCGGAGCTTCGGTCCCTGCTAAACCCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTGCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTAcc                                                                                                                                          <  2:318885/245‑1 (MQ=255)
                                                                                                  tGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGCCGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGcc                                                                                                                                       >  1:224134/1‑251 (MQ=255)
                                                                                                    ccccAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTCCCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCtt                                                                                                                                     <  1:153335/250‑1 (MQ=255)
                                                                                                       cAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑CGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGGGCG‑GCCCCGCAGTTCTGGACGATATCTGTCGCGTTGCCATTATCGTTGACTGTCGGCGTACCCCCTTCAt                                                                                                                                  >  2:564488/1‑251 (MQ=255)
                                                                                                          gTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTc                                                                                                                                    >  1:529552/1‑246 (MQ=255)
                                                                                                          gTTTTGCATTCCCCCCCCCTTGCAAACCGGGGGGGCGCG‑GCCTTCGTTCCTTGTTAATCCAAAAAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTc                                                                                                                                    <  2:529552/246‑1 (MQ=255)
                                                                                                           ttttGCATTGCCCCAACCTTGCCAACCGTGGAGTCCCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑GTGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGAt                                                                                                                              >  1:184731/1‑251 (MQ=255)
                                                                                                            tttGCATTGCCCCACCCTCCCCAACGGTGGAGCCGCG‑AGTTTCGGTACCTTCTAATCCCAATAAATGACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGACTTCATTGAt                                                                                                                              <  1:302101/250‑1 (MQ=255)
                                                                                                             ttGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGACATTTTCCATCCACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑CTGATTACCGCTGATTTCCATACCGTTAGCCGTGCTCGTCTGGAGTGCG‑GCAACGCTGATCTGGACGATATCTGTCGCTTTGCCCTTATCTTTGACTGTCCGCGGCCCGCCTTCCTTGATCt                                                                                                                            >  2:253829/1‑251 (MQ=255)
                                                                                                                     ccccACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑Gttttt                                                                                                                   >  1:375513/1‑251 (MQ=255)
                                                                                                                      cccACCCTTGCCAACCATGGAGTCGCG‑AGCTTCGGTCCCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGTAAATGGAAAAAGTACCGTAATA‑ATGAGTACCGCTGATTTCAATACCGTTAGCCTGGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGCCTGTAACTGTCGCTTTGCCATTATCGTTGACTGTAAGCGGACCGCCTTCATTGAGCTTTTGtttt                                                                                                                    >  2:692296/1‑250 (MQ=255)
                                                                                                                           ccTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGcc                                                                                                             >  1:240391/1‑251 (MQ=255)
                                                                                                                              tGCCAACCGTGGAGTCGCG‑AGTTGCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCct                                                                                                            <  1:493116/249‑1 (MQ=255)
                                                                                                                                        ggggccgcg‑AGCTCCGGTACCTGTTAATACCAAAAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTg                                                                                                <  1:794912/246‑1 (MQ=255)
                                                                                                                                         ggggtcgctagtgtcggctccTGCTAATACCAACAATTCCCCGCCATTTTCCAGCAACATATTGGGCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGGACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTg                                                                                                <  2:460619/232‑1 (MQ=255)
                                                                                                                                                             cTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATCGCCGGAAATGGAAAAAGTACCGTACTG‑ATACGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTccgccg                                                                            >  1:833022/1‑251 (MQ=255)
                                                                                                                                                               gCTAATCCCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTcc                                                                                                        <  1:422117/221‑1 (MQ=255)
                                                                                                                                                               gCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTcc                                                                                                        >  2:422117/1‑221 (MQ=255)
                                                                                                                                                                   aTACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCCTTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGtt                                                                       >  2:185660/1‑250 (MQ=255)
                                                                                                                                                                     aCCAATAAATTCCCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGATCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAAt                                                                    <  1:760217/251‑1 (MQ=255)
                                                                                                                                                                     aCCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCCTTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAACCTGTTTGTATTCCGCCGTCGTTAAt                                                                    >  2:213081/1‑251 (MQ=255)
                                                                                                                                                                      ccAATAAATTCCCGCCATTTTCCAGCAACCTATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTCACTGTCAGCGTACCGCCTTCCTTGATCTTT‑GTTTTTTTTGACTCCACGTTAGCTGATACTATTTTTATTCCGCCATCGTTAAtt                                                                   >  2:157937/1‑251 (MQ=255)
                                                                                                                                                                        aaaaaaTTACCGCCATTTTCCAACAATATATTTGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTCCCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAAttgt                                                                 <  2:228255/248‑1 (MQ=255)
                                                                                                                                                                         aTAAATTTCCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAAttgtt                                                                <  1:599937/251‑1 (MQ=255)
                                                                                                                                                                           aaaTTACCGCCATTTTCCATCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCTCCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGGAACTGTTTGTATTCCGCCGTCGTTAATtgttgt                                                              <  2:333810/251‑1 (MQ=255)
                                                                                                                                                                           aaaTTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATtgttgt                                                              <  1:644936/251‑1 (MQ=255)
                                                                                                                                                                                 ccGCCATTTTTCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCTTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCAtt                                                        <  2:391052/251‑1 (MQ=255)
                                                                                                                                                                                  cGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTc                                                       >  1:212120/1‑251 (MQ=255)
                                                                                                                                                                                      attttccagcactagatggttCGCTTAATTGCCGGAAATGGAAAAAGTACCTTATTG‑ATGAGTCCCGCTGAGTTCAATACCTTTAGCCTTGCTCGTCTGGAGAGCG‑TCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTGGACTGTCAGTGTACCGCCTTCATTGATCGTT‑TTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCca                                                   <  2:475591/231‑1 (MQ=255)
                                                                                                                                                                                      aTTTGCCAGCAACATATTGGTCGCTAAATTGCTGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCTCCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACGGTTTGTATTCCGCCGTCGTTAATTGTTGGCTCATTCGCca                                                   <  2:240391/251‑1 (MQ=255)
                                                                                                                                                                                       ttttCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCcac                                                  >  2:108090/1‑251 (MQ=255)
                                                                                                                                                                                           ccAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAACCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCATCGTACCGCTTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGACACAccc                                              >  1:863804/1‑251 (MQ=255)
                                                                                                                                                                                               cAACATTTTGGTCGCTAAATTGCCGGAATTGGAAACAGTACCGTACTG‑ATGAGTACCGCTTATTTCAATCCCGTTAGCCGTGCTCGTCTGGATAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCACACCCTCGa                                          <  2:261320/251‑1 (MQ=255)
                                                                                                                                                                                                aaCATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGCAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCACACCCTCGAc                                         >  1:749219/1‑251 (MQ=255)
                                                                                                                                                                                                            gCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCACACCCTCGACAATTTGTTCAcc                             >  1:843430/1‑251 (MQ=255)
                                                                                                                                                                                                            ccgaaagtGCCGAAAATTGAAAAAGTACCGTACTG‑ATGAGGACCGCTGAGTTCAATACCGTTATCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCACACCCTCGACAATTTGTTCAcc                             <  2:166820/244‑1 (MQ=255)
                                                                                                                                                                                                             cTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAACCTGTTTGTATTCCGCCGTCGTTAATTGTTGTCGCATTCGCCACACACTCGACAATTTATTCAccg                            >  2:607856/1‑251 (MQ=255)
                                                                                                                                                                                                               aaaTTGCCGGAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTGAATTGTTGTCTCAGTCGCCACACCCTCGACAATTTGTTCAccgcc                          <  1:229824/251‑1 (MQ=255)
                                                                                                                                                                                                                     caggAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCACACCCTCGACCATTTGTTCACCGCCGGTgcg                    >  2:478851/3‑249 (MQ=255)
                                                                                                                                                                                                                       ggAAATGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGAG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCACACCCTCGACAATTTGTTCACCGCCGGTGAGcg                  >  1:33974/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaaTGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCACACCCTCGACAATTTGTTCACCGCCGGTGAGcgtc                >  1:447744/1‑251 (MQ=255)
                                                                                                                                                                                                                          aaTGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCACACCCTCGACAATTTGTTCACCGCCGGTGAGcgtcg               >  1:146737/1‑251 (MQ=255)
                                                                                                                                                                                                                           aTGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGCAACTGTTTg                                                                                      >  2:287111/1‑179 (MQ=255)
                                                                                                                                                                                                                           aTGGAAAAAGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTg                                                                                      <  1:287111/179‑1 (MQ=255)
                                                                                                                                                                                                                               aaaaaGTACCGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCACACCCTCGACAATTTGTTCACCGCCGGTGAGCGTCGTGCCt          >  1:277867/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaGTACCGTACTG‑ATGGGTCCCGCTGATTTCAAGCCCGTTAGCCGGGCTCGTCTGGAGAGCG‑CCCCCGCTGTTCTGGACGATATCTGCCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCTTTAATTGTTGTCTCATTCGCCACCCCCTCGACAATTTGTTCACCGCCGGTGAGCGTCGTGCCTGTc       <  1:316490/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aaGTACCGTACTG‑ATGATTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGGGCG‑GCACCGCTGTTCTGGACGATATCTGTCGTTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCACCCCCTCGACAATTTGTTCACCGCCGGTGAGCGTCGTGCCTGTc       <  1:108090/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ccGTACTG‑ATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCACACCCTCGACAATTTGTTCACCGCCGGTGAGCGTCGTGCCTGTCGCAGt  >  1:230757/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
GACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTG‑GCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCG‑AGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTG‑GTGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCG‑GCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTT‑GTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCACACCCTCGACAATTTGTTCACCGCCGGTGAGCGTCGTGCCTGTCGCAGT  >  NC_000913/2779823‑2780302

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