Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 43,172 G→A T650T (ACC→ACT)  unknown ← Cytochrome c heme lyase subunit CcmF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*143,1720GA100.0% 158.1 / NA 55T650T (ACC→ACT) unknownCytochrome c heme lyase subunit CcmF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/30);  total (25/30)

CGGTTTGCCCTGTTTAAACAGCGTTTCATCATAGGTTTTATCCGATAATAAATCCTGCAACGTCTTCGGCGGAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACGGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGTTTGTAATACAAGCGCAACGCCCAGGAATTGTCGTCTAATTTTTCACCTAACGCCACATACAAGTCACGGCTGAAACCCCAATCGATCGCGGCTTCGGTCATCGGCATTTTGCTCACGGTGTATAAGCGTTTTT  >  1/42931‑43413
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
cGGTTTGCCCTGTTTAAACAGCGTTTCATCATAGGTTTTATCCGATAATAAATCCTGCAACGTCTTCGGCGGAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTa                                                                                                                                                                                                                                           <  1:314013/251‑1 (MQ=255)
        ccTGTTTAAACAGCGTTTCATCATAGGTTTTATCCGATAATAAATCCTGCAACGTCTTCGGCGGAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAAt                                                                                                                                                                                                                                   <  1:876478/251‑1 (MQ=255)
           gTTTAAACAGCGTTTCATCATAGGTTTTATCCGATAATAAATCCTGCAACGTCTTCGGCGGAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAAtttt                                                                                                                                                                                                                                >  1:650213/1‑251 (MQ=255)
                       tttCATCATAGGTTTTATCCGATAATAAATCCTGCAACGTCTTCGGCGGAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATa                                                                                                                                                                                                                    <  1:660687/251‑1 (MQ=255)
                            tcatAGGTTTTATCCGATAATAAATCCTGCAACGTCTTCGGCGGAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAAcgg                                                                                                                                                                                                                >  1:241191/1‑250 (MQ=255)
                                 ggTTTTATCCGATAATAAATCCTGCAACGTCTTCGGCGGAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATa                                                                                                                                                                                                                    <  1:753511/241‑1 (MQ=255)
                                    tttATCCGATAATAAATCCTGCAACGTCTTCGGCGGAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAAc                                                                                                                                                                                                       <  1:42484/251‑1 (MQ=255)
                                              aataaATCCTGCAACGTCTTCGGCGGAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAACTTTACTAAACCGATAACGGCGGTCGAACAGGCataata                                                                                                                                                                                             <  1:173405/251‑1 (MQ=255)
                                                    tCCTGCAACGTCTTCGGCGGAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCccgc                                                                                                                                                                                       <  1:160043/251‑1 (MQ=255)
                                                                  cggcggAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAaac                                                                                                                                                                         <  1:50458/251‑1 (MQ=255)
                                                                    gcggAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAaaca                                                                                                                                                                        <  1:372427/250‑1 (MQ=255)
                                                                     cggAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAaacaac                                                                                                                                                                      <  1:719331/251‑1 (MQ=255)
                                                                                     accaACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATcc                                                                                                                                                     <  1:768738/251‑1 (MQ=255)
                                                                                                              acgacaattcctgggcgttGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGTTTGTAATAc                                                                                                                              <  1:291991/232‑1 (MQ=255)
                                                                                                                tcatcaCCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑cc                                                                                                                                                           <  1:14025/218‑1 (MQ=255)
                                                                                                                tcatcaCCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCCCCCAGCACCATAAAAAACCCACCAAA‑CCATATCCAGCGGGTAAAGGGGTTGTAATACaaa                                                                                                                           >  1:323490/1‑249 (MQ=255)
                                                                                                                          gTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCCCCCAGCACCCTaaaa                                                                                                                                                                         >  1:264499/1‑194 (MQ=255)
                                                                                                                              aTTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGTTTGTAATACAAGCGCAACGCCCAGGaa                                                                                                            <  1:849201/251‑1 (MQ=255)
                                                                                                                                 gCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAaac                                                                                                                                                                         <  1:835650/188‑1 (MQ=255)
                                                                                                                                 gCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGTTTGTAATAc                                                                                                                              <  1:882070/230‑1 (MQ=255)
                                                                                                                                 gCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAAGCATAATATCCCGCCGAGCACCATAAAAAAACCACCCAG‑CCATATCCAGCGGGTAAAGGGGTTGTAAAAAAAACGCAAAGCCCAGGaaatg                                                                                                         >  1:1581/1‑248 (MQ=255)
                                                                                                                                      gTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCaa                                                                                                                                                               <  1:309463/193‑1 (MQ=255)
                                                                                                                                             cACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCCAGCACCATAAACAACCCACCAAA‑CCAAATCCCGCGGAAAAAAGGGTTGTAATAAAAACGCAACCCCCCGGaacttgctcttataaa                                                                                              >  1:616174/1‑236 (MQ=255)
                                                                                                                                                       gCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAAACCACCAAA‑CCCTATcc                                                                                                                                                     >  1:490189/1‑185 (MQ=255)
                                                                                                                                                          gTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCataata                                                                                                                                                                                             <  1:10224/143‑1 (MQ=255)
                                                                                                                                                            aaCCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAAAAGGCAAAATATCCCGCCCAGCACCATAAACAACCCACCAAG‑CCATATCCCGCGGGTAAAGGGGTTGTAATAAAAACGCAACGCCCCGGAATTGTCGTCTAATTTTTCACCTAACCCcaca                                                                              >  1:328250/1‑251 (MQ=255)
                                                                                                                                                                   aCGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGTTTGTAATACAAGCGCAACGCCCAGGAATTGTCGTCTAATTTTTCACCTAACGCCACATACAAGt                                                                       <  1:738503/251‑1 (MQ=255)
                                                                                                                                                                              aatTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCCCCCAGCACCATAAAAAACCCACCAAA‑CCCTATCCAGCGGATAAAAGGGTTGTAATACAAGCGCAACGCCCAGGAATTGTCGTCTAATTTTTCACCTAACGCCACATACAAATCACGGCTGaaa                                                            >  1:669357/1‑251 (MQ=255)
                                                                                                                                                                                          atatAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAaacaaacc                                                                                                                                                                    >  1:86997/1‑133 (MQ=255)
                                                                                                                                                                                          atatAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCCCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGTTTGTAATACAAACGCAACGCCCCGGAAATgt                                                                                                        >  1:295439/1‑195 (MQ=255)
                                                                                                                                                                                           tataATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAaac                                                                                                                                                                         >  1:110159/1‑130 (MQ=255)
                                                                                                                                                                                                  tCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTAc                                                                                                                                                                                                                                          <  1:326094/58‑1 (MQ=255)
                                                                                                                                                                                                            tATCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCataata                                                                                                                                                                                             >  1:221270/1‑93 (MQ=255)
                                                                                                                                                                                                             aTCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGTTTGTAATACAAGCGCAACGCCCAGGAATTGTCGTCTAATTTTTCACCTAACGCCACATACAAGTCACGGCTGAAAccc                                                         <  1:185144/223‑1 (MQ=255)
                                                                                                                                                                                                              tCAATACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCCCCCTAAACAACCCACCAAG‑CCATATCCAGCGGAAAAAGGGTTTGTAATACAAGCGCAACGCCCAGGGATTGTTGTCTAATTTTTCACCTAACGCCACATACAAGACACCGCCGAAACCCCAATCGATCGCGGCTTCCGTCATCGGCAt                            >  1:754697/1‑251 (MQ=255)
                                                                                                                                                                                                                 aTACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAaacaac                                                                                                                                                                      <  1:510803/111‑1 (MQ=255)
                                                                                                                                                                                                                 aTACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCCCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGGTAAAGGGGTTGTAATACAAGCGCAACGCCCAGGAATTGTCGTCTAATTTTTTACCTAACGCCACATACAAGTCACGGCTGAAACCCCAAACGATCGCGGCTTCGGGCATCGGCAtttt                         >  1:558968/1‑251 (MQ=255)
                                                                                                                                                                                                                 aTACTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCCCCGAGCACCATAAACAAACCACCAAG‑CCATATCCAGCGGGTAAAGGGGTTGTAATACAAGCGCAACGCCCCGGAATTGTCGCCTAATTTTTTACCTAACGCCACACACAAGTCACGGCTGAAACCCCAATCGATCGCGGCCTCGGTCATCGGCAtttt                         >  1:654459/1‑251 (MQ=255)
                                                                                                                                                                                                                    cTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGTTTGTAATACAAGCGCAACGCCCAGGaa                                                                                                            <  1:247742/165‑1 (MQ=255)
                                                                                                                                                                                                                    cTCGATGATTCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGTTTGTAATACAAGCGCAACGCCCAGGAATTGTCGTCTAATTTTTCACCTAACGCCACATACAAGTCACGGc                                                                 <  1:545244/208‑1 (MQ=255)
                                                                                                                                                                                                                            ttCTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAAAAGGCATAATATCCCGCCGAGCACCCTAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGGTTGTAATACAAGCGCAACGCCCAGGGATTGTCGTCTAATTTTTCACCTAACGCCACATACAAGTCACGGCTGAAACCCCAATCGATCGCGGCTTCGGTCATCGGCATTTTTCTCCCGgtgt              >  1:181101/1‑251 (MQ=255)
                                                                                                                                                                                                                             gcTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAAAAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGGTTGTAATACAAGCGCAACGCCCAGGAATTGTCGTCTAATTTTTCACCTAACGCCACATACAAGACACGGCTGAAACCCCAATCGATCGCGGCTTCGGTCATCGGCATTTTGCTCACGGTGta             >  1:230973/2‑251 (MQ=255)
                                                                                                                                                                                                                              cTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTTGAACAGGCCTAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCAAATCCACCGGATAAAGGGTTTGTAATACAAGCGCAACGCCCAGGACTTGTCGTCTAATTTTTTACCTAAAGCCACATAAAAGTCACGGCTGAAACCCCAATTCATCGCGGCTTCGGCCCTCCGCAATTTGCTCAccgcgcat            >  1:106235/1‑244 (MQ=255)
                                                                                                                                                                                                                              cTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTTGAAAAAGCAAAATATCCCCCCGAGAACCATAAAAAAACCCCCAAGCCCAAACCCA‑CGGGTAAAAGGGTTGGAATAAAAACGCAACCCCCAGGGATTTTCGGCTAATTTTTTACCTAAAGCCCCATACAAATCACGGCTGAAACCCCAAACGAACGCGGCCTCGGGCATCGGCCTTTTGCTCACGGGGtat            >  1:154269/1‑251 (MQ=255)
                                                                                                                                                                                                                              cTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCCTAAACAACCCACCAAA‑CCCTATCCCGCGGGTAAAGGGTTTGTAATACAAACGCAACGCCCAGGAATTGTCGTCTAATTTTTCACCTAACGCCACATACAAATCACGGCTGAAACCCCAATCGCTCTCGGCTTCCGTCATCGGCCTTTTTCTCACGGGGtat            >  1:829361/1‑251 (MQ=255)
                                                                                                                                                                                                                              cTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAggg                                                                                                                                        <  1:601431/127‑1 (MQ=255)
                                                                                                                                                                                                                              cTATAAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCAAAATATTCCGCCCAGCACCATAAACAACCCACCAAG‑CCAAATCCAGCGGATAAAGGGGTTGTAAAACAAGCGCAACCCCCCGGAATTGTCGGCCAATTTTTTACCTAAAGCCACCTACAAGACACCGCTGAAACCCCAAACGATCGCCGCTTCCGTTATCGGCATTTTGCTCACCGTGtat            >  1:598420/1‑251 (MQ=255)
                                                                                                                                                                                                                               tataAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAatat                                                                                                                                                                                            <  1:872268/75‑1 (MQ=255)
                                                                                                                                                                                                                               tataAGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAaac                                                                                                                                                                         <  1:713639/94‑1 (MQ=255)
                                                                                                                                                                                                                                   aGGGGTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGTTTGTAATACAAGCGCAACGCCCAGGAATTGTCGTCTAATTTTTCACCTAACGCCACATACAAGTCACGGCTGAAACCCCAATCGATCGCGGCTTCGGTCATCGGCATTTTGCTCACGGTGtata           <  1:200136/247‑1 (MQ=255)
                                                                                                                                                                                                                                     gggTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGGATAATATCCCGCCGAACACCATAAAAAAACCACCAAG‑CCATATCCAGCGGGTAAAAGGGTTGTAATAAAAGCGCAACGCCCAGGGATTGTCGTCTAATTTTTTTACTAACGCCACATACAAATCAccgc                                                                 >  1:253146/1‑188 (MQ=255)
                                                                                                                                                                                                                                     gggTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGTTTGTAATACAAGCGCAACGCCCCGGAATTGTCGTCTAATTTTTTACCTAACGCCACATACAAGTCACGGCTGAAACCCCAATCGATCGCGGCTTCGGTCATCGGCCTTTTGCTCACGGTGTATAAGCggt     >  1:12435/1‑249 (MQ=255)
                                                                                                                                                                                                                                     gggTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGTTTGTAATACAAGCGCAACGCCCAGGAATTGTCGTCTAATTTTTCACCTAACGCCACATACAAGTCACGGCTGAAACCCCAATCGATCGCGGCTTCGGTCATCGGCATTTTGCTCACGGTGTATAAGCg       <  1:685955/249‑1 (MQ=255)
                                                                                                                                                                                                                                       gTTACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAaac                                                                                                                                                                         <  1:177122/86‑1 (MQ=255)
                                                                                                                                                                                                                                        ttACTGCACAGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGCAGAACAGGCATAATATTCCCCCGAGCACCATAAACAAACCACCAAG‑CCATATCCAGCGGATAAAGGGCTTGGAAAACAAACGCAACGCCCAGGAAAATACGTCCAATTTTTTACCTAACGGCACATACAAGCCAAGGCTGAAACCCCAATCGGACGCGGCCTCGGTCCTCGGCATATTGCCCACCGTGTTAAAACGttttt  >  1:381727/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
CGGTTTGCCCTGTTTAAACAGCGTTTCATCATAGGTTTTATCCGATAATAAATCCTGCAACGTCTTCGGCGGAACCGGCTTACCCACCAACGCCGATTCCAATGCTTTCGGATCATCACCCAGTGCATTGCGTTGTAATTGCACCAGAAAGGCGGTAACCAAAACGAAGAAAATAATTAAAGGTAAATATAATTTCTTTTTCATTATCAATACTCGATGATTCTATAAGGGGTTACTGCACGGTTTCTTTACCTAAAATTTTACTAAACCGATAACGGCGGTCGAACAGGCATAATATCCCGCCGAGCACCATAAACAACCCACCAAG‑CCATATCCAGCGGATAAAGGGTTTGTAATACAAGCGCAACGCCCAGGAATTGTCGTCTAATTTTTCACCTAACGCCACATACAAGTCACGGCTGAAACCCCAATCGATCGCGGCTTCGGTCATCGGCATTTTGCTCACGGTGTATAAGCGTTTTT  >  1/42931‑43413

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