Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,285,320 +A intergenic (‑107/+93) yjcH ← / ← acs DUF485 domain‑containing inner membrane protein YjcH/acetyl‑CoA synthetase (AMP‑forming)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,285,3201.A97.6% 292.5 / ‑0.6 84intergenic (‑107/+93)yjcH/acsDUF485 domain‑containing inner membrane protein YjcH/acetyl‑CoA synthetase (AMP‑forming)
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (43/39);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

AGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TT‑TACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGT  >  NC_000913/4285078‑4285559
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
aGGCGATCAGTAAAAGAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGACG‑AATGTAGGCCGGATAAGGCG‑TTATAcgccg                                                                                                                                                                                                                                          <  2:257988/251‑1 (MQ=255)
 ggCGATCAGTAAAATAAAGCCGCTATAAACTGCCAGCATAATAATCGACAGGATGGAGGCAAACCGTTGCAGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCACTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATAcgccgc                                                                                                                                                                                                                                         >  2:453887/1‑251 (MQ=255)
 ggCGATCAGTAAAATAAACCCGATATAAACTGCCAGCCTAATCATCGCCAGGATGGTGGCAAACCGCTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCCGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTATCCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATAcgccgc                                                                                                                                                                                                                                         >  2:444964/1‑251 (MQ=255)
   cGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCAt                                                                                                                                                                                                                                       >  1:99643/1‑251 (MQ=255)
     aTCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCagcc                                                                                                                                                                                                                                     <  2:75365/251‑4 (MQ=255)
        aGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTACTGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑CATGTAGGCCGGATAAGGCG‑TTATACGCAGCCTCCGGc                                                                                                                                                                                                                                  >  2:330929/1‑251 (MQ=255)
         gTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCTACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACGCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGca                                                                                                                                                                                                                                 <  2:284398/251‑1 (MQ=255)
         gTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGc                                                                                                                                                                                                                                  >  1:96625/1‑250 (MQ=255)
            aaaTAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAATCCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGcaatc                                                                                                                                                                                                                              <  1:336277/251‑1 (MQ=255)
                 aaGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGc                                                                                                                                                                                                                         >  2:149829/1‑251 (MQ=255)
                       atatAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATg                                                                                                                                                                                                                   >  1:284643/1‑251 (MQ=255)
                         ataAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATTGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTAGACGCCGCATCCGGCAATCAATGCCTGATGc                                                                                                                                                                                                                  <  2:443750/250‑1 (MQ=255)
                         ataAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACAGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATCAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCg                                                                                                                                                                                                                 >  1:467275/1‑251 (MQ=255)
                           aaaCTGCCAGCATAATAAGCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTATGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGAc                                                                                                                                                                                                               <  1:473483/251‑1 (MQ=255)
                              cTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCCATCAATGCCTGCTGCGACGCt                                                                                                                                                                                                            >  2:209540/1‑251 (MQ=255)
                                gCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                                                                                                                                                                                          >  2:390377/1‑251 (MQ=255)
                                     cATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGGTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCGGTA‑TACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGt                                                                                                                                                                                                     >  2:369551/1‑251 (MQ=255)
                                       taataaTCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                                                                                                                                   >  1:239084/1‑251 (MQ=255)
                                               gACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCTCATTGTCTTCTATCCGCTGATAAATCGTACCAGTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATGCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCTTCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                                                                                                                                                           <  2:325745/251‑1 (MQ=255)
                                                  aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCt                                                                                                                                                                                         >  2:440296/1‑250 (MQ=255)
                                                  aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCt                                                                                                                                                                                         >  1:131871/1‑250 (MQ=255)
                                                    gATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAc                                                                                                                                                                                       >  1:148106/1‑250 (MQ=255)
                                                          gtcAAACCGTTGCCGTTTGCCGACTAACTCCCTGAAATGCGCAGTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCCCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCGGATGCGACGCTGTCGCGTCTTATCAGGCCTTCAAACCg                                                                                                                                                                                 <  1:254342/248‑1 (MQ=255)
                                                                  gTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACTCCTCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                                                                                                                               <  2:314178/228‑1 (MQ=255)
                                                                  gTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                                                                                                                               >  1:314178/1‑228 (MQ=255)
                                                                   ttGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGcc                                                                                                                                                                                          <  1:154910/232‑1 (MQ=255)
                                                                   ttGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGcc                                                                                                                                                                                          >  2:154910/1‑232 (MQ=255)
                                                                    tGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTACCATTCATCAAAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCAGCTCACTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCAGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGCCTCg                                                                                                                                                                      >  1:457631/1‑251 (MQ=39)
                                                                          ttttCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCGCAGATTCTCCAGAGGTACTGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGcctccg                                                                                                                                                                >  2:435750/1‑246 (MQ=39)
                                                                              cGACTAACTCCCTGAAATGTGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACCGCTTCTCCAGAGGTAATGTCGGGCTTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGAAGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACAGACTCGCATCCGGCa                                                                                                                                                             >  2:434770/1‑250 (MQ=38)
                                                                              cGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCa                                                                                                                                                             <  1:253701/250‑1 (MQ=255)
                                                                                    aCTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAAGCGTACCATTCATCACAGATTCTCCAGAGGGAATGTAGGGATTGGTTTAATTCCCGCTCCCTTAAGGGAGAAGGTTAACGCTCGGGTAACCCTGGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGAGGCGACGCTGTCGCGTCTTATCAGTCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTgg                                                                                                                                                      <  2:300353/251‑1 (MQ=255)
                                                                                       cccTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTGAATGCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAACCGTTACCGACTCGCATCCGGGAATTGTggg                                                                                                                                                     <  2:465970/249‑1 (MQ=255)
                                                                                       cccTGAAATGCGCATTGTCGTCTATCCGCTGATAAATCGTACCATTCAACACAGATTCTCCAGAGGTAATTTAGGGATTGTTTTAATTCCCGCCCCCTGATGGGAGAAGGGTAACGCTCGGGTATCCCGGGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCAGCCGGCAATCAAGGCCTGATTCGACGCTGTCGTGTCTTATCAGGCCTACAAACCGTTACCGACTCGCACCCGGCAATTGGGGGtt                                                                                                                                                   <  2:362542/251‑1 (MQ=255)
                                                                                                    aTTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCCATCAATGCCTGCTGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATcgc                                                                                                                                      >  2:67666/1‑251 (MQ=255)
                                                                                                       gTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGAt                                                                                                                                   >  1:360992/1‑251 (MQ=255)
                                                                                                            ctATCCACTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCt                                                                                                                              >  1:84912/1‑251 (MQ=255)
                                                                                                                   cTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATACCCTGCTtctc                                                                                                                       >  2:256339/1‑251 (MQ=255)
                                                                                                                    tGATAAATCGTACCATTCATCACAGATTCTCCCGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTctct                                                                                                                      >  1:281521/1‑251 (MQ=255)
                                                                                                                    tGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTctct                                                                                                                      >  2:30172/1‑251 (MQ=255)
                                                                                                                         aatcgtaccATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCa                                                                                                                   >  1:25264/8‑249 (MQ=255)
                                                                                                                          atcgtaccATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAagc                                                                                                                <  1:154685/245‑1 (MQ=255)
                                                                                                                             gtaccATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGATAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATCCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGCCTCGCATCCGGCAATTGTGGGTTACGCTGGCATCGCGATAGCCTGCTTCTCTTCAagcagc                                                                                                             >  2:179078/4‑251 (MQ=39)
                                                                                                                             gtaccATTAATCACAGATTCTCCAATGGTACTGTAGGGATTGATTTAATTACCGCTCCCTTAAGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑TATGTAGGCCGGATAAGGCGCTA‑TACGCCGCACCCGGCCATCCATGCCTGATGCGCCGCTGTCGCGTCTTATCAGGCCTACAAACCGTTCCCGACTCGCATCCGGAAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTATTCAaccagc                                                                                                             >  2:462483/4‑251 (MQ=255)
                                                                                                                               accATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGAGATAGCCTGCTTCTCTTCAAGCAGCtt                                                                                                           >  1:38399/2‑251 (MQ=255)
                                                                                                                                    tcatcaCAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCTATCAATGCCTGATGCGCCGCTGTCGCGTCTTATCATGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCg                                                                                                       <  1:323981/250‑1 (MQ=255)
                                                                                                                                        cacaGATTCGCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCGCCCTTATGGGAGAAGGTTAACGCTCGGGTATCCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATTCCTGATGCGACGCTGTCGCGTCTGATCAGGCCGTCAAACCGTTACCGACGCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCGCTTCAAGCAGCTTCTCGACTAc                                                                                                  <  1:299487/251‑1 (MQ=255)
                                                                                                                                         acaGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACg                                                                                                 >  2:401183/1‑251 (MQ=255)
                                                                                                                                               tctcCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTAGGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTGCTCTTCAAGCAGCTTCTCGACTACGCCAgg                                                                                            <  1:434770/250‑1 (MQ=255)
                                                                                                                                               tctcCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTACCTTATGGGAGAAGGTTAACGCTCGGGCAACCCTTGCCG‑AATGAAGGCCGGATAAGGCG‑TTATACGCCGCATCCGACAATCAATGCCTGATGCGACGCTGTCGCGTATTATCAGGCCTACAAACAGTTACCGACTCGCATACGGCAATTGTGGGTTACAAGGGCATCGCGACAGCCTGCTTCTCTTCAAGCAGCTTCTTGACTACGCCAggc                                                                                           >  1:311521/1‑250 (MQ=255)
                                                                                                                                                  ccAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGTAATCAAGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAATCGTTGCCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCGTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCg                                                                                        <  1:290178/251‑1 (MQ=255)
                                                                                                                                                        gTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGCTAAGGCG‑TTATACGCCGCATACGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCCAGc                                                                                  >  1:303914/1‑251 (MQ=255)
                                                                                                                                                           aTGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTc                                                                               <  1:121940/251‑1 (MQ=255)
                                                                                                                                                            tGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCg                                                                              >  2:274116/1‑251 (MQ=255)
                                                                                                                                                             gTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACGACGCCAGGATCGGCAAGCGTCg                                                                              <  1:141309/250‑1 (MQ=255)
                                                                                                                                                               aGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAgg                                                                           <  1:453887/251‑1 (MQ=255)
                                                                                                                                                               aGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAg                                                                            <  2:284763/250‑1 (MQ=255)
                                                                                                                                                                gggATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTACCCCTTGCCG‑AATGGAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGt                                                                          <  2:128418/251‑1 (MQ=255)
                                                                                                                                                                         ttAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCa                                                                  <  1:181961/250‑1 (MQ=255)
                                                                                                                                                                         ttAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCa                                                                  <  1:400248/250‑1 (MQ=255)
                                                                                                                                                                         ttAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCa                                                                  <  1:38988/250‑1 (MQ=255)
                                                                                                                                                                          tAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAgg                                                                >  2:5237/1‑251 (MQ=255)
                                                                                                                                                                            aTTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGtt                                                              <  2:99643/251‑1 (MQ=255)
                                                                                                                                                                            aTTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGtt                                                              <  1:30172/251‑1 (MQ=255)
                                                                                                                                                                                ccGCTCCCTTAGGGGAGAAGGTTAACGCGCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTg                                                          <  1:248936/251‑1 (MQ=255)
                                                                                                                                                                                 cGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTgg                                                         <  2:107429/251‑1 (MQ=255)
                                                                                                                                                                                      cctgaTGGGAGAAGGTTAACGCTCGGGTAACCCTGGCCG‑AATGGAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATcg                                                    <  2:284643/247‑1 (MQ=255)
                                                                                                                                                                                        ttATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGCGGTATCGCCCAGGTTGCTGGTATcgcc                                                  >  1:192462/1‑251 (MQ=255)
                                                                                                                                                                                          aTGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATcgccc                                                 <  1:250388/250‑1 (MQ=255)
                                                                                                                                                                                                  aGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGGTGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGCATCGCCCAGGTTGCTGGTATCGCCCGCCGCAAtt                                        >  2:21916/1‑251 (MQ=255)
                                                                                                                                                                                                     ttAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTg                                     >  2:454785/1‑251 (MQ=255)
                                                                                                                                                                                                             cGGGTACCCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGGGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAAt                             <  1:440296/251‑1 (MQ=255)
                                                                                                                                                                                                             cGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGaa                              <  1:67666/250‑1 (MQ=255)
                                                                                                                                                                                                             cGGGTAACCCTTGCCG‑AAGGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAGCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGa                               <  1:390377/249‑1 (MQ=255)
                                                                                                                                                                                                              ggggtaCCCTTGGCGGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAAACAATGCCTGATGCGCCGCTGGCGCGTCTTATCAGGCCTACAAACCGTTCCCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGGATCGCCCGCCGCAATTTTGCGCAGAAt                             <  2:308020/246‑1 (MQ=255)
                                                                                                                                                                                                                     ccTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATTCGCCGCCTCCGGCAATCACTGCCTGATGCGCCGCTGACGCGTCTTATCAGGACTACAAACCGTTACCGCCTCGCATCCGGCAATTGCGGGTTACGATGGAATCGCGAGAGCCTGCTTctct                                                                                                                      >  2:266483/1‑154 (MQ=255)
                                                                                                                                                                                                                     ccTTGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTGACCGACTCGCATCCGGCAATGGTGGGTTACGATGGCATCGCGATAGCCTGCTTctct                                                                                                                      <  1:266483/154‑1 (MQ=37)
                                                                                                                                                                                                                        tGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATc                                                                                                                                        >  2:149965/1‑133 (MQ=255)
                                                                                                                                                                                                                        tGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATc                                                                                                                                        <  1:149965/133‑1 (MQ=255)
                                                                                                                                                                                                                        tGCCG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATa                   >  2:213183/1‑250 (MQ=255)
                                                                                                                                                                                                                          ccG‑AATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAAtt                >  1:267261/1‑251 (MQ=255)
                                                                                                                                                                                                                           tg‑AATGTAAGACGGATAAGGAG‑TTATCCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAAACGTTACCGACTCGCATCCGGCTATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGACTACGGCGCATAAttt               >  2:52038/2‑251 (MQ=38)
                                                                                                                                                                                                                                       cGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCgg    >  1:185871/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGt  >  1:254091/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
AGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCG‑AATGTAGGCCGGATAAGGCG‑TT‑TACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGT  >  NC_000913/4285078‑4285559

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