Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 2,056,174 C→T T3T (ACC→ACT)  unknown → UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase (EC 6.3.2.‑)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*12,056,1740CT100.0% 154.4 / NA 54T3T (ACC→ACT) unknownUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase (EC 6.3.2.‑)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/24);  total (30/24)

CCGCTCAATTCTCCTTTAGCATGTGGGTATTCCGCTTGTTTTTCTACGATAAATTGACCTAATGTTTTCATCTTTTATTCCTTAACTTCATTAATGAATTTACTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACCAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAACGCGGCAATCCCTGCATCGAGTATATGTTA  >  1/2055938‑2056420
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
ccGCTCAATTCTCCTTTAGCATGTGGGTATTCCGCTTGTTTTTCTACGATAAATTGACCTAATGTTTTCATCTTTTATTCCTTAACTTCATTAATGAATTTACTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCatat                                                                                                                                                                                                                                          <  1:145558/251‑1 (MQ=255)
            ccTTTAGCATGTGGGTATTCCGCTTGTTTTTCTACGATAAATTGACCTAATGTTTTCATCTTTTATTCCTTAACTTCATTAATGAATTTACTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTg                                                                                                                                                                                                                              >  1:173832/1‑251 (MQ=255)
                      gtgGGTATTCCGCTTGTTTTTCTACGATAAATTGACCTAATGTTTTCATCTTTTATTCCTTAACTTCATTAATGAATTTACTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCat                                                                                                                                                                                                                                            <  1:290383/227‑1 (MQ=255)
                       tgGGTATTCCGCTTGTTTTTCTACGATAAATTGACCTAATGTTTTCATCTTTTATTCCTTAACTTCATTAATGAATTTACTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTa                                                                                                                                                                                                                   >  1:400089/1‑251 (MQ=255)
                          gTATTCCGCTTGTTTTTCTACGATAAATTGACCTAATGTTTTCATCTTTTATTCCTTAACTTCATTAATGAATTTACTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATgg                                                                                                                                                                                                                >  1:177360/1‑251 (MQ=255)
                                             aCGATAAATTGACCTAATGTTTTCATCTTTTATTCCTTAACTTCATTAATGAATTTACTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCg                                                                                                                                                                                              <  1:364458/250‑1 (MQ=255)
                                                                ttttCATCTTTTATTCCTTAACTTCATTAATGAATTTACTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCa                                                                                                                                                                           <  1:81608/250‑1 (MQ=255)
                                                                    cATCTTTTATTCCTTAACTTCATTAATGAATTTACTCACTTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTaa                                                                                                                                                                      >  1:351059/1‑251 (MQ=255)
                                                                     aTCTTTTATTCCTTAACTTCATTAATGAATTTACTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGTTATCATGTAAc                                                                                                                                                                     <  1:454720/251‑1 (MQ=255)
                                                                       cTTTTATTCCTTAACTTCATTAATGAATTTACTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAg                                                                                                                                                                   <  1:260524/251‑1 (MQ=255)
                                                                                                aaTTTACTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCata                                                                                                                                                                                                                                           >  1:773579/1‑154 (MQ=255)
                                                                                                aaTTTACTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACccac                                                                                                                                           >  1:503004/1‑250 (MQ=255)
                                                                                                     aCTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGa                                                                                                                                     >  1:590520/1‑251 (MQ=255)
                                                                                                               gAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGt                                                                                                                                                   <  1:668981/227‑1 (MQ=255)
                                                                                                                  cTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTg                                                                                                                                                                                                       <  1:39187/172‑1 (MQ=255)
                                                                                                                          attatAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAAt                                                                                                                 >  1:633075/1‑250 (MQ=255)
                                                                                                                                   gCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGaaa                                                                                                       >  1:42093/1‑251 (MQ=255)
                                                                                                                                                   aTTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCatat                                                                                                                                                                                                                                          <  1:739590/104‑1 (MQ=255)
                                                                                                                                                   aTTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGa                                                                                       >  1:800859/1‑251 (MQ=255)
                                                                                                                                                   aTTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGa                                                                                       >  1:571383/1‑251 (MQ=255)
                                                                                                                                                   aTTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGa                                                                                       >  1:672821/1‑251 (MQ=255)
                                                                                                                                                           gaTTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATAc                                                                                                 <  1:550988/233‑1 (MQ=255)
                                                                                                                                                              tCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTc                                                                                                                            <  1:656954/203‑1 (MQ=255)
                                                                                                                                                              tCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGt                                                                            >  1:721545/1‑251 (MQ=255)
                                                                                                                                                               ccTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGtt                                                                           >  1:398658/1‑251 (MQ=255)
                                                                                                                                                                cTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGagcc                                                                                    >  1:7966/1‑241 (MQ=255)
                                                                                                                                                                    ttCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACaa                                                                       >  1:477421/1‑250 (MQ=255)
                                                                                                                                                                        gtgATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTAc                                                                         >  1:394987/1‑244 (MQ=255)
                                                                                                                                                                           aTTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAgcgc                                                                >  1:496154/1‑250 (MQ=255)
                                                                                                                                                                                     caTTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATggg                                                                                                                                                                                 <  1:194971/127‑1 (MQ=255)
                                                                                                                                                                                     caTTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGt                                                      >  1:407377/1‑250 (MQ=255)
                                                                                                                                                                                     caTTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTc                                                     <  1:24066/251‑1 (MQ=255)
                                                                                                                                                                                      aTTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCa                                                    <  1:179671/251‑1 (MQ=255)
                                                                                                                                                                                                      gAAATAAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCAt                                    >  1:314392/1‑251 (MQ=255)
                                                                                                                                                                                                       aaaTGAACCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGc                                       <  1:638445/247‑1 (MQ=255)
                                                                                                                                                                                                            aaCCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAAcg                              <  1:9075/251‑1 (MQ=255)
                                                                                                                                                                                                            aaCCGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAAcg                              >  1:273958/1‑251 (MQ=255)
                                                                                                                                                                                                              ccGGTATATAAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACccac                                                                                                                                           <  1:242071/140‑1 (MQ=255)
                                                                                                                                                                                                                   atataAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGggc                                                                                                                                                                                                              >  1:712409/1‑68 (MQ=255)
                                                                                                                                                                                                                   atataAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCGAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAACGCGGCAAt                       >  1:880064/1‑251 (MQ=255)
                                                                                                                                                                                                                   atataAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGacac                                                                                                                                                         >  1:823452/1‑121 (MQ=255)
                                                                                                                                                                                                                   atataAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATAc                                                                                                 <  1:597866/177‑1 (MQ=255)
                                                                                                                                                                                                                   atataAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAACGCGGCAAt                       >  1:406643/1‑251 (MQ=255)
                                                                                                                                                                                                                   atataAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAACGCGGCAAt                       >  1:679226/1‑251 (MQ=255)
                                                                                                                                                                                                                   atataAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAACGCGGCAAt                       >  1:805079/1‑251 (MQ=255)
                                                                                                                                                                                                                     ataAAAGGACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGagcc                                                                                    >  1:365712/1‑188 (MQ=255)
                                                                                                                                                                                                                           ggACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTAccc                                                                                                                                             <  1:855238/125‑1 (MQ=255)
                                                                                                                                                                                                                            gACTCAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAACGCGGCAATCCCTGCAt               <  1:355016/250‑1 (MQ=255)
                                                                                                                                                                                                                                cAACATGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCga                                                                                                                                                                                                    <  1:483887/65‑1 (MQ=255)
                                                                                                                                                                                                                                     tGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAACGCGGCAATCCCTGCATCGAGTATATg     <  1:630032/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAACGCGGCAATCCCTGCATCGAGTATATg     <  1:283359/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tGACAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAACGCGGCAATCCCTGCATCGAGTATATg     <  1:224294/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAACGCGGCAATCCCTGCATCGAGTATATGTTa  >  1:510524/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cAACTAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAACGCGGCAATCCCTGCATCGAGTATATGTTa  >  1:159421/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
CCGCTCAATTCTCCTTTAGCATGTGGGTATTCCGCTTGTTTTTCTACGATAAATTGACCTAATGTTTTCATCTTTTATTCCTTAACTTCATTAATGAATTTACTCAATTCTGAACTTATTTTATTATAACAGCTTCAATTTTGAATGATTAAAGAGATTCCTATTTCTGTGATTCACAACACATTTTTATTGGAAACAGAAATAAACCGGTATATAAAAGGACTCAACATGACAACCAAACATCTTCATATTCTCGGCATTTGCGGCACATTTATGGGCGGCATTGCGATGATCGCCAAGCAAATGGGCTATCATGTAACAGGCTCCGACACCAACGTTTACCCACCAATGAGCACTTTTCTTGCCGAAAATCACATTGAAATCATACCGCACTTTGAGCCAAGCCAGTTACAACCAGCGCCGGATATGGTCATCGTGGGAAATGCCATGAAACGCGGCAATCCCTGCATCGAGTATATGTTA  >  1/2055938‑2056420

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