Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,576 T→C 10.4% intergenic (+21/+220) ybiT → / ← ybiU ABC‑F family putative regulatory ATPase/DUF1479 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913857,5760TC10.4% 429.7 / 30.8 147intergenic (+21/+220)ybiT/ybiUABC‑F family putative regulatory ATPase/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base T (61/69);  new base C (6/9);  total (68/79)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-04

Only 100 of 148 total aligned reads displayed.

CTGTTTGGTAAGTTAATGATGCAGAAGCCGAACATTCTGATCATGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTA  >  NC_000913/857304‑857850
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                  
cTGTTTGGTAAGTTAATGATGCAGAAGCCGAACATTATGATCATGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGAACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTGACACGACCGTGAGTTCGTAAGCGCCCTGGCGACACGCATTCTGGAAATCAACCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATa                                                                                                                                                                                                                                                                                   <  2:516017/274‑1 (MQ=255)
           gTTAATGATGCAGAAGCCGAACATTCTGATCATGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGTTTACCTGCGTAGTAAAGGGAACGAGTAAGTGATTATCCGTCACCCCGAAACTGGCgggggg                                                                                                                                                                                                                                                                        >  4:680078/1‑272 (MQ=255)
               atgatgCAGAAGCCGAACATTCTGATCATGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcgg                                                                                                                                                                                                                                                                    <  2:680080/274‑1 (MQ=255)
                                          aTGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGCTGACAACTGCAAAATAGCCTg                                                                                                                                                                                                                                        >  2:294533/1‑275 (MQ=255)
                                          aTGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAAGGCAAAATAGCCTg                                                                                                                                                                                                                                        >  2:73018/1‑275 (MQ=255)
                                               cGAACCGACCAATCACCTGTGTATGGAATCCATTACGTCGCTGAACATGGCACTGTAACTGTAGCAGGGCACGCTGATATTTGTTTCACAAGACCGTGAGTGCGTAAGCTCCCTGGAGACCCGCATTCTAGACATCACCCCGGAACGCGTGATAGACTTTAGCGGTAAATACGAAGATTACCAGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGACACCCCGATACTGTCGGGGTGACGGCGGTTGATGACAAATGCAAAATAGACTGATgcg                                                                                                                                                                                                                                   <  4:294532/275‑1 (MQ=255)
                                                    cGACCAACCCCCTGGATATGGAAGCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACg                                                                                                                                                                                                                              <  2:108299/275‑1 (MQ=255)
                                                                 ggATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGCGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTATTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTCACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGTCAGGCTTATCAGGCCTa                                                                                                                                                                                                                 >  2:710770/1‑275 (MQ=255)
                                                                       ggAATCCATTGAGGCGCTGAACATGGCACTGGAACTGTTTCAGGGGACGCTGATATTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCAGTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGt                                                                                                                                                                                                           <  4:677705/275‑1 (MQ=255)
                                                                                tgtgTCGTTGAACGTGGGACTGGTACGGTATCAGGGCACCCTGATCTTTGTTTAACACGACCGTGAGTGCGTAAGGTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTAGCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCaa                                                                                                                                                                                                   <  4:73017/271‑1 (MQ=255)
                                                                                                                aGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTACGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGAt                                                                                                                                                                  <  1:22340/275‑1 (MQ=255)
                                                                                                                 gggCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtga                                                                                                                                                                                                                                                             >  2:215329/1‑183 (MQ=255)
                                                                                                                     aCGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGAc                                                                                                                                                                                                                                                         <  4:215328/183‑1 (MQ=255)
                                                                                                                      cGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAcc                                                                                                                                                                                                        >  2:413773/1‑231 (MQ=255)
                                                                                                                      cGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGANTATCCGTCACCCCGATCATGTCGGGGTGACGGCGGAGGATGACAAATGAAAAATAGCCTGATGCGTTACGCTTATCAGGCNAACCCAGTACCTCACATACATTGAANTTTCACGCATTTGTATGCCTGATaaatt                                                                                                                                                             >  3:144506/1‑271 (MQ=255)
                                                                                                                          gATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCa                                                                                                                                                                                                    <  4:413771/231‑1 (MQ=255)
                                                                                                                             cTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATa                                                                                                                                                                                                                                                                                   >  8:672019/1‑149 (MQ=255)
                                                                                                                              tttGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAc                                                                                                                                                                                                                                                                                  >  8:358633/1‑149 (MQ=255)
                                                                                                                               ttGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATAc                                                                                                                                                                                                                                                                                  >  5:261414/1‑148 (MQ=255)
                                                                                                                                    tCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgg                                                                                                                                                                                                                                                                            <  8:35252/149‑1 (MQ=255)
                                                                                                                                    tCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCg                                                                                                                                                                                                                                                                             <  5:256836/148‑1 (MQ=255)
                                                                                                                                     cacaCGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCggg                                                                                                                                                                                                                                                                           <  8:6207/149‑1 (MQ=255)
                                                                                                                                     cacaCGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgg                                                                                                                                                                                                                                                                            <  5:697264/148‑1 (MQ=255)
                                                                                                                                     cacaCGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgg                                                                                                                                                                                                                                                                            <  5:597593/148‑1 (MQ=255)
                                                                                                                                       cacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCgggg                                                                                                                                                                                                                                                                          <  6:495489/148‑1 (MQ=255)
                                                                                                                                       cacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGt                                                                                                                                                                                                                                                                         <  7:603073/149‑1 (MQ=255)
                                                                                                                                         cGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTg                                                                                                                                                                                                                                                                        <  6:542551/148‑1 (MQ=255)
                                                                                                                                          gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAc                                                                                                                                                                                                                                                                      <  7:127862/149‑1 (MQ=255)
                                                                                                                                            ccGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTg                                                                                                                                                                                                                                                                        >  2:376096/1‑145 (MQ=255)
                                                                                                                                            ccGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcg                                                                                                                                                                                                                                                                     <  6:779844/148‑1 (MQ=255)
                                                                                                                                               tGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCgggg                                                                                                                                                                                                                                                                          <  6:169907/140‑1 (MQ=255)
                                                                                                                                               tGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCgggg                                                                                                                                                                                                                                                                          >  5:169907/1‑140 (MQ=255)
                                                                                                                                                   ttCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTGCGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtg                                                                                                                                                                                                                                                              >  6:261609/1‑148 (MQ=255)
                                                                                                                                                   ttCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtg                                                                                                                                                                                                                                                              <  6:261414/148‑1 (MQ=255)
                                                                                                                                                    tCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtga                                                                                                                                                                                                                                                             <  5:88544/148‑1 (MQ=255)
                                                                                                                                                              ccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTGTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACGGTCGGGGTGACGGCGGCTGATGACAAATGc                                                                                                                                                                                                                                                   >  5:591129/1‑148 (MQ=255)
                                                                                                                                                                tGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCaa                                                                                                                  >  2:346841/1‑275 (MQ=255)
                                                                                                                                                                  gCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCGGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAAt                                                                                                                 <  2:154834/274‑1 (MQ=255)
                                                                                                                                                                       ccGCATTCTGGAAATCACCCCGGATCGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGa                                                                                                                                                                                                                                                                       >  7:764252/1‑119 (MQ=255)
                                                                                                                                                                       ccGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCt                                                                                                                                                                                                                                         >  8:486276/1‑149 (MQ=255)
                                                                                                                                                                        cGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTg                                                                                                                                                                                                                                        <  7:486662/149‑1 (MQ=255)
                                                                                                                                                                                     taaCCCCGGAACGGGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATGCTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGAt                                                                                             <  1:144507/273‑1 (MQ=255)
                                                                                                                                                                                         cccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGc                                                                                                                                                                                                                             <  7:850670/143‑1 (MQ=255)
                                                                                                                                                                                            ggAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGc                                                                                                                                                                                                                    <  8:822234/149‑1 (MQ=255)
                                                                                                                                                                                                cgggTGATGGAATTTTGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGGGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATcc                                                                                  <  4:196897/272‑1 (MQ=255)
                                                                                                                                                                                                cgcgTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTa                                                                                                                                                                                                                 >  6:210607/1‑148 (MQ=255)
                                                                                                                                                                                                   gTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCa                                                                                                                                                                                                             >  7:133256/1‑149 (MQ=255)
                                                                                                                                                                                                       tCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTa                                                                                                                                                                                                          <  5:763107/148‑1 (MQ=255)
                                                                                                                                                                                                          aCTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTa                                                                                                                                                                                                                 <  5:128734/138‑1 (MQ=255)
                                                                                                                                                                                                          aCTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTa                                                                                                                                                                                                                 >  6:128734/1‑138 (MQ=255)
                                                                                                                                                                                                            tttAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAc                                                                                                                                                                                                         <  7:653107/144‑1 (MQ=255)
                                                                                                                                                                                                            tttAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAc                                                                                                                                                                                                         >  8:653107/1‑144 (MQ=255)
                                                                                                                                                                                                             ttAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATc                                                                                                                                                                                                                        >  2:604281/1‑128 (MQ=255)
                                                                                                                                                                                                                 cGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGc                                                                                                                                                                                                                    <  4:604279/128‑1 (MQ=255)
                                                                                                                                                                                                                        tACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAAt                                                                                                                                                                                         >  5:294788/1‑148 (MQ=255)
                                                                                                                                                                                                                         aCGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAAtt                                                          <  2:351860/274‑1 (MQ=255)
                                                                                                                                                                                                                          cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCa                                                                                                                                                                                                             >  4:280988/1‑126 (MQ=255)
                                                                                                                                                                                                                            aaGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTa                                                                                                                                                                                                                          <  5:813406/111‑1 (MQ=255)
                                                                                                                                                                                                                            aaGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTa                                                                                                                                                                                                                          >  6:813406/1‑111 (MQ=255)
                                                                                                                                                                                                                            aaGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTg                                                                                                                                                                                      <  6:117281/147‑1 (MQ=255)
                                                                                                                                                                                                                              gATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAc                                                                                                                                                                                                         <  2:280989/126‑1 (MQ=255)
                                                                                                                                                                                                                               aTTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAgg                                                                                                                                                                                  <  5:452727/148‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACttt                                                                                                                                                                             <  6:294788/148‑1 (MQ=255)
                                                                                                                                                                                                                                       cGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCCATATATTGAATCTGCAGGACTTTGTAGGCCGGCTATGTCGTTAACGCCGCATCCGGCATTAACTAAGCGCATTTGTCTGCAAGCTGAGTAACGTCTAATCGTTTTTCTTAATCACTCATCAATTCACTTACATTAATTCTCAAACATCCATg                                            >  2:4667/1‑274 (MQ=38)
                                                                                                                                                                                                                                         tagtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCAttt                                                                                                                          >  4:587633/1‑194 (MQ=255)
                                                                                                                                                                                                                                         tagtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTCATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGcc                                          >  2:237148/1‑274 (MQ=255)
                                                                                                                                                                                                                                          agtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGc                                                                                                                                                                       <  5:622570/148‑1 (MQ=255)
                                                                                                                                                                                                                                             aaagggttctgttAAGTGATAATCCGTAAGCCCGATACTGTCGGGGTGACGGCTGCTAATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCCCTTCCATTAATTCTCAAACATCCAGGCCCAAcg                                     <  4:237147/263‑1 (MQ=38)
                                                                                                                                                                                                                                             aaaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCa                                                                                                                      <  2:587635/194‑1 (MQ=255)
                                                                                                                                                                                                                                              aagaggcagAGTAAGTGATTATCGTTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGACCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAAcg                                     <  2:406034/266‑1 (MQ=255)
                                                                                                                                                                                                                                                gggATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATaa                                                                                                                                                                <  7:544755/149‑1 (MQ=255)
                                                                                                                                                                                                                                                gggATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATa                                                                                                                                                                 >  6:298308/1‑148 (MQ=255)
                                                                                                                                                                                                                                                 ggATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGAc                                                                                                                                                                                                                                                         >  2:247333/1‑59 (MQ=255)
                                                                                                                                                                                                                                                   aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                                                                                                                                                             <  8:169029/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                                                                                                                                                             <  7:358633/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     cGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAAt                                                                                                                                                                                                                                                     <  4:247332/59‑1 (MQ=255)
                                                                                                                                                                                                                                                         tAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAc                                                                                                                                                       <  7:672019/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          aaGTGATTATCCGTCACCCCGATACTGTCGGGGTGGCGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCaa                                                                                                                                                                                                   <  4:354951/104‑1 (MQ=255)
                                                                                                                                                                                                                                                           aGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGc                                                                                                                                                                                                                             >  4:537206/1‑77 (MQ=255)
                                                                                                                                                                                                                                                              gATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgccg                                                                                                                                                   <  5:298308/148‑1 (MQ=255)
                                                                                                                                                                                                                                                               aTTATCCGTCACCCCGATACTGTCGGTGTGACGGCGGCTGATGACAACTGCAAAATAGCCTGATGCGCTACGCTTCTCAGGCCTACGCAATACCTGCACTATATTGAATCTGCATGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGCGTAACGTCCAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAACCACCCATGCCCACCGCTCGCTTACCGTGAATg                    >  4:45668/1‑274 (MQ=38)
                                                                                                                                                                                                                                                               aTTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                                         <  2:537208/77‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGc                                                                                                                                                                                                                             >  4:816036/1‑71 (MQ=255)
                                                                                                                                                                                                                                                                  aTCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGACCGGATTATGCGTTACCGCCGCATCATGCATTAACAAAGCGCCTTTGTCAGCAATTTGAGTAACGTCGAATCCATTTTCTAAATCCCACATCAATTCACTTCCATTAATTCCAAAACATCCATGccc                                         >  4:406032/1‑250 (MQ=37)
                                                                                                                                                                                                                                                                     cGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                                         <  2:816038/71‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatata                                                                                                                                                                                               >  4:320672/1‑95 (MQ=255)
                                                                                                                                                                                                                                                                       tCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCa                                                                                                                                              >  4:36923/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                        cACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCgg          <  4:346839/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aCCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa                                                                                                                                       >  7:602952/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         aCCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCCTATCAAGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTg                                                                                                                         >  2:162714/1‑163 (MQ=255)
                                                                                                                                                                                                                                                                         aCCCCGATACTGTCGGGGTGACGGCGACTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCAATTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGACTTTTGAGGTCgg          >  2:474489/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                         aCCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAtt                                                                                                                                        >  5:900722/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          ccccGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACt                                                                                                                                                                               <  5:252108/108‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCa                                                                                                                                          <  2:36924/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc                                                                                                                                     >  7:705902/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           cccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc                                                                                                                                     >  7:482962/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            ccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc                                                                                                                                     >  5:436350/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                              gATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTg                                                                                                                                                                            <  8:833576/107‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTa                                                                                                                                                                          >  6:732898/1‑109 (MQ=255)
                                                                                                                                                                                                                                                                                tACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTa  >  2:210882/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                  
CTGTTTGGTAAGTTAATGATGCAGAAGCCGAACATTCTGATCATGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTA  >  NC_000913/857304‑857850

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