Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,489,889 A→C 15.9% V59G (GTG→GGG)  yjgR ← DUF853 family protein with NTPase fold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,489,8890AC15.9% 163.7 / 3.6 79V59G (GTG→GGG) yjgRDUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base A (30/34);  new base C (2/12);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 3.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

AGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAA  >  NC_000913/4489664‑4490112
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                               
aGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGACTGCCCGACGCCGGTCAGGTCGCCTTTcacc                                                                                                                                                                                                                               >  1:505745/1‑227 (MQ=255)
  cACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGCTATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCCCGACGCCGGTCAGCTCGCCTTTCACAt                                                                                                                                                                                                                              >  2:643902/1‑227 (MQ=255)
  cACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTCCCAGTCATTTACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCCAACCCGTTCCTGCCTGCGCGACCCCGGTCAGATAGCCTTTcacaa                                                                                                                                                                                                                              >  1:40789/1‑226 (MQ=255)
  cACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATACCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATCTTTTTAAGCCTTGCGAGCAGTTTTTCCGAAACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTaacct                                                                                                                                                                                                                              >  2:657739/1‑225 (MQ=255)
   cccccatgttGCACACCTTTGGGCTTCAACAGGCGTGCCAGCAACAGCGGCCCCAGGTCCGAAACCGGCGCCCGCACCGGATGGCCTTTCCCGCCAAAGATATCCCCCACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACAt                                                                                                                                                                                                                              <  1:93212/218‑1 (MQ=255)
   cacccAGATTGCCCTTCGTTGAGATTCCACAGTCGAGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTCCTCGCCAAAGATACCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACAt                                                                                                                                                                                                                              <  1:106170/223‑1 (MQ=255)
   aCGCCAGATTGCCCATCGTTGAGGTTCCACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCTCCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACAt                                                                                                                                                                                                                              <  2:186117/226‑1 (MQ=255)
         gATTGCACATCGTTGAGATTCAACAGGCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATa                                                                                                                                                                                                                      <  2:617861/228‑1 (MQ=255)
                   cGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCAAACCACCACCGGATTGGCATGCGGTTGCCCGGCCTTGACGCCGATATTTTTAAGCCTTGCGCGCAGTTTTTCCGACACCGTTCCTTCCTGCGCGACCCCGGTCAGATCGCCTTTCCGATCAGCCATAAAAAccccca                                                                                                                                                                                                            >  1:649398/1‑224 (MQ=255)
                      tgggaTTCAACGGGCGTGCCCGCCACAGCGGCCCCAGGTCCGAACCCGTCGCCCGCTCCGGTTGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGc                                                                                                                                                                                                         <  1:549950/225‑1 (MQ=255)
                      tGAGATTCAACAGTCGGGCCAGCAACAGCGGCCCCAGACCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGc                                                                                                                                                                                                         <  2:545178/228‑1 (MQ=255)
                        agaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCCTAAACACCGGAACGccc                                                                                                                                                                                                       >  2:315337/1‑227 (MQ=255)
                        agaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCg                                                                                                                                                                                                       >  2:470923/1‑228 (MQ=255)
                             cAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTTCCCGACGCCGGTCAGATCGCCTTTCCCATCAGCCATAAACACCGGCACGcccaatt                                                                                                                                                                                                   >  1:316857/1‑222 (MQ=255)
                             cAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGCCCCCGGTCAGCTCGCCTTTCACACAAGCCATAACCACCGGCACGCCGCttt                                                                                                                                                                                                   >  2:33910/1‑227 (MQ=255)
                             cAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAAACCCGGCACCCCGCttt                                                                                                                                                                                                   >  2:630188/1‑227 (MQ=255)
                              aaccggcGTGCCAGCAACAGCGGCCCCAGGTCCGAAACCTTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCCCCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACGCCGGCACGCCGATTTcc                                                                                                                                                                                                 <  2:40789/222‑1 (MQ=255)
                                  gcggcgttcgcgaacgCGGCCCCAGACCCGAAACCGTCGCGCGCGCCGGATGGCCTTTCTCGACAAAGATATCCCACACCACCAACGGAGTGGCGGGCGGTTGCCAGTCATTGAAGCCGATATGTTTAAGCCGTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCTGGCCGCCGGTCAGATCGCCTTTCCCATCAGCCATAAACACCGGCACGCCGATATCCGACa                                                                                                                                                                                             <  1:612500/213‑1 (MQ=255)
                                     gTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATg                                                                                                                                                                                          >  2:72672/1‑228 (MQ=255)
                                        ccAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGAAACCGTTCCTGCCTGCGCGCCGCCGGTCCGATCGCCTTTCACATCAGCCATAAACACCGGCAGCCCGATCTCCGCCCATGACt                                                                                                                                                                                       >  2:4954/1‑228 (MQ=255)
                                           ctacgactgggcccAGGCCCTAACCCCTCGCCCCGACCGGATGGCCTTTTCCCGCAAAGGTGACGCCAGCCACCACCGTATTGGCATTCGGTTGCCATACATTGACGCCGTTATTTTTAAGCTTTGCGAGCAGTTTTACCGACACCGTCCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAACCACCGGTACGCCGATTTCCGACACTGCCTCCg                                                                                                                                                                                    <  1:336465/217‑1 (MQ=255)
                                                 gCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGtt                                                                                                                                                                              >  2:341112/1‑228 (MQ=255)
                                                        cAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTg                                                                                                                                                                          <  2:446019/225‑1 (MQ=255)
                                                          gATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACCCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACAGGCACGCCGATTGCCGACCATGCCTCCGCCAGATTTTGCAGc                                                                                                                                                                      >  2:160640/1‑227 (MQ=255)
                                                            tCCGAAACCGTCCCCCCGCCGGTGTGGCCTTTCCCCCCAAAAAACTCCCACCCCACCACCGGATGGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTGCCGACCCCGTTCCTGCCTGCGGGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTa                                                                                                                                                                   <  1:72672/228‑1 (MQ=255)
                                                              ccaaaCCGCCGCCCGCACCGGTGTGCCTTGCGCCCCAAAGATATCCCACCGCCCCACCGGATTGGCGTGCGGTTGCCAGTCATTGACGGCGATATTTTTAAGCCTTTCGAGCAGTTTTTCCGACCCCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTGTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAAc                                                                                                                                                                 <  1:185926/226‑1 (MQ=255)
                                                              cGAAACCGTCGCCCGCACCGGATGCCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTa                                                                                                                                                                   <  2:46373/226‑1 (MQ=255)
                                                              cGAAACCGGCGCCCGCCCCGGATGGCCGTTCTCGCCAAAGATGTTCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTa                                                                                                                                                                   <  2:289480/226‑1 (MQ=255)
                                                                    cctCGCCCGCCACCGGGTGGCCTTTGCCGCCAAGGATACCCCACACAGCCCCCGGGTTGCCAGGGGGTTGGCAGTGATTGAGCCCGGTAGTTTTAAGCCTTCCAACCAGTTTGTCCAACACCGGCCCTGCCTGGGGGACGCTGGTCAGATTGCCTTTCTCATCAGCCATAAACCCCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGttt                                                                                                                                                            <  1:188926/226‑1 (MQ=255)
                                                                      tcgctcgcACGGGAAGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGATGCCAGTCATTGACGCCGATAACTCTAAGCATTGCGAGCAGCTTCTCCGACACCGTTCCCCCCGCCGCGACCCCTGTCACATCCCCTTTCACTTCAGCCATAAACCCCGGCCCGACGACTTCCCACAATCACCCCCCCAGGTTTTCCACCGTAACGGCTTTa                                                                                                                                                          >  1:42080/1‑227 (MQ=255)
                                                                        gCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCCCACCACCACCGGCTTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTGCCTGCCTGCGCGCCGCCGGTCAGAGCGCATTTCACCTCAGCCATCAACACCGGCCCCCCGCTTTCCGACAATACACCCCCCAGTTTTTGCCGCGTAACGGTTTTAccc                                                                                                                                                       >  1:61032/1‑228 (MQ=255)
                                                                                          cTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTAATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGCCAGATCGCCTTTCACATCCGCCATAAACACCGGCACCCCGACTTCCGACAATGACCCCGCCAGTTTTACCACGGTACCGGTTTTACCCCTCCCCGTTggtcccg                                                                                                                                       >  1:137520/1‑222 (MQ=255)
                                                                                              ttccccAAAGATCTCCCCACCCACCTCCGTGTTGGCTTGCTGTTGCCAGTCATTGACGCCGATTTTTTTAAGCCTTGCGAGCAGTTTTCCCGACCCCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAACCACCGGCACGCCGATTTCCGCCAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAAt                                                                                                                                   <  1:339373/224‑1 (MQ=255)
                                                                                              cgcccaaaaagtatCCCACCCCCCCCACGGGTGTGCAGGCGGTTGCCCGGCAGTGACCGCGAATGTTTTAAGCCTTGCGAGCAGTTTTTCCGCCTCCGTCCCTGCCGGCGCGCCGCCGGTCAGATCTCCGTTCCACTCAGCCATAAGCACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAAt                                                                                                                                   <  1:80217/217‑1 (MQ=255)
                                                                                                 gCCAAAGATATCCCACCCCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCa                                                                                                                                 <  2:496854/227‑1 (MQ=255)
                                                                                                      atatatCCCACACCGCCGCCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCATTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCAt                                                                                                                           <  2:649398/226‑1 (MQ=255)
                                                                                                      aggtatCCCACGCCCCCACCGACTGTGCAGGCGGTTGCGAGTAGTCGCCTCCGACGTTTTTAAGCCTTGCGAGCATGTTTTCCGCCGCCGGTCCGGCCTGCGCGCCGCCGGTCGAATAGCCTTTCCCAACAGCCATAAACACCGGCAGGCCGATTTCCGACAATGACTCCACCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCAt                                                                                                                           <  1:45755/225‑1 (MQ=255)
                                                                                                            ccccCACCCCCCCCGGATTGCCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGt                                                                                                                     <  2:505745/228‑1 (MQ=255)
                                                                                                                caccaccacCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTcaca                                                                                                                                                                                                                               >  2:88141/1‑118 (MQ=255)
                                                                                                                 accaccaccGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACAt                                                                                                                                                                                                                              <  1:88141/118‑1 (MQ=255)
                                                                                                                 accaccaccGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAACCACCGGCACGCCGCTTTCCCACAATGACTCCCCCAGTTTCTCCAGCGTAACGGCTTTACCCGGCCCCGTTGCGCAGGGCATCCCCCCATCAccgacaccc                                                                                                                >  1:115759/1‑220 (MQ=255)
                                                                                                                 accaccaccGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAACCACCGGCACGCCGATTTCCGACAATGCCTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCCTCCCCGTTGCGCCGGTAATCAGCCCATTACCGTTAGcc                                                                                                                >  1:108586/1‑228 (MQ=255)
                                                                                                                 accaccaccGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATCACGGTTAGcc                                                                                                                >  2:62117/1‑228 (MQ=255)
                                                                                                                 accaccaccGGATTGGCATGCGGTTGCCAGTCATTGAAGCCGATATTTTAAGCCTTGCGAGCAGTTTTTCCGCCACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGCCAATGACTCCCCCAGTTTTTGCAGCGTAACGGTTTTACCGGCCCCCGTGCGCCGGAAATCACCCCATGAAGGGTAGCCAt                                                                                                              >  1:460493/1‑228 (MQ=255)
                                                                                                                     ccaccGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGCTTTCCGACAATGACTCCGCCAGTTTTTGCAGGGTAAGTGGTTTACCCATCCCCGTTGCGCCGGAAACCACCCAATTCAGGTTACCCTTTc                                                                                                            >  1:105355/1‑228 (MQ=255)
                                                                                                                     ccaccGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACCCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGCCAATGACTCCGCCCGTTTTTGCAGCGTAACGGTTTTACCCCCCCCCGTTGCCCCGGCACTCAGccccttcagggtagcccttt                                                                                                            >  1:410180/1‑211 (MQ=255)
                                                                                                                      caccGTATTGGCATGCGCGTCCCAGTCATTAACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCCGATCGCCTTTCACATCCGCCATAAACACCCTCACCCCGATTTCCGACAATGACTCCGACACTTTTGCCAGCGTAACGGTTTTACCCGTCACCTTTTCCCCGCTAACCACACCTTGACGCTTAGACCTTcc                                                                                                           >  1:407455/1‑228 (MQ=255)
                                                                                                                        ccGGATTGGCTGGCGGTTGCCCGTCATTGACCGCCATATTTTTAAGCCTTGCGATCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCCCATCAGCCATAAACACCTGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGGAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCg                                                                                                         <  2:581279/228‑1 (MQ=255)
                                                                                                                             ttggcagtgggTTACCGGTCCTTGACGCCGATTTTTTTAAGCCTTGCGAGCAGTTTATCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGt                                                                                                    <  2:325505/219‑1 (MQ=255)
                                                                                                                             ttGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAgc                                                                                                    >  2:457759/1‑227 (MQ=255)
                                                                                                                             ttGGCATGCGGTGGCCGGGCATTGACGCCGATATTTTTAAGCCTTTGCAGCAGTTTTTCCGACACCGTTCCTGCCTGCGAGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGt                                                                                                    <  2:426148/228‑1 (MQ=255)
                                                                                                                               ggCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTaa                                                                                                  >  2:492891/1‑228 (MQ=255)
                                                                                                                                  ttccgtttGCCCGTCATTGCTGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTCCCTGCCTGCGCGCCGCCGGTCAGATCCCCTTTCCCATCACCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTaaaac                                                                                               <  1:415763/222‑2 (MQ=255)
                                                                                                                                     cGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACa                                                                                             >  2:415763/1‑227 (MQ=255)
                                                                                                                                        tttccAGTCATTGCCGCCTATATTTTTAAGCCCTGTGAGCAGTTTTTCCGCCACCGTTCCTGCCTGCGCGACCCCGGTCAGATCGCCTTTCACATCCGCCATAAACCCCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTc                                                                                         <  2:315683/225‑1 (MQ=255)
                                                                                                                                         tgcacgtattgcaccccGATTTTTTCAAGCCTGGGGAGCTGTTTCTCCGACCCCGTCCCTGCCTGCGCGACGCCGGTCAGGCCGCCTTTCCCATCAGCCATAAACACCGGCACGCCCATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTc                                                                                         <  1:657739/212‑1 (MQ=255)
                                                                                                                                          gCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCtg                                                                                       >  2:592886/1‑228 (MQ=255)
                                                                                                                                              tcatctgcccacgttttttttAGACTTTGGGGGCAGGTTTTCCTGTCCCTTTCCTTCCTCGGGGCCGCCGGTCGCATCGCCTTTCCCATCAGCCATAAACACCGGCCCGCCAATTGCCGACAATGACGCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTc                                                                                   <  1:437537/212‑1 (MQ=255)
                                                                                                                                              gTCATTGCCCCCGATTTTTTTAAGCCTGGCGAGAGGTTTTCCCGCCCCCGTCCCTGCCTGCGGGACGCCGGTCAGACCGCCTTTCCCATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTc                                                                                   <  1:470923/228‑1 (MQ=255)
                                                                                                                                                       ccggATTTTTTTAACCCTGGAGAGCATTTGTGCCCACACGGTTCCTGCCTGCGCGACGCCGGTCAGGTCGCCTTTCACATCAGCCATAAACCCCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTACCGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCAGGACGGTTAGCCATTCCCGGGAGTAAAAACAGTTCTGTGTCCGGCGTGc                                                                           <  1:371492/224‑1 (MQ=255)
                                                                                                                                                         cgtttttttaAAGCCTTGCGGCGCGTTTTCCCGACGCCGTTCCTCGCTGGGAGCCGCCGGTCAGTGCGCCTTTCCCATCAGCCAAAAACCCCGGCACGCCGATTGCCGACAATGGCTCCGCCAGTTTGTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCggg                                                                        <  1:4954/223‑1 (MQ=255)
                                                                                                                                                         cGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCACCCCATGTCGGTTAGCCCTTCCCGGCCGTAAAAACAATTCTGTGTCCGGCGTGCggg                                                                        >  2:257423/1‑228 (MQ=255)
                                                                                                                                                           atatTTTTAAGCCTTGCGATCAGTTTTTCCGAACCCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACGCCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCgg                                                                         <  2:367679/225‑1 (MQ=255)
                                                                                                                                                                  gaacccTTGCGCGCAGTCTTCCCGACACCGTTCCCGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGTCAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACa                                                               <  2:134421/224‑1 (MQ=255)
                                                                                                                                                                   cagccttgccgtgcgTTTTGCCGACCCCCTTCAGGCCTCGGGGGTGCCGGACCGATCGGTTTTCACATCGGCCATAAACACCGGCACGCTGATGTCCGGAACTGACTCCGCCAGGTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAg                                                              <  1:160640/214‑1 (MQ=255)
                                                                                                                                                                       cTTGCAAGCCTTTTTTCCGACACCGTTCCTGCCGGCGCCACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGt                                                          <  1:315337/228‑1 (MQ=255)
                                                                                                                                                                        ttGCGAGCAGTTTTTCCGGCCCCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACCCCGGCACGCCGATTCCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGtt                                                         <  2:215574/228‑1 (MQ=255)
                                                                                                                                                                        tgtcGAGCAGTTTTTCCGACACCGTACCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGtt                                                         <  2:2412/225‑1 (MQ=255)
                                                                                                                                                                        gtccgggCAGTTTTTCCGAAACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGtt                                                         <  2:234138/222‑1 (MQ=255)
                                                                                                                                                                        ctGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAACCCCCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGtt                                                         <  2:620317/227‑1 (MQ=255)
                                                                                                                                                                         tGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTTTGTCCGGCGTTCGGGCAATTAACAGGGGtta                                                        >  2:140073/1‑227 (MQ=255)
                                                                                                                                                                          gtgAGCAGTTTTCCCGACCCTGTGCCCGCCTGCGGGCCCCCGGTCAGTTCCCCTTCCACATCGGCCATAAACACCGTCACGCCGATTTCCGACAATGACTCCGCCAGTTTGTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCa                                                       <  1:33910/226‑1 (MQ=255)
                                                                                                                                                                          gggagcgttttttCCGACCCCGTTCCTGCCTGCGCGCCGCCGGTCAGATCGCCTTTCACATCAGCCAGAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCa                                                       <  1:643902/220‑1 (MQ=255)
                                                                                                                                                                                  tttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGCCAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTTAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATgag                                               >  2:280565/1‑228 (MQ=255)
                                                                                                                                                                                                 tcccgcccggccgACGCCGGTCAGATCGCCTTACCAATCAGCCATAAACGCCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTc                                <  2:221188/217‑1 (MQ=255)
                                                                                                                                                                                                 tCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGt                                                                                                    <  2:360793/160‑1 (MQ=255)
                                                                                                                                                                                                       cgtcggaccGCCGGTCGAAGCGCTATTCCCATCAGCCATAAACACCGGCACTCCGAGTGCCGACAATGGCTCCGCAGGTTTTTGCAGCGTAACGGGTTTACCCGTCCCCGGTGCGCCGGTAATCAGCCCAGGGCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCTCTCATGAGATATCCTCCAGTTTcctgcct                          <  1:211296/220‑1 (MQ=255)
                                                                                                                                                                                                        tgcGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATACTCCAGTTTcctgcctg                         >  2:101025/1‑228 (MQ=255)
                                                                                                                                                                                                          cgcgACGCCGGTCAGATCGCCTTTCACATCAGCCATACACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAACCAGTTCTGTGTCCGGCGTGTGGGCATTTAACAGGGGTCGACTCATAGGATACAACCCACGTTCCTGCCAGGa                       >  1:25313/1‑228 (MQ=255)
                                                                                                                                                                                                          cgcgACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGGGCGGGCAATTAACAGGGGTTCACTCATTAGATATCCTCCAGTTTCCTGCCTGGa                       >  2:68209/1‑228 (MQ=255)
                                                                                                                                                                                                                   ggTCAGATCGCTTTTCCCATCAGCCATAAACACCGGCACGCCGAGTGCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGAAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCa              <  2:105355/228‑1 (MQ=255)
                                                                                                                                                                                                                         aTCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCTTGCCTGGAGGATAGGCACAGTcc        >  2:251725/1‑228 (MQ=255)
                                                                                                                                                                                                                                ttCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTGTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCaaa  <  2:137520/227‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                               
AGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAA  >  NC_000913/4489664‑4490112

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